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Preface 


This book is based on a course on the differential geometry of curves and surfaces at 
Technische Universitat Berlin in the spring term of 2020. The 13 chapters roughly 
reflect the 13 weeks of that term. 

The pioneers of differential calculus like Newton, Bernoulli, and Euler immedi- 
ately applied their ideas to questions about curves and surfaces. In 1673, Newton 
defined the curvature « of a plane curve, and in 1691, Jacob Bernoulli characterized 
elastic plane curves (cf. [30], [25], or [3] for a historical overview), i.e. curves 
that minimize the bending energy f «? ds among all curves held fixed at their end 
points. In 1859, Kirchhoff showed that the tangent vector of an elastic curve follows 
the motion of the axis of a spinning top [19]. Even today, many applications of 
differential geometry of curves in other sciences (ranging from the coiling of DNA 
strands (cf. [37]) to the modeling of hair for computer-generated imagery (cf. [4])) 
are centered around elastic curves. Our approach to curve theory emphasizes its 
connections to the calculus of variations and we will explore elastic curves quite 
thoroughly. 

There is also a dynamic aspect of curve theory, where deformations of curves in 
time are studied. In 1906, Da Rios, a student of Levi-Civita, derived an evolution 
equation [33], the so-called filament flow, for space curves that models the motion 
of vortex filaments in a fluid (Sect. 5.3). In 1932, Levi-Civita wrote the equations 
satisfied by filaments that do not change shape under this flow [24]. In 1991, Langer 
and Perline showed that the possible shapes of such filaments are given by elastic 
curves [31], a fact that had escaped Levi Civita. Already in 1972, Hasimoto had 
shown that the filament flow is a so-called Soliton equation [15]. Even today this 
insight remains a source of ongoing inspiration for curve theory (see [10] for a 
survey). 

While minimizing the length of a curve results in straight line segments, 
minimizing the area f det of a surface with a given boundary curve leads to a 
rich class of surfaces, the so-called minimal surfaces (Sect. 12.4). Already in 1744, 
Euler proved that the catenoid minimizes area among all surfaces of revolution with 
prescribed boundary circles (cf. [13]). Minimizing area while fixing the enclosed 
volume leads to surfaces with constant mean curvature H (Sect. 12.5). For a surface, 
the analog of the bending energy f «? ds of a curve is the so-called Willmore 
functional f H? det (Sect. 13.1). In the context of surfaces, the analog of an elastic 
curve is a so-called Willmore surface, whose equation we derive in Sect. 13.2. It 
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was a major milestone in differential geometry when in 2012 Marques and Neves 
proved the so-called Willmore conjecture (cf. [26]), which states that for any torus 
in IR? the Willmore functional has to be at least 27°. 

A pervasive theme in differential geometry is the interplay between curvature 
and topology. In Sect. 3.4, we will show that the integral f « ds of the curvature 
of a closed plane curve y equals 27 times an integer, the so-called tangent winding 
number of y. In Sect. 3.6, we follow Whitney and Graustein who proved in 1937 that 
this integer characterizes the connected components of the space of all closed plane 
curves [45]. In the context of surfaces, the analog of this result is the Gauss-Bonnet 
theorem for closed surfaces [8], which we prove in Sect. 10.2. 

The only prerequisites for this book are the calculus of several variables 
including the transformation formula for integrals, the Picard-Lindelof theorem for 
ordinary differential equations, and Green’s theorem from vector calculus. Neither 
manifolds nor results from functional analysis are needed. Variational problems 
under constraints are accessible with these prerequisites because our definition of 
a critical point under constraints (Definition 2.19) is slightly stronger than the usual 
one. Similarly, our ability to discuss the genus of closed surfaces without diving into 
algebraic topology can be traced back to our definition of a compact domain with 
smooth boundary in IR? (Definition 6.1). Our definition is intuitive and equivalent to 
the standard one, but proving this equivalence would need serious additional work. 


Berlin, Germany Ulrich Pinkall 
January 2021 
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Curves in R” 


Differential Geometry studies smoothly curved shapes, called manifolds. One- 
dimensional shapes are called curves and two-dimensional shapes are called 
surfaces. In this chapter we look at curves in n-dimensional Euclidean space. The 
basic properties of curves in R” (length, tangent, bending energy) were explored 
right after the invention of calculus by Newton, Bernoulli and Euler. 


1.1 What is a Curve in R”? 


Since many interesting curves (for example a figure eight) have self-intersections, it 
is not a good idea to define a curve as a special kind of subset in R”. Intuitively, a 
curve is something that can be traced out (“parametrized”) as the path of a moving 
point (cf. Fig. 1.1). 


Definition 1.1 


A curve in R” is a smooth map y: [a,b] — R” such that its velocity vector 
y’ (x) never vanishes, i.e. 


y'(x) #0 
for all x € [a, b]. 


> Remark 1.2 If M C R” is an arbitrary subset, then a map f: M — R is called 
smooth (or C™) if there is an open set U C R” with M C U and an infinitely often 
differentiable map f : U — RÉ such that f = f |m (cf. Appendix A.1). Instead of 
a closed interval [a, b] one could also allow an open or semi-open interval (or even 
a finite union of intervals) as the domain of definition for a curve. The only problem 
that would arise is that then the integral of a smooth function would not always be 
defined. For all of our applications we can stick to closed intervals. 
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4 1 Curves in IR” 


a b R” 


Fig. 1.1 A curve can be described as the trajectory of a particle moving in space. The particles 
position at time x is given by y (x) 


Definition 1.3 


A curve y: [a,b] > R” is called closed if y can be extended to a smooth map 
y: R— R” with period b — a, which means 


yx + (b-a) = p(x) 
for all x € R. 


Example 1.4 


(i) The quarter circle is a curve: 


. i i 2 = t 
r: |- 4] R. y= ( — 


(ii) Another version of the quarter circle is also a curve: 


y: E 3] > R?, yœ) = ena? 


sin x 


(iii) The full circle 


y: [0,27] > RÊ, p(x) = ea 


sin x 


is a closed curve with period 27. It can be extended to 


y: R—> RŽ, y(x) = Pa r 


sin x 
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y'(@) 


Fig. 1.2 A circle (left), the Cartesian leaf (middle) and Neil’s parabola (right) 
(iv) The Helix is a curve: 
y: [a,b] > RŽ, yx) = | sinx 


(v) The Cartesian leaf (see Fig. 1.2) is a curve: 


3 
. 2 _ (x —4x 
vila b> Rx = (52) 


so that 


1, _ (3x7 -4 
ro=( 2x ). 


(vi) Neil’s parabola (see Fig. 1.2) is given by 


3 
y: [a,b] > R, y(t) = a f 


It is not a curve if 0 € [a, b], because at t = 0 


ro _ (0 
o= (g): 


For the purposes of geometry, the speed with which we run through a curve does 
not really matter, nor does the particular time interval [a, b] that we use for the 
parametrization. However, we will always assume that our curves are oriented, so 
we want to keep track of the direction in which we run through the curve. This 
means that we are only interested in properties of a curve that do not change under 
orientation-preserving reparametrization (see Fig. 1.3): 


6 1 Curves in R” 


Fig. 1.3 A reparametrization 
of a curve is given by a 
strictly increasing function 
with nowhere vanishing 
derivative which maps [c, d] 
onto [a, b] 


Definition 1.5 


Let y: [a,b] > R” and y: [c,d] — R” be two curves. Then + is called an 
orientation-preserving reparametrization of y if there is a bijective smooth 
map 9: [c, d] > [a, b] such that g'(x) > 0 for all x € [c,d] andy = y og. 


Example 1.6 
For the two curves y from Example 1.4 (i) and 7 from Example 1.4 (ii) we have y = y o ọ with 


Q: [š =| > [ vo 4] g(x) = cosx. 
> Remark 1.7 Orientation-preserving reparametrization is an equivalence relation 
on the set of curves in R”. Although we are ultimately only interested in properties 
shared by all curves in the same equivalence class, we will always work with a 
particular representative curve y. 


1.2 Length and Arclength 


The most simple numerical quantity that can be assigned to a curve as a whole is its 
length. 


Definition 1.8 


Let y : [a, b] — R” be a curve. Then the function 
v: [a,b] > R, t |y’@| 


is called the speed of y and 


b 
L(y) =f v 


is called the length of y. 


1.2 Length and Arclength 7 


The length of a curve does not change under reparametrization: 


Theorem 1.9 

Suppose y: [a,b] > R” and y: [c,d] — R” are two curves such that 
VY = y og for some diffeomorphism : [c,d] — [a, b]. Then y and y have 
the same length. 


Proof. By the substitution rule, we have 


d d b 
cy) = f Ww oo'l= f iv owe’ = f Iv’ = £0). 


a 


Example 1.10 


(i) For the half circle y : [0,7] > R?, 
cos x 
yœ) = ea 


we have |y’| = 1 and therefore L(y) = x. 
(ii) The line segment y: [a, b] > RŽ, 
x 
ya) = (5) 


has length L(y) = b—a. 


Definition 1.11 


A rigid motion of R” is a map g: R” — R” of the form 
g(y) = Ay +b 
where A € O(n) is an orthogonal matrix and b € R” is a vector. 
Rigid motions are those transformations of the ambient space IR” which preserve 
distances between points. Two shapes that differ only by a rigid motion are said to 


be congruent. Matching the physical intuition for curves as trajectories of a particle 
moving in space, the length of a curve is invariant under rigid motions: 


8 1 Curves in R” 


Theorem 1.12 
Let y: [a,b] > R” be a curve and g: R” — R” a rigid motion. Then 


L(goy)= L(y). 


Proof. For y = g o y we have y = Ay + b and y’ = Ay’. Therefore, 


b b 
c= f avi= f y= £0). 


Definition 1.13 


Let y : [a, b] — R” be a curve. Then the function 


t 
s: [a,b] > R, s(t) = L (Y lian) =l ly'| 


is called the arclength function (or arclength coordinate) of y. 


In most situations however, the arclength function s itself is less useful than its 
derivative, the speed s’ = v = |y’|. Using only v, not s, we can define the derivative 
with respect to arclength: 


Definition 1.14 


Let y : [a, b] > R” be a curve and v = |y’| its speed. Let g: [a, b] > RÝ bea 
smooth function. Then we define the derivative with respect to arclength of g 
as the function 


dg g 
ds” v 


and the integral over arclength of g as 


b b 
f gds el gv. 
a a 


> Remark 1.15 Once we have learned about 1-forms in Sect. 7.2, we will be able 
to interpret ds as a 1-form on [a, b] and dg as quotient of 1-forms, just as it had 
been the dream of Leibniz. For now, they are just R*-valued functions on [a, b]. 
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Theorem 1.16 

The arclength fucnction s: [a,b] —> R of a curve y: [a,b] > R” is an 
orientation-preserving diffeomorphism of the interval [a, b] onto the interval 
[0, L] where L = L(y). The reparametrization 


y: [0, L] > R”, y=yos! 


has unit speed, i.e. \y’| = 1. 


> Remark 1.17 It is common in the literature on curves to routinely assume that 
the curves under consideration have unit speed, usually expressed by saying that 
they are “parametrized by arclength”. We will not do this here, for the following 
reasons: 


(i) Making use of the derivative with respect to arclength defined in 1.14 gives 
us the same elegant formulas as they arise in the context of unit speed curves, 
without actually changing the parametrization. 

(ii) When dealing with one-parameter families t + y, of curves of varying length, 
one cannot assume that all curves y; are parametrized by unit speed. Therefore, 
in this situation one has to resort anyway to formulas that remain valid for 
arbitrary curves. 

(iii) In the context of surfaces, there is no obvious analog for the unit speed 
parametrization of a curve. Therefore, habitual reliance on unit speed 
parametrizations makes the theory of surfaces look more different from the 
theory of curves than it actually is. 


1.3 Unit Tangent and Bending Energy 


Definition 1.18 


For a curve y : [a, b] > R”, the normalized velocity vector field 


dy y 
ds |y’| 


T: [a,b] > 8), T = 


is called the unit tangent field of y. 
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Next to the length, the most important numerical quantity that can be assigned to a 
curve as a whole is its bending energy: 


Definition 1.19 


Let T be the unit tangent field of a curve y: [a, b] —> R”. Then 


1 f?]dT dT 
B(y) =- —,—)d 
(y) 5S (F aig 


is called the bending energy of y. 


The bending energy is invariant under orientation-preserving reparametrization. 
The name comes from the following physical picture: 

Consider a rod manufactured out of some elastic material in the shape of a thin 
cylinder of length L and radius €. Then we bend the cylinder into the shape of 
a curve y of length L. While doing this, we make sure that we do not force any 
twisting on the cylinder, for example we place the cylinder in a hollow tube with 
shape y, leaving it free to untwist itself within the tube (see Fig. 1.4). Then, in the 
limit of € — 0, the energy needed to bring the initially straight rod into its new 
shape will be proportional to B(y). 

In later sections we will find out what curves we obtain if we hold a curve fixed 
near its end points but otherwise let it minimize bending energy (cf. Fig. 2.3). We 
also will find a way to deal with twisting. 


Fig. 1.4 A rod is bent into the shape of a curve. Then it is fixed in its position by a porcelain case 
within which it can untwist while staying in shape 
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Variations of Curves 


Many important special curves y arise by minimizing a certain variational energy 
E(y). For example, E(y) could be a linear combination of length and bending 
energy, in which case the curve is called an elastic curve. We are not only interested 
in minima but also in unstable energetic equilibria, possibly under constraints 
like fixing the curve near its end points. In this chapter we develop the basics of 
Variational Calculus. In particular, this allows us to explore elastic curves. Beyond 
straight lines and circles, these are the most important special curves in R”. 


2.1 One-Parameter Families of Curves 


On many occasions we will have to deal not only with individual curves y : [a, b] > 
R” but with whole one-parameter families t > y; of curves. 


Definition 2.1 


Let g;: [a,b] > R* be a smooth map, defined for each t € [fo, t1]. Then the 
one-parameter family of maps [fo, t1] > t + gr is called smooth if the map 


[a, b] x [tH] > RY, Œ, > g(x) 
is smooth (as always, in the sense of Remark 1.2). 
Given a smooth one-parameter family 
t> (gr: [a,b] > B®), t € [rot] 
of maps, also 
teg 
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is a smooth one-parameter family of maps g/: [a,b] > RŽ. The same holds for 
t > g; where g;: [a, b] > Ré is defined as 


: _ a 
B(x) := Be | ej 8x (x). 


The dot and prime derivatives are just partial derivatives, so they commute by 
Schwarz’s theorem: 


Theorem 2.2 
For a smooth one-parameter family of maps t +> gr, where g;: [a,b] > R* 
we have 


(2Y = GY. 


In our context, one-parameter families of maps will mainly arise as variations of 
a single map g: [a,b] > RE: 


Definition 2.3 


A smooth one-parameter family tf + g; of maps from M to R* is called a 
variation of a smooth map g: M —> R* if to < 0 < t and go = g. In 
this context, we will also use the notation 


& := ŝo. 


Our main interest is in variations of curves y: [a, b] —> R” (and the associated 
variations of derived quantities like the unit tangent or the length): 


Definition 2.4 


For a variation t > y of a curve y: [a, b] > R” the map 
Y:=y: [a,b] > R” 
is called its variational vector field. 


Suppose we have a smooth one-parameter-family t > yr of curves y: [a, b] > 
R”, meaning that y/(x) # 0 for all x € [a,b] and all t € [t,t]. Then we 
can think of this family Gust for the purpose of intuition, no need for further 
formal definitions) as a smooth map from [fo, tı] into in the space M of all curves 
y: [a,b] > R”. The vector ý% € C™([a, b], R”) can then be thought of as the 
“velocity vector” of that map at time f (see Fig. 2.1). 
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Fig. 2.1 A variation of a 
curve y can be interpreted as 
a map into the space M of all 
curves y : [a, b] > R” 


C™([a, b], R”) 


> Remark 2.5 Throughout this whole book we will treat C°° (a, b], R*) (and its 
analog in the context of surfaces) only as a vector space, based on notions from 
Linear Algebra. So, for example, we will indeed use the Euclidean inner product 


b 
lg, hi) = f te, hi 


but we will never put any topology on C% ([a, b], R*). This means that you will 
get confused if you try to interpret what we say based on notions from Functional 
Analysis. These notions have important applications in Differential Geometry, but 
they are not used at all in this book. 


2.2 Variation of Length and Bending Energy 


Given a variation t + y, of a curve y: [a,b] — R”, we want to determine 
£ | 10 EYD) and 4 | 1=0 B(y,). We first compute the time derivative of the inte- 
grands of these integrals: 


Theorem 2.6 

Let t +> y be a variation with variational vector field Y : [a,b] > R” ofa 
curve y: [a,b] — R” with speed v = ds and unit tangent field T. Then the 
variation of ds is given by 


o dY 
ds = == I eS, 
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Proof. Differentiating the equation v; = |y/| with respect to ¢ at t = 0 we obtain 


Wid / y’ Y 
pa Pari =(2 y) = (r)a. 
v ds 


oO 


Before we proceed to compute the rate of change for the bending energy integrand, 
note that (unlike the situation for partial derivatives), for a one-parameter family 
t œ> yı the derivative with respect to t does not commute with the derivative with 
respect to arclength: 


Theorem 2.7 

Let t œ> y be a variation with variational vector field Y : [a, b] > R” of a 
curve y : [a, b] > R” with speed v = ds. Then for any one-parameter family 
t +> g; of functions gr: [a, b] > R* with go =: g we have 


dg\ dg dY A dg 
ds} ds ds’ ds" 


Proof. By Theorem 2.6, 


dg = (£ eY (aT), _ dg _faY 5 \ ae 
ds v v we © as ds’ | ds’ 


Theorem 2.8 

Given a variation t œ> y; with variational vector field Y : [a,b] — R” ofa 
curve y: [a,b] > R” with speed v = ds, the corresponding variation of the 
bending energy density is 


1 /dT aT\ g = d&Y dT 3/d¥ dT dT ; 
aiam Tl Jom Wae a 2a ei 


Proof. Applying Theorem 2.7 to g = y we obtain 


. dY (5 ) 
T = — — | —, T)T. 
ds 
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Using this, Theorem 2.6, the fact that (T, T) = 1 implies (Z, T) = 0 and 
Theorem 2.7 with g = T we obtain 


1 /dT dT 4 dT\ dT 1/dT dT\/[dyY 
a Gere aera ds = ; ds + , ,T ds 
Ge a) ) (F) a) Ate alte 


dT |dY dT dT 1/dT dT\/dY 
( r) Jasta }( .t)as 


ds ds’ ds’ ds 2\ds’ ds ds 


dT dT 1/dT dT\|/dY 
= ; ds ; ,T }ds 
(g 3 Age aha 


dY dT dT 1/dT dT\/dyY 
LT ’ ds ; ,T ds 
(z ka =} Ate malts 


d&Y dT 3 /dT dT\|/dY 
Jas al \{ .T)as 


ds?’ ds ds’ ds | \ds 


Oo 


The proof of the following theorem is based on applying integration by parts 
repeatedly. 


Theorem 2.9 

Given a variation t œ> y; with variational vector field Y : [a,b] > R” of 
a curve y: [a,b] — R”, the corresponding variations of the length and 
bending energy are 


b 
ons hr = f (r, Tas 
t=0 a S 


a a ye aT dT\,, 
dtl T Was’ ds ds? | 2\ds’ ds 


dt 


b 


a 


+f pee dT &@T fee dT dT\dT i 
L ere ds’ ds? 2\ds’ ds} as}) 
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Proof. By Theorem 2.8, 


d dT dT i 
a= S GETA) 
[( &@Y dT\ 3/dY dT dT 
= , wi ; ds 
d ds?’ ds 2\ds ds’ ds 
-f d |dY dT dY @T 3d YT) dT dT 
~ Ja \ds\ds’ ds ds’ ds? 2ds \ ds’ ds 
3/,, aT dT dT she dT dT i 
teaho a Um ao 
b d ([dY dT @T\ 3 dT dT 
=] , Y, (Y, T) , ds 
a ds ds ds ds? 2 ds ds 
+f (t aT 4 dT dT pie dT dT a 
’ j S 
z ds? ds’ ds? ds 
dY dT dT dT dT 4 
= , Y, +5 T 
ds’ ds ds? 
+f . dT, dT T PE. dT dT j 
p ’ ds? ds’ ds? ds’ d k 


(mi 
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2.3 Critical Points of Length and Bending Energy 
Variations of curves (as defined in Definition 2.3) are needed in order to define 
and determine those curves that represent equilibria of geometrically interesting 


variational functionals. Functionals are certain real-valued functions on the space 
M of all curves y : [a, b] — R”, that was already introduced in Sect. 2.1. 


Definition 2.10 


Suppose we have a way to assign to each curve y: [a, b] —> R” a real number 
E(y). Then € is called a smooth functional if for every smooth one-parameter 
family 

trey, te [t,t] 
of curves y: [a, b] > R” the function 


is smooth. 


In many circumstances, we want to consider only variations of y : [a, b] > R” 
that keep the curve fixed near the boundary of the interval [a, b] (see Fig. 2.2). 


Definition 2.11 


Let y : [a, b] — R” a curve. Then a variation 
trey, tell 


of y is said to have support in the interior of [a, b] if there is € > O such that 
for all x € [a,a + €] U [b — €, b] we have 


y(x) = y(x) forall t € [t,t]. 


Fig. 2.2 A variation y; of a curve y with support in the interior 
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Now we can make precise what we meant by an equilibrium of a variational 
energy: 


Definition 2.12 


Let € be a smooth functional on the space of curves y: [a,b] — R”. Then a 
curve y: [a,b] > R” is called a critical point of £ if for all variations t œ> y; 
of y with support in the interior of [a, b] we have 


d 
— E =0. 
dt |o (vı) 


We denote the space 
fY :[a, b] > R” smooth | Y lfa a+5]U[b—5,b] = 9 for some ô > 0} 


of all functions Y: [a,b] — R” with support in the interior of [a, b] (Defini- 
tion A.4) by C8 (a, b), R”). 


Theorem 2.13 

For every vector field Y : [a,b] — R” along a curve y: [a,b] > R” there 
is a variation t œ> y; with variational vector field Y. If Y has support in the 
interior of |a, b], then also the variation t +> y; can be chosen in such a way 
that it has support in the interior of [a, b]. 


Proof. The proof of Theorem 2.13 is left as an exercise. o 


> Remark 2.14 In the case of the length functional, instead of using variations with 
support in the interior we could have used variations that fix both end points. For 
other variational problems (that involve higher derivative), additional derivatives 
(not only the function value) of y would have to be clamped to fixed values at the 
end points. On the other hand, variations with support in the interior will work all 
the time, with equivalent results. 


Theorem 2.15 (Fundamental Lemma of the Calculus of Variations) 
On the vector space C™ ({a, b], R”) equipped with the inner product 


b 
(f, 8) =| (f, 8) 


(continued) 
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Theorem 2.15 (continued) 
only the zero vector is in the orthogonal complement of Cy° ((a, b), R”): 


CF ((a, b), R")~ = {0}. 


Proof. Suppose that f € C™ ([a, b], R”) would be non-zero but in the orthogonal 
complement Ce ((a, b), R”)+. Then there would be xọ € [a, b] such that f (xo) Æ 
0. Choose ô > O such that [xg — ô, x9 + 6] C (a, b) and (f(x), f (xo)) > O for all 
x € [xo — ô, xo + ô]. Construct a smooth bump function (cf. Appendix A.2) 


g € Cy (xo — ô, xo + 8), R”) C CH? Cla, b), R”) 


such that g > O and g(xo) = 1. Then (f,g) # 0, which implies f ¢ 
Cx (a, b), R”)+, a contradiction. o 


Now we are in the position to determine the critical points of the length functional: 


Theorem 2.16 

A curve y: [a,b] > R” is a critical point of the length functional £ if and 
only if its unit tangent field T : [a, b] > R” is constant, i.e. if y parametrizes 
a straight line segment. 


Proof. By Theorem 2.9 and 2.15, y is a critical point of £ if and only if for all 


Y € C9 ((a, b), R”) we have 
(ET) 
Y, —))=0. 
ds 


By Theorem 2.15 this is the case if and only if af = 0. o 


Definition 2.17 


A curve y: [a,b] > R” is called a free elastic curve if it is a critical point of 
the bending energy functional B. 


22 2 Variations of Curves 


An almost identical proof as the one of Theorem 2.16 gives us 


Theorem 2.18 
A curve y: [a,b] > R” is a free elastic curve if and only if its unit tangent 
field T: [a,b] > R” satisfies 


dT dT at 3 JdT dT\ dT 
ae seh : =0 
ds? ds ds 2\ds ds] ds 
or equivalently 
ae d'y d*y\d 3 [dy d-y\ d? 
y +3 x Y Y a K Y Y =A 
ds? ds?’ ds? | ds 2 \ds?’ ds? | ds? 


Solving the fourth order differential equation for y that appears in Theorem 2.18 
with suitable initial values will give us unit speed parametrizations of free elastic 
curves. In the next chapter we will explore in more detail the geometric conse- 
quences of this differential equation. 


2.4 Constrained Variation 


In the context of many variational problems that arise in applications, general 
variations might violate some constraints that are imposed by the nature of the 
problem at hand. For example, thin elastic wires (for most practical purposes) have a 
fixed length. This means that here we should minimize bending energy only among 
those curves (held fixed near their boundary) that have a prescribed length. 

This kind of problem is known under the name of optimization under con- 
straints. Here we will work with a definition of a critical point under constraints 
that is slightly stronger than the standard one. The usual definition would replace 
the condition 4 | 1=0 E=0 by the requirement that E (y+) is independent of t. 


Definition 2.19 


Let £, € be two smooth functionals on the space of all curves 7: [a,b] > R”. 
Then a curve y : [a, b] > R” is called a critical point of € under the constraint 
of fixed £, if for all variations t > yr of y with support in the interior of [a, b] 


d x 
—| E=0 
dt |;=0 
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implies 


d 
dt 


t=0 


Both for the length functional € = C and for the bending energy € = B we 
know (Theorems 2.8 and 2.9) how to express the infinitesimal variation of E that 
corresponds to a variational vector field Y : [a, b] > IR” with support in the interior 
of [a, b] as an integral 


d 


b 
E= Y,G 
dt t=0 i r) 


for some smooth map G; : [a, b] > R”. If a formula like the one above holds, Gy 
is called the gradient of the energy € at y. 


Theorem 2.20 

Let E€, E be two smooth functionals on the space of all curves y: [a, b] > R”. 
Suppose we have a way to associate to each curve y: [a,b] — R” smooth 
maps 


Gre Gy: [a,b] > R” 


such that for all variations t +> y; of y with support in the interior of [a, b] 
we have 


d b 

— e= y,G 
de |p fe y) 
d z eE 
— e= AGNE 
dt |o fe y) 


Then y is a critical point of E under the constraint of fixed E if and only if 
there is a constant à € R such that 


à is called a Lagrange multiplier for the constraint of fixed È. 


Proof. We apply Theorem 2.21 below to the case where H = C®™ ([a, b], R”), 
V = Ch ((a, b), R”) and U = RG,. Then y is a critical point of € under the 
constraint of fixed € if and only if G; is orthogonal to all Y that are simultaneously 
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in V and orthogonal to U, i.e 
G, € (V NUH =U =RG,. 
E 


The theorem below is pure linear algebra, no Functional Analysis is involved. The 
formulation is such that it can also be applied to a situation where there are constraint 
functionals €),..., Ex instead of a single functional E. 


Theorem 2.21 

Let H be a (possibly infinite dimensional) vector space with inner product 
(.,.). Let V C H be a subspace such that Vt = {0} and U C H finite 
dimensional. Then 


(UtnV)t =U. 


Proof. As for all x € UŁAN V it holds that (u, x) = 0 for all u € U, the inclusion 
U c (Ut N V)+ is immediate. In order to show that also (Ut N V)+ c U we 
choose an orthonormal basis {u1, .. . , un} of U and define the map 


n 
P: H> U, xP S (x, uiui. 


i=l 


It is not hard to check that P defines an orthogonal projection of H onto U, i.e. 
P? = P, P* = P andimP = U. Now for u € U and h € H it holds 


(u, h) = (P (u), h) = (u, P(h)). 


Therefore we have U N P(V)+ c V+ = {0}, hence P(V) = U. So there are 
U1, ..., Un E€ V such that P(v;) = u;i. We now define the map 


n 


Q: H >V, xme XOG, vj) uj 


i=l 


which is symmetric (i.e. Q* = Q) and satisfies im Q C V and P o Qly = idy. 
Therefore, for x € (Ut N y)+ and v € V: 


(x — P o Q(x), v) = (x, v — Qo P(v)) =0, 
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1 
W 


Fig. 2.3 Elastic curves are everywhere 


since v — Q o P(v) € UL NV. Thus x — P o Q(x) € V+ = {0} and therefore 
x=PoQ(x)eU. o 


Definition 2.22 
A curve y: [a,b] — R” is called a torsion-free elastic curve if it is a critical 


point of bending energy under the constraint of fixed length (Fig. 2.3). 


Theorems 2.8, 2.9 and 2.20 together allow us to characterize torsion-free elastic 
curves by a differential equation: 


Theorem 2.23 
A curve y: [a,b] > R” is a torsion-free elastic curve if and only if there is a 
constant à € R such that its unit tangent field satisfies 


T dT dT 3 {dT dT\ dT ay 
ds3 ds’ ds? 2\ds’ ds| ds ds 
or, equivalently, 
ae, ay @y\ dy 3f}d*y dy\ dy ay = 
ds4 ds?’ ds? | ds  2\ ds?’ ds? | ds? ds ` 


The constant x is called the tension of y. 
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2.5  Torsion-Free Elastic Curves and the Pendulum Equation 


By Theorem 1.16 every curve in R” admits a reparametrization y: [0, L] > R” 
with unit speed. Then for any function g: [0, L] —> RE the derivative with respect 
to arclength is just the ordinary derivative: 


Theorem 2.24 

A curve y: [0, L] > R” with unit speed is torsion-free elastic with tension 
à if and only if its unit tangent field T : [a, b] > S”! solves the equation of 
motion 


TO TT —a= (a, T\P 


of a spherical pendulum with unit mass and some gravity vector a € R” and 
à equals the total energy of the pendulum: 


A= T TA =a T 


NI = 


Proof. Let y: [0, L] — R” be a torsion-free elastic curve with tension à and with 
unit speed. Then, by Theorem 2.23 


3 3 Í 
0= T” +3(T', TOT + zT T5T' -AT = (r+ zT T^T -a7') ; 
i.e. if there is a constant vector a € R” such that 
A 3 / i. 
T +P. TT ar =a. 


Looking at the component orthogonal to T on both sides of this equation gives us 
the first of the two equations that we want to prove. Taking the scalar product with 
T and using 


1 
O= S(T, Pi ate Pag Te 


we obtain the second equation. Conversely, if T: [0, L] —> S”! solves the 
pendulum equation 


T" —(T",T)T =a — (a, T)T, 
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then it is easy to verify that the total energy à defined by 

1 / / 

X= ~(T’,T’) — (a, T) 

2 

is constant and 
" 3 / / 
oO 


Figure 2.4 shows planar torsion-free elastic curves that lie in a plane. They arise 
from pendulum motion on a circle, whereas a pendulum motion on S? gives a 
torsion-free elastic curve in R3 as seen in Fig. 2.5. 
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Fig. 2.4 Trajectories of a pendulum (drawn in blue color). Below each of these trajectories the 
corresponding torsion-free elastic curve is shown 
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Fig. 2.5 Trajectory of a pendulum on S? (drawn in blue color), together with the corresponding 
torsion-free elastic curve in R? 
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Curves in R2 


Curves in the plane IR? are special in several respects: For a closed plane curve y 
an enclosed area A(y) can be defined, providing another geometric functional in 
addition to length and bending energy. Unlike the situation in higher dimensions, 
the geometry of an arbitrary unit speed plane curve y : [a, b] > IR? is captured in a 
smooth real-valued curvature function x : [a,b] > R. We prove our first theorem in 
Global Differential Geometry: The integral of the curvature of a closed plane curve 
is 277n where n is an integer, called the tangent winding number of y. Two closed 
plane curves can be smoothly deformed into each other if and only if they have the 
same tangent winding number. 


3.1 Plane Curves 


The case of curves y: [a,b] —> IR? is special because R? comes with a 
distinguished linear map J: IR? — R?, the 90°-rotation in the counterclockwise 
(positive) direction: 


woe -CO 


Here are some properties of J that are easy to check: We have J? = —I and J is 
orthogonal as well as skew-adjoint, i.e. for all vectors X, Y € RÊ? we have 
(JX, JY) = (X,Y) 
(JX,Y) = —-(X, JY). 
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Furthermore, the determinant function det on IR? can be expressed in terms of J and 
the scalar product: 


(JX, Y) = det(X, Y). 


If y: [a,b] > R? is a curve and T: [a,b] > R? is its unit tangent, then 4° is 


ds 
orthogonal to T and therefore proportional to JT: 


Definition 3.1 


Let y : [a, b]  R? be a curve and T: [a, b] — R? its unit tangent. Then the 
unique function « : [a, b] > R such that 


dT a 
eee: 
ds 


is called the curvature of y. 


More explicitly, 


dT 1 1/1 i det(y’, y” 
K=(JT, = Jy’, y Se ) 
ds v v wv ly’ 


We convince ourselves that the curvature is independent of the parametrization: Let 
y = y og bea reparametrization of a plane curve y. Then 


Y=g'-y'og 
and, since gy’ > 0, |7’| = |v’ o o|. Hence 
y=" -yloet+¢y'-y" 0g) 
and therefore 


det(y’, 7") — @ PP - dety’, y") _ 
In’) Ply’? 


k= 


The curvature « of a straight line segment vanishes since we can parametrize the 
segment going from a point p € R? to a point q € RÊ by 


x—a 
y: [a,b] > R?, x> p+>—(G-p) 
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and therefore y” = 0, so that « = 0. A circular arc 


y: [a,b] > RÊ, x > Ea 


r sinx 


of radius r has constant curvature k = i. If we restrict attention to unit speed 


curves y : [0, L] —> IR’, the curvature function «: [0, L] —> R determines y up to 
orientation-preserving congruence: 


Theorem 3.2 (Fundamental Theorem of Plane Curves) 


(i) For every smooth function x: [0, L] —> R there is a unit speed curve 
y : [0, L] > R? with curvature x. 

(ii) If y,y: [0, L] = R? are unit speed curves with the same curvature 
function x, then there is an orthogonal (2 x 2)-matrix A with det A = 1 
and a vector b € R? such that 


ý = Ay +b. 


Proof. For (ii), denote by T, T the unit tangent fields of y and y and take for A the 
orthogonal (2 x 2)-matrix A with determinant one for which AT (0) = T (0). Then 
both T and 

Î := AT 


solve the linear initial value problem 


Y(0) = T(0) 
Y’=kJY 


and therefore, by the uniqueness part of the Picard-Lindeléf theorem, we must have 
T = T. Then 


(p — Ayy =I (Tf —T) =0, 
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which proves (ii). For (i), define aw: [0, L] ~ Rand T, y : [0, L] > R2 by 


a(x) := f K 
0 
c= i) 
sina 
y(x) := f T; 
0 


Then |T| = 1 and y’ = T, so y is a curve and T is its unit tangent field. 
Furthermore, T’ = « JT and therefore y has curvature x. o 
3.2 Area of a Plane Curve 
Let y : [a, b] > R? be a curve such that det(y, y’) > 0 and the map 

f: 0,1] x [a,b] > R?, f(t, x) = ty (x) 


is a bijective map onto a subset S C RŽ. Then the derivative f'(t, x) at the point 
(t,x) € (0, 1) x [a, b] satisfies 


det f’(t, x) = t det(y (x), y’(x)) > 0 


and using the transformation formula of integrals it is not difficult to show that the 
area of S is given by 


b 1 1 b 
area(s) = [1 = f i=/ f det f' = =f det(y, y^). 
s f(0,11x[a,b]) a Jo 2 Ja 


For the curve y shown on the left of Fig. 3.1, the above formula correctly yields 
the area enclosed by y and the line segments from the origin to y (a) and y (b). It 


Fig. 3.1 For the curve on the 
right, the position vector from 
the origin to y (x) covers 
some areas multiple times 


a 


b 
| 


A, 
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therefore seems reasonable to use this formula in order to define an area for arbitrary 
curves y: [a,b] > R?: 


Definition 3.3 


The sector area of a curve y : [a, b] — R? is defined as 


1 b 
Ay) = 5S det(y, y’). 


The curve on the right of Fig. 3.1 illustrates the consequences of this definition. 
There the position vector from the origin to y(x) covers some areas multiple 
times. However, for some of these times (where y, as seen from the origin, moves 
clockwise) the contribution to the covered area, as it is computed by the above 
formula, is negative. 

The sector area A(y) depends on the origin in IR?, which means that it changes 
if we apply a translation p +> p — v to y. Therefore, at first sight the sector area 
does not look like a good geometric invariant for curves. However, this dependence 
disappears as soon as we restrict attention to closed curves, or consider differences 
between the sector areas of curves that share the same endpoints (see Fig. 3.2): 

Let v € R° be a vector and y : [a, b] — R? a curve. Then we define a modified 
sector area Ay(y) as the sector area of the curve y translated by the vector v: 


Ay(y) := A(y + y). 


Fig. 3.2 Independent of the choice of origin in IR”, the sector area of the curve on the left of the 
above picture equals the area of the blue region minus the area of the orange region. Similarly, 
the difference of the sector areas of the two curves on the right equals the area of the blue region 
between them 
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Theorem 3.4 
Let y: [a,b] > R? be a closed curve and yı, ¥2: la, b] > R2 two curves 
with yı (a) = y2(a) and yı (b) = y2(b). Then, for any vector v € R2 we have 


Av(y) = A(y) 
Ay(2) — Avi) = AQ2) — A(1)- 


Proof. Because y is closed, we have 


1 fP 1 fP 
Avy) Ay = 5 f det ery) dy’) 


a 


1 b 
=- | det(v,y’ 
Al et(v, y) 


1 b 
=- det(v, y) 
>| et(v, y) 


: t( y JA 
= det(v, 
2 a 
= 0. 
By the same arguments we obtain 
(Ay(2) = Av(11)) = (A2) = A(1)) 
= 1 f ae yok f ae y 
= et(v, et(v, 
7 7 y2 7 : yı 


= det(v, » — vw)? 
=0. 


oO 


In particular, we expect that for variations with support in the interior of [a, b] of 
a curve y : [a,b] —> R?, the corresponding variation of sector area is independent 
of the choice of origin: 
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Theorem 3.5 
Let t œ> y, be a variation with support in the interior of [a, b] of a curve 
y: [a, b] > R’. Then 


b 
E| Amw=- / (Y, Ty’). 


p=) 


Proof. Since Y vanishes at the endpoints, we have 


d 


dt 


Aly) 


t=0 


1 fÈ no Lpf TE 
3S decry +5 f det(y, (y’)’) 


1 fÈ 1 f? 
=; dey) +5 f det(y, Y’) 
Dafa 2 Ja 
1 fÈ 1 f? 
=- | de(Y,y)—- | det(y’,Y 
S et(Y, y^) al et(y', Y) 


b 
= f det(Y, y’) 
a 


b 
= -f (Y, Jy’). 


As a consequence, the sector area functional by itself does not have any critical 
points. On the other hand, minimizing length among all curves with the same 
endpoints and the same sector area is possible: 


oO 


Theorem 3.6 

A curve y: [a,b] > R? is a critical point of length under the constraint of 
fixed sector area if and only if its curvature x is constant, i.e. if and only if its 
image lies on a circle or a straight line. 


Proof. By Theorems 2.9 and 2.20 y is a critical point of length under the constraint 
of fixed sector area if and only if there is a constant A € R such that 


M-Jy') = -T! =-KJy’. 
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3.3 Planar Elastic Curves 
For a unit speed curve y : [0, L] > IR? with unit tangent T and curvature « we have 


T'=KJT 
T" = —=k°T +K'JT 


T” = —3kk'T + (k" — )IT 


The bending energy of a plane curve is also called its total squared curvature. This 
is because for a unit speed plane curve y as above we have 


1 f? oe me 
Biy)=-= T',T'\d =; ds. 
o) sf jds=3 |] ds 
By Theorem 2.23, y is an elastic curve with tension À if and only if 
m / ” 3 / f f / 
O=T  +3(T,T ee ane — AT 
ie 
= (KS + Aa) IT 


which means 


3 
hid K 
K ETAREN, 


This differential equation can be interpreted as the equation of motion 
ƏV 
k” + —(k)=0 
ðk 
for a particle with unit mass moving on the real line subject to the potential energy 
V(k) = Ta + AK? 
As expected (and as is easy to verify by taking the derivative) the total energy 
1 1\2 
BE: 5) + Vik) 


is constant. In particular, we see that along for each solution the potential energy is 
bounded from above by E. In Fig. 3.3 we see examples that should be compared to 
the shapes of the corresponding curves that were shown in Sect. 2.5. 
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Fig. 3.3 The potential wells for various values of À. For a solution « of the equation of motion, 
V(x) — E is always non-positive. The values of « that satisfy this condition are indicated in blue 


Fig. 3.4 A curve which is a critical points of the total squared curvature with constrained length 
and the sector area 


If we look for critical points of the total squared curvature while constraining not 
only the length but also the sector area, by Theorem 3.6 we arrive at the differential 
equation for x: 


K? 
E +e + = 0. 


The closed curve in Fig. 3.4 is such a critical point. 
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3.4 Tangent Winding Number 


Definition 3.7 


For a curve y: [a, b] > R? with curvature « the integral 


b 
f k ds 
a 


is called the total curvature of y. 


In this section we will prove that for a closed curve in RÊ the total curvature is 
an integer multiple of 27: 


Theorem 3.8 
Ify: [a,b] > RŽ is a closed curve with curvature x, then there is an integer 


n € Z such that 
b 
f Kds = 2n: 
a 


n is called the tangent winding number of y. 


Proof. Define a: [a, b] > R by 


a(x) := œo + f k ds 
a 
where œo is chosen in such a way that 
T (a) = (cos ao, sin œọ). 
As in the proof of Theorem 3.2, we conclude 
T = (cosg, sing). 
Since y is closed, we have T (b) = T (a), which means 


(cosa(b), sina(b)) = (cos a (a), sing (a)). 
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Fig. 3.5 The path on S! of the unit tangent can be visualized more clearly if it is drawn slightly 
outside of the unit circle 


90800 


Fig. 3.6 A list of representatives for every homotopy class of plane curves 


Therefore, there is an integer n € Z such that 


b 
J k ds = a(b) — a (a) = 27n. 
o 


As is clear from the above proof, the tangent winding number counts how often the 
unit tangent T (x) turns around the unit circle S ' as x runs from a to b (see Fig. 3.5). 
Figure 3.6 shows that all integers n € Z arise as the tangent winding number of 
some curve in R?. 


3.5 Regular Homotopy 


The following two sections will deal with the question: “Given two curves y, y in 
R”, is it always possible to smoothly deform y into y through intermediate curves?” 
For convenience, we assume that y and y have the same parameter interval. 


Definition 3.9 


A regular homotopy between two curves y, 7: [a, b] > R” is a one-parameter 
family t > y, of curves w: [a,b] —> R”, defined for t € [0,1], such that 
yo = y and yı = F. If there exists such a regular homotopy, y and y are called 
regularly homotopic. 
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Regular homotopy is an equivalence relation on the set of curves y: [a,b] > 
R”: Reflexivity and symmetry are easy and for transitivity we make use (see 
Appendix A.2) of a smooth function A: [0, 1] — [0, 1] such that 


iz 0, forx € [0,€] 
1, forx € [l-e, 1]. 
If now t t+» y is a regular homotopy between y and y and t > j a regular 
homotopy between y and y then 


Va(2t), fort € |o, l 


2 
Ynar-1), fort € (3. 1 


is a regular homotopy from y to 7. One can think of regular homotopies as smooth 
paths in the space of all curves y: [a,b] — R”, the equivalence classes under 
regular homotopy are the path-connected components of this space. Indeed, this 
space is connected, as we will prove below for the case n = 2. Using the curvature 
function for curves in R” that will be introduced in Sect. 4.2, it would not be difficult 
to modify the proof and show that any two curves y: [a,b] —> R” are regularly 
homotopic. 


Theorem 3.10 
Any two curves y, 7: [a, b] > R? are regularly homotopic. 


Proof. By transitivity, we can construct the desired regular homotopy in steps. As a 
first step we use a regular homotopy to achieve that y has length b — a: 


= (1-142) 
y= LY) y 


Therefore, without loss of generality we may assume that the original curve y 
already has length b-a. Then we can use a regular homotopy to achieve that y has 
unit speed: using the inverse of the arclength function s: [a,b] —> [0, L] of y 
(Definition 1.13), we define a regular homotopy 


w(x) = y = Ax +t s7"). 
So we can assume without loss of generality that y has unit speed. Now we use 
a regular homotopy in order to translate the starting point of y to the origin and 


achieve y (a) = 0: 


“%=U-hy@t+y. 
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Similarly, we can rotate y to achieve that the unit tangent 


nosa) 


of y at the starting point becomes the first standard basis vector e; of R? (Fig. 3.7): 


= ( cos((1 — 1)B) sin((1 — >) 
” = N- sim((1 — 1)B) cos((1 — 1)8) 


We apply the same normalizations to y. Now we consider the linear interpolation 
Ke = (1 — t)k + tk 


between the curvature functions «x of y and « of y and define the desired regular 
homotopy from y to 7 by 


Fig. 3.7 The initial regular homotopy that brings the curve into a standard position and size 
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3.6 Whitney-Graustein Theorem 


Definition 3.11 


A regular homotopy through closed curves between two closed curves 
y, 7: [a,b] > R? is a regular homotopy t +> y; between y and 7 such that for 
all t € [0, 1] the curve y; is closed. If there exists such a regular homotopy, y 
and y are called regularly homotopic through closed curves. 


Let us start with an example that will be needed below. Recall that y : [a, b] > 
R” is called closed if y = ?|{a,o) for some periodic smooth map y: R —> R” with 
period b — a. A simple way to make a new closed curve y: [a,b] :—> R” out of 
such a curve y is by a so-called parameter shift, which depends on a number t € R: 


7X) = 9-7). 


This closed curve y is regularly homotopic through closed curves to y, a suitable 
regular homotopy being t > y, with 


V(x) = P(x — tr). 


Like regular homotopy in Sect.3.5, regular homotopy as closed curves is an 
equivalence relation on the set of closed curves in R? and the equivalence classes 
can be thought of as the connected components of this space. This time however, 
the whole space is not connected: 


Theorem 3.12 ([45]) 
Two closed curves y, Ÿ : [a, b] > RR? are regularly homotopic through closed 
curves if and only if they have the same tangent winding number. 


In Fig. 3.8 we see an example of a regular homotopy through closed curves. 


Proof. Suppose there is a regular homotopy as closed curves t +> y; between y and 
y. Denote by ds; = |y;| and x; the speed and the curvature of y;. Then the tangent 
winding number 


1 b 


= 2 Ky ds; 
2m Ja 


nt 


is an integer for all £ € [0,1] and it depends continuously on t. Therefore it is 
constant and no = nı means that y and y they have the same tangent winding 
number. 
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Fig. 3.8 A sequence of curves from a regular homotopy between the elastic figure eight curve 
traversed twice (top left) and the elastic figure eight curve traversed only once (bottom right) 


Conversely, suppose that y and y they have the same tangent winding number. 
As in the proof of Theorem 3.10 we can assume without loss of generality that y and 
y both have unit speed. By Lemma 3.13 below and the fact that parameter shifts can 
be accomplished by regular homotopy as closed curves, we may also assume that 
the curvature functions « and K are either constant or linearly independent. As in 
the proof of Theorem 3.10, we can apply another regular homotopy through closed 
curves to achieve y (a) = 0 and y’(a) = e; (see Fig. 3.7). The same can be assumed 
for y. 

Let then t +» y be the regular homotopy between y and y constructed at the end 
of the proof of Theorem 3.10. The only problem is that for the intermediate curves 
yı might not be closed. We are going to repair this by modifying y; to a closed curve 
Yı as follows: 


= =a b 
Fats) = mtx) — Ef n 


The only fact that needs to be checked is that y/ (x) 4 0 for all x € [a, b]. 
Suppose we would have 


_ 1 
0= HG) =Te)- f T 
b—a Ja 
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and therefore 


1 
1 = |T (x)| = bza 


b 1 b 
T, S) IT] =1. 
Í t b—a t 


The inequality sign in the above formula must be an equality, and this implies that 
T; is constant, i.e. 


O=K = (l-t) +tk. 


This would imply that « and « are linearly dependent as functions, which by our 
assumptions means that « and x are constant. Since both coefficients in the previous 
equation are positive, this would imply x = K = 0, which is impossible for closed 
curves. Oo 


As a consequence of Theorem 3.12, every closed curve in IR? is regularly homotopic 
through closed curves to one of the curves in the following list in Fig. 3.6. 


We conclude this chapter with the Lemma that was needed in the proof of the 
Whitney-Graustein theorem: 


Lemma 3.13 
Let k: R — R be a non-zero periodic function such that for allt € R the 
functions k and x +> k(x — T) are linearly dependent. Then k is constant. 


Proof. Given our assumptions, there is a smooth function à: R —> R such that for 
all x € R we have 


K(x —T) =A(T)K(X). 
Differentiation with respect to t at t = 0 yields 
K' (x) = dV’ (O)K (x) 


The only non-zero periodic functions that satisfy such a differential equations are 
the constant functions. o 


Suppose we have a diffeomorphism g: M —> R? where 
M := {x € R? | |x| < 1} 


is the unit disk in R?. Then we can define a closed curve y: [a,b] > R? in sucha 
way that the Fig. 3.9 becomes a commutative diagram. 
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Fig. 3.9 A diffeomorphism 
g from the unit disk M into 
IR? and the corresponding 
boundary loop y 


If a closed curve bounds a region in IR? that can be mapped onto the unit disk by 
a diffeomorphism, its tangent winding number is one or minus one: 


Theorem 3.14 

In the setup of Fig. 3.9, the tangent winding number of y is £1, where the plus 
sign applies if and only if g preserves orientation, i.e. if det g'(x) > 0 for all 
xeM. 


Proof. Already in the proof of Theorem 3.12 we saw that applying a scale or a 
rotation to y does not change the regular homotopy class of y. Therefore, without 
loss of generality we may assume 


g’(O)ei = e. 
For t € [0, 1] let A; be the 2 x 2-matrix such that 


A;e, = ej 
Ar g'(e2) = (1 — t)g'(e2) + te 
where the plus sign is chosen if and only if det g’(0) > 0. Then the matrix A; is 
invertible for all ż and the one-parameter family t +> yp = Ary of closed curves is 


a regular homotopy, so after replacing g with A; o g we can assume without loss of 
generality that 


g'(0) =I. 
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Now define for r € (0, 1] closed curves y, : [0, 277] —> R? by 


_ cos(rx) 
w= 8 i) i 


For small € > 0 the curve ye is close the parametrization 


= cos(x) ) 
*F \ + sin(x) 


of to the unit circle, so the tangent winding number of ye is +1. On the other hand, 
y = yı is regularly homotopic to ye, and therefore also the tangent winding number 
of y is +1 (Fig. 3.8). o 
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Parallel Normal Fields 


For curves y: [a,b] — R” there is an analog «: [a,b] > IR”~'of the curvature 
function of a plane curve. In the context of unit speed curves, this function « 
determines y up to an orientation-preserving rigid motion of IR”. Before we can 
define «x, we have to study parallel normal vector fields along a curve in R”. 


4.1 Parallel Transport 


Definition 4.1 


Let y: [a,b] — R” be an immersion with unit tangent field T: [a,b] —> R”. 
Then a smooth map Z: [a, b] > R” is called a normal field for y if 


(Z(x), Tœ) =0 


for all x € [a, b]. The (n — 1)-dimensional linear subspace T (x)*+ is called the 
normal space of y at x. 


If Z: [a,b] > R” is a normal field for y, then we can split its derivative Z’ into 
its tangential part and its normal part: 


Z’=,~T+W 
where à: [a, b] —> R is a smooth function and W is another normal field. It turns 
out that A(x) can be computed from Z(x) alone, without taking the derivative of Z: 


differentiating the expression (Z, T} = 0 we obtain 


à = (Z', T) =—-(Z,T’). 
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The scalar product (Z’, Z) = 3(Z , Z)’ measures how the length of Z changes along 
y. A component of Z’ orthogonal to Z and T indicates a rotation of Z around the 
tangent T. If Z has constant length and there is no such twisting, Z is called parallel: 


Definition 4.2 


A normal field Z: [a, b] — R” along a curve y : [a, b] > R” with unit tangent 
field T : [a, b] —> R” is called parallel if there is a function 4: [a, b] —> R such 
that 


Z' =kT. 


There is a parallel normal field Z for every immersion y and all such fields come 
in an (n — 1)-parameter family: 


Theorem 4.3 
Given a vector W € T(a)+ in the normal space of a curve y: [a,b] > R” 
at a, there is a unique parallel normal field Z: [a, b] > R” of y such that 


Z(a) = W. 


If Z, Y are two parallel normal fields along y, their scalar product (Z, Y) is 
constant. 


Proof. If Z is a parallel normal vector field along y with Z(a) = W, then 
differentiating the equation (Z, T) = 0 yields (Z’,T) + (Z, T’) = 0 and, using 
Z' = —(Z,T’)T, we see that Z solves the linear initial value problem 
Z(a)=W 
Z' =(Z,T)T' —(Z,T')T. 

By the Picard-Lindelof theorem, such a solution is unique, which proves the 
uniqueness part of our claim. For the existence part, let Z be the solution of the 
above initial value problem. For any further solution Y of the above differential 
equation we have 

(Z, Yy =(Z',Y)+(Z,Y') 

= ((Z,T)T' —(Z,T'\T, Y) + (Z, (Y, T)T' — (Y, T’)T) 

0. 
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and therefore the scalar product (Z, Y) is constant. In particular, Y = T is such a 
solution, so (Z(a), T(a)) = O implies (Z, T} = 0. Therefore Z is a normal field, in 
fact a parallel one. o 


If y: [a,b] > R” is a curve and W is a vector in T(a)+, for every x € [a, b] we 
can use the parallel normal field Z with Z(a) = W to “transport” W to a normal 
vector Z(x) € T(x)+. This parallel transport map 

P(x): T(a)t > T(x)t+ 
is obviously linear, and by Theorem 4.3 it is in fact orthogonal, i.e. it preserves 
scalar products. Moreover, each normal space T(x)+ carries an orientation in the 
sense that a basis Wi, ..., W„—1 of T (x)+ is called positively oriented if 


det(W,,..., Wn—1, T(x)) > 0. 


If Wi,..., Wn—-1 is a positively oriented basis of T(a)t and Z),..., Zn,—1 are 
parallel normal fields with Z;(a) = Y; then the map 


x > det(Z1(x),..., Zn-1(x), T(x)) 
is continuous and never zero. Therefore, for all x € [a, b] we have 
det(Z, (x), ..., Zp_1(x), T(x) > 0 
and the map P(x) is orientation-preserving. We summarize this as follows: 


Definition 4.4 


Given a curve y: [a,b] — R” and x € [a,b], the orientation-preserving 
orthogonal map P(x): T(a)+ — T(x)+ defined above is called the parallel 
transport from the normal space T (a)+ to the normal space T(x)+. 


By Theorem 4.3, each vector Z(x) of a parallel normal field has the same length. 


Therefore, we can use parallel normal fields Z in order to displace a curve y by a 
fixed distance € = |Z|, without introducing unnecessary twisting: 


Definition 4.5 


If Z is a parallel normal field along a curve y : [a, b] —> R” and the derivative of 
y=yv+Z 


vanishes nowhere, then the y is called a parallel curve of y. 
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For a curve y : [a, b] — R” the continuous (but not necessarily smooth) function 
dT 
— | : [a,b] > R 
| T; [a, b] 

is called the absolute curvature of y. If € > 0 is such that 


1 dT 
-> max | Zo 
€ ds 


elabi} 


and Z is a parallel normal field with |Z| = e then by the Cauchy-Schwarz inequality 


we have 

Z T 

i 1+ z T))r#o 
ds ds 


Therefore, if we pick a vector W € T (a) with sufficiently small norm and define 
Z as the parallel normal field Z with Z (a) = W, then y + Z will be a parallel curve 
for y. 

As an application, we always visualize a curve in R? by thickening it, which 
means that we chose a suitable collection of W € T(a)+ with small length and 
draw the union of the corresponding parallel normal fields. Most of the time we use 
a small circle centered at the origin in T (a)*, but different choices (as in Fig. 4.1) 
are also possible. 


Fig. 4.1 A closed curve in the normal space T (a)+ is used to build a thickened version of y by 
parallel transport 


4.2 Curvature Function of a Curve in IR” 51 


4.2 Curvature Function of a Curve in R” 


We saw in Sect.3.1 that, up to rigid motions of RÊ, the geometry of a unit 
speed curve y: [a,b] — R? is completely determined by its curvature function 
k: [a, b] > R. Here we will define a similar curvature function « : [a, b] > R”! 
for any unit speed curve y: [a, b] > R”. To define «(x), we use parallel transport 
to transfer the normal vector T’ (x) € T(x)+ to the normal space T (a)*. Afterwards 
we use an orthonormal basis of T (a)+ in order to identify T (a)* with R”7!. 


Theorem 4.6 

Let y: [a,b] — R” be a curve with unit tangent T and parallel trans- 
port maps P(x): T(a)+ —> T(x). Then there is a unique smooth map 
Y: [a,b] > T(a)+ such that for all x € [a, b] we have 


dT 
Pœ) = ———(). 
ds 
W is called the Hasimoto curvature of y. 


See Sect. 5.3 for the details on Hasimoto’s contribution. The Hasimoto curvature 
determines y uniquely: 


Theorem 4.7 

Given a point p € R”, a unit vector S € R” and a smooth map Y: [a,b] > 
T(a)+, there is a unique unit speed curve y: [a,b] > R” such that y (a) = 
p, y' (a) = S and W is the Hasimoto curvature of y (see Fig. 4.2). 


Fig. 4.2 The Hasimoto —= 
curvature W of a curve y N 
indicated as a blue curve in Ta)" | y 


T(a)t 
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Proof. First we prove uniquess of y. Let y : [a, b] > R” be a curve with the desired 
properties. Choose an orthonormal basis Wj, ..., W,-1 of T (a)* such that 


det(W,..., Wn_-1, T(a)) = 1 
and define «K1, .. . Kn—1 by 
Y =k Wi +... + Kn-1Wn-1. 


Let Z1, ..., Zn—1 be the parallel normal fields along y such that Z; (a) = Wj for 
all j e {1,...,n — 1}. Then 


(Zi, tiwy Zn-1; T): la, b] => R"*” 
solves the initial value problem 


(Zi, ae iy Zn-1; T)(a) = (Wi, Sass Wr-1, S) 


n—-1 


(Z1,..., Zn—1, TY = | «T,..., KnT, -9 «Zj 
j=l 


and is therefore uniquely determined by p, S and W. In particular, T is uniquely 
determined and so is 


re yas f r. 


For existence, we can use the above initial value problem to define the vector fields 
(Z1,..., Zn-1, T). Atx = a these vectors are orthonormal and their pairwise scalar 
products solve the system of linear differential equations 


(Zi, Ziy = «i(T, Zj) + «j(Zi, T). 


We can interpret this as an initial value problem for the functions (T, Z;), (T, T) and 
(Zi, Zj). The functions (T, Z;) = 0, (T, T) = 1, (Zi, Zj} = ôij solve this initial 
value problem, and by Picard and Lindelöf such a solution is unique. Therefore, 
(Z1,..., Zn—1, T) stay orthonormal. So by integration of T we obtain a unit speed 
curve y: [a,b] > R” with y (a) = p and y'(a) = S. Zi, ...Zn—1 are parallel 
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normal fields along y with Z;(a) = Wj. Because we already know that T” = 
— Xi KjZj, this implies that Y is indeed the Hasimoto curvature of y. o 


In the above proof we used a basis of T (a)* in order to turn Y into an R”~!-valued 
function x. This function is the promised analog of the curvature function of a plane 
curve: 


Definition 4.8 


Let y: [a,b] — R” be a unit speed curve with unit tangent T and Hasimoto 
curvature W. Let W1, ..., W,-1 be a positively oriented orthonormal basis of 
T (a)*. Then the function 


K1 


Kn—1 
defined by 
Y = ki Wi +... + Kn—-1 Wn—1 
is called a curvature function of y. 


In the case n = 2 the positively oriented orthonormal basis of T(a)+ mentioned 
in the above definition is unique, and therefore each plane curve has a unique 
curvature function x: [a, b] > R! = R, which is the one we already encountered 
in Sect. 3.1. It is clear from its definition that for any n the function « is at least 
unique up to a rotation of R”~!: 


Theorem 4.9 

If, R: [a,b] > R”! are curvature functions of the same curve y : [a, b] > 
R”, then there is an orthogonal ((n — 1) x (n — 1))-matrix A with det A = 1 
such that 


k = Ak. 


On the other hand, as in the case of curves in R?2, for every curvature function 
k: [a, b] > R”! there is a corresponding curve y: [a,b] > R” and y is unique 
up to post-composition with an orientation preserving rigid motion of IR”. Also the 
following theorem is a direct consequence of Theorem 4.7: 
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Theorem 4.10 

Given a smooth function x: [a,b] — RR""!, there is a unit speed curve 
y: [a,b] = R” for which k is a curvature function. The curve y is unique 
up to an orientation preserving rigid motion of R”, which means that if y is 
another curve having k as a curvature function, then there is an orthogonal 
(n x n)-matrix A with det A = 1 and a vector b € R” such that 


ý = Ay +b. 


4.3 Geometry in Terms of the Curvature Function 


Let y: [a,b]—> R” be a unit speed curve with unit tangent field T and 
W,..., Wn—1 a positively oriented orthonormal basis of T(a)t. Let Z1,..., Zn—1 
be the corresponding parallel normal fields along y with Z;(a) = Wj. Then we can 
describe every normal field Y along y in terms of a function y: [a, b] > R” as 


Y = y1 Zi +... Yn-1Zn-1 


| | yı 
= | Zi = Zy-1 : 
| | Yn-1 
=: Ny 
where for x € [a, b] the matrix N(x) has the vectors Z1(x),..., Zn—1(x) € R” as 


its column vectors. In terms of the curvature function « introduced in Definition 4.8 
the derivative of Y can be expressed as 


Y’=(k,y)T+Ny’. 
In particular, for Y = T’ we obtain 


T’ = -Nk 
T” = —(k,«)T — Nk’ 
T” = —-3(k,«’)T +N ((k, K)K — K") 


Now we are able to generalize the results we obtained in Sect. 3.3 for plane curves: 
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Theorem 4.11 
A unit speed curve y : [a,b] + R” is torsion-free elastic if and only if there 
is a constant à € R such that its curvature function k satisfies 


2 
te tae =0. 


Proof. By Theorem 2.23, y is torsion-free elastic if and only if there is a constant 
à € R such that 


3 
0=T" +3(T',TT + zT T’)T' —AT' 
Ik]? 
=-N (e + Ea +1) ; 


oO 


Here are further examples of how the geometry of y is reflected in the properties of 
K: 


Theorem 4.12 


Let k: [a,b] — R”! be a curvature function of a unit speed curve 
y: [a,b] > R”. Then: 


(i) k = 0 if and only if the image of y lies on a straight line. 
(ii) k is a non-zero constant if and only if the image of y lies on a circle. 
(iii) The image of k lies in a hyperplane through the origin of R"! if and 
only if the image of y lies in a hyperplane of R”. 
(iv) The image of k lies in a hyperplane of R”™! that does not pass through 
the origin if and only if the image of y lies in a hypersphere of R” (see 
Fig. 4.3). 
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Fig. 4.3 The curvature 
function « of a curve on S? 
lies on a straight line which 
does not pass through the 
origin 


R2 


Proof. Claim (i) is obvious, since the image of a curve lies on a straight line if and 
only if its unit tangent T is constant. If the image of « lies in a hyperplane through 
the origin of R”~!, there is a unit vector a € R”~! such that (a, «) = 0. Then 
(Nay = —(k,a) =0 
so there is a fixed vector n € R” such that Na = n. We have 
(n, y) = (Na, T) =0 
and therefore the image of y is contained in a hyperplane with normal vector n. The 


proof of the converse is left to the reader. This establishes (iii). For (iv), suppose that 
there is a unit vector a € R”! and a number r > 0 such that 


Then 
(y —rNa) =T —r(k,a)T =0 
so there is a fixed point m € R” such that 
y—-rNa=m 
and we have 
ly —m| =r. 
Therefore, the image of y lies on the hypersphere with center m and radius r. Again, 


the proof of the converse is left to the reader and we have established (iv). For (ii) 
we use induction on n based on (iii), starting at n = 2 where we use (iv). oO 
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Curves in R? 


For a closed curve y : [a,b] + R? with unit tangent field T, we can use a parallel 
normal vector field to transport a normal vector W; € T(a)* at the starting point of 
y to anormal vector Wp € T(b)+ at the end point. If y is closed, the angle T (y) 
between W, and W; is called the total torsion of y. A notion of total torsion can also 
be defined for curves in IR? that are not necessarily closed. Therefore, for curves in 
R3, total torsion provides another geometric functional besides length or bending 
energy. Critical points of a linear combination of length, total torsion and bending 
energy are needed for modelling the physical equilibrium shapes of elastic wires in 
R3. 


5,1 Total Torsion of Curves in R3 


Let us focus now on curves y: [a,b] —> R?. In this case we can visualize the 
parallel transport of normal directions introduced in Sect.4.1 as approximately 
implemented by a chain of so-called “constant velocity joints” (cf. [35]). Such joints 
are able to transport normal directions in an angle-preserving manner. Rotating the 
initial vector Z (a) of a parallel normal field by an angle œ will make the final vector 
Z(b) rotate by the same angle a (see Figs. 5.1 and 5.2). 


Definition 5.1 


For a curve y : [a, b] + R? the orthogonal linear map 
P: T(a) > T(b)*, XB Z(b) 


where Z: [a,b] —> R? is the parallel normal field along y with Z(a) = X is 
called the normal transport of y. 
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Fig. 5.1 Two curves built from constant-velocity joints. The total torsion is zero for the curve on 
the left, but not for the one on the right. The red line indicates a parallel normal field 


Fig. 5.2 A trefoil knot cannot be built from constant-velocity joints. Because of the angle-defect 
due to the total torsion, the joints would not close up 
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After having chosen a pair W = (Wa, Wp) of unit vectors Wy € T(a)* and 
Wp € T(b)+ we can describe the normal transport P by an angle: 


Definition 5.2 


Let y : [a, b] > R? be curve with unit tangent T and normal transport P. Then, 
given a pair W = (Wg, Wp) of unit vectors Wa € T(a)+ and Wp in T(b)+, the 
unique angle 


Tw € R/anz 
with 
P(Wa) = (cos Tw) Wp + (sin Tw) T(b) x Wy. 


is called the total torsion (cf. Sect. 5.4) of the curve y with respect to W, and 
Wp. 


For a closed curve we have T(a) = T(b) and we can always choose W(b) = 
W (a). The total torsion 7w then becomes independent of the choice of Wz, so in 
this case we can drop the subscript W and denote the total torsion of a closed curve 
y by T (y). Let us determine the infinitesimal variation of the total torsion Tw if we 
vary the curve y as well as the unit vectors Wa and Wp: 


Theorem 5.3 

Lett +» y be a variation with variational vector field y = Y of a curve 
y: [a,b] > R?. Lett > W,(t) € T(a)* and t œ> Wp(t) € T, (b)+ be two 
smooth families of unit vectors. Then the total torsion Tw (t) of y; with respect 
to 


W(t) = (Wa (t), W(t) 


satisfies 


b 
T Tw (t) = (Wa, T (a) x Wa) — (Wp, T (b) x Wp) -f det(T, T’, T). 
t=0 a 


In terms of Y, the above integral can be expressed as 


b ae dT b b dT 1 
det(T, T, T) = (Y, T x — } — yI x | = ; 
HS Ihe dhe ds 


a 
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Proof. Let Z, be the parallel normal field along y; with Z; (a) = Wa (t). In particular, 


Zo =: Z is a parallel normal field and we have 
Z' = —(Z,T')T. 
Taking the time derivative of the equation 
Z, (b) = cos(Tw (t)) Wb (t) + sin(Tw(t)) T (b) x Wp (t) 
at t = 0 yields 
(Z(b), T(b) x Z(b)) = Tw + (Wy, T(b) x Z(d)). 
From Theorem 5.4 
b b. 
f det(T, T’, T) = -f (Z,T x ZY 
a a 


= —(Ż(b), T (b) x Z(b)) + (Ż (a), T (a) x Z(a)) 


= —İw — (Wp, T (b) x Z(b)) + (Wa, T (a) x Z(a)). 


The second claim is a consequence of 


s dY 
det(T, T’, T) = det (7. TE =) 
ds 


dT .\' dT\N' 
= det | T, —, Y | — det | T, | —),Y], 
ds ds 


where we used that T = gY — (Z, T)T (cf. proof of Theorem 2.8). 


oO 


The following Theorem was used in the above proof and will be needed also in 


Sect. 5.3: 


Theorem 5.4 


Let t œ> y be a variation of a curve y : [a,b] > R? with unit tangent field 
T. Lett +» Z, be a smooth one-parameter family of maps such that Z, is a 


parallel unit normal field along y;. Then 


(Zeb x Z} = (Pil x T^. 
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Proof. The prime derivative commutes with the dot derivative, hence 
(Z,T x Zy =((Z')',T x Z)+(Z,T' x Z) 
= (-(z, TT, T x z)+ (Z,T' x Z). 
Because (Z;, Z;) = 1 and (Z;, T;) = 0, for all t, we have 
Ż=(Ż,T)T+(Ż,Tx DT xZ 
=-—(Z,Î)T+(Ż,T x Z)T xZ 
and therefore we can continue the previous calculation of (Z, T x Z)' as follows: 
(Z,T x ZY =(Z,T)T x Z,T')— (Z, TT x Z,T) 
=(Zx(T x Z),T xT’) 
=(T,T xT’). 
| 


In Sect.4.3 we described normal vector fields Y along a curve 
y: [a,b] — R” in terms of functions y: [a,b] — IR”~!. Given a parallel unit 
normal field Z along a curve y : [a, b] —> R? with unit tangent 7, in terms of this 
correspondence, any normal field Y can be written as 


Y=Ny=y1Z+yoT xZ. 


Then the function (Z, T x Z) in Theorem 5.3 also appears if, given a variation of y 
and Z, we want to know the time derivative of the above formula: 


Theorem 5.5 

Lett œ> y be a variation of a curve y : [a,b] > R? with unit tangent field 
T. Lett œ> Z; be a smooth one-parameter family of maps such that Z, is 
a parallel unit normal field along y, and t +> y; a smooth family of maps 
yz: [a, b] > R2. Then 


(Noy = (Ny, TT (54+ (Z,T x Z)Jy) 
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Proof. Taking into account the time derivatives of the equations that tell us 
T,Z,T x Z are orthonormal, we obtain 


(Ny) = N+ yı UŽ, T)T + (Z,T x Z)T x Z) 
+ yo((T x Z), T)T + ((T x Z), Z)Z) 


= —(Ny, ÈT +N (5+ (ZF x Z)Jy) , 


5,2 Elastic Curves in R? 


The torsion-free elastic curves studied in the Sects. 2.4 and 2.5 were critical points 
of bending energy under the constraint of fixed length. For general elastic curves 
in R? also the total torsion is constrained (see Fig. 5.3). Note that for a variation 
t +» yı with support in the interior of [a, b] of a curve y: [a, b] > R? with unit 
tangent T all parallel transport maps P, are defined on the same vector space T(a)+. 
Therefore it makes sense to consider the derivative a | 0 P;. 


Fig. 5.3 Elastic curves obtained by minimizing bending energy under the constraint of fixed 


length and fixed total torsion, for various values of the total torsion constraint: 0, 2x, Sx, gt 


3, 20 
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Definition 5.6 


A curve y: [a,b] — R? is called an elastic curve if it is a critical point of 
the bending energy 6 under the constraint of fixed length £ and fixed normal 
transport P. 


During a variation t œ> y, of a curve y: [a,b] — R? with constant support in 
the interior of [a, b], after choosing unit vectors Wa € T(a)+ and Wp € T(b)*+ 
(independent of t), we can measure the normal transport along y; as the total torsion 


angle Tw (y;). Theorem 5.3 then will tell us the infinitesimal variation of the normal 
transport, in a way that does not depend on the choice of Wa and Wp. 


Theorem 5.7 
Let y: [a,b] > R? be a unit speed curve with unit tangent T. Then the 
following are equivalent: 


(i) y is an elastic curve. 
(ii) There are constants à, u € R such that 


T” — (T”,T)}T + a T^T'— uT x T" —aT’ =0. 
(iii) There are constants à, u € R and a constant vector a € R? such that 
T” + a T'\T —pT x T’-AT +a=0. 
(iv) There is a constant u € R and a constant vector a € R? such that 
T” —(T" T)T +a—(a,T)T -uT xT =0. 
(v) There is a constant u € R and constant vectors a, b € R? such that 
y xy" = uy" +a x y +b. 
(vi) There are constant vectors a, b € R? such that 


y"=(axy+b)xy’. 


Proof. Theorem 2.21 was formulated in such a way that it is capable of dealing 
with several constraints, so that it is possible to prove Theorem 2.20 with more 
than just a single constraint. So more constraints than just the length are possible in 
Theorem 2.23. Therefore, the equivalence of (i) and (ii) can be shown following the 
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same arguments that lead to Theorem 2.23. Here, taking derivatives of the equation 
(T, T) = 1 gives 


iT" T} = —3(T, 7"), 


The equivalence of (ii) and (iii) follows from the equality 
3 + 
(r+ zT TT —pT xT -ar) 

A A 3 / / / MN / 

=T" — (T A ie 


which can again be verified by taking derivatives of the equation (7,7) = 1. (iv) 
is just the component of (iii) orthogonal to T, so it follows from (iii). In order to 
show that (iv) implies (iii), we have to show that (iv) implies that there is a constant 
à such that also the equation obtained by taking the component of (iii) parallel to T 
is satisfied if (iv) holds, which is indeed the case: 


3 / 1 1 

(e. T) T zT T’) T (a, r)) = Gr. T’) F (a, r)) 

= (T" +a,T’) 

= (uT x T',T') 

=0. 

To prove that (iv) is equivalent to (v), note first that, the left-hand side of (iv) being 
orthogonal to T, (v) is equivalent to the equation obtained from (iv) by taking the 
cross product with T: 


TxT"”"+Txa+uT'=0 
or 
0=(TxT'+yxa+uTY. 


Therefore, (iv) is equivalent to (v). Taking the cross product of the equation in (v) 
with T yields 


—y"=-(axy)xy’=bxy’ 


which is equivalent to (vi). To show that (vi) implies (v), note that the component 
orthogonal to T of the equation in (v) is equivalent to (vi). This means that we have 
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Fig. 5.4 Kirchhoff showed that, as an elastic curve is traversed with unit speed, its tangent vector 
T follows the motion of the axis of a gyroscope. The photograph with long-time exposure was 
reproduced from [41] with permission from Javier Villegas 


only to show based on (vi) that the scalar product with T of the sum of the terms 
without the uT term is constant. This is indeed the case: 


—(T,y x a+b) = (T, T’) =0. 
oO 
In 1858 Gustav Kirchhoff realized (cf. [19]) that in the form (ii) or (iii) the equations 
for T describe the motion of the axis of a heavy symmetric top (or gyroscope). The 
vector a describes the direction of gravity and u is related to the spinning speed of 


the gyroscope (see Fig. 5.4). Figure 5.5 illustrates an interesting special case of the 
characterization given in part (vi). of Theorem 5.7: Assume b = 0, a 4 0 and 


(e3, a) = (e3, y (a)) = (e3, y'(@)) = 0. 
Then equation in (vi) can be written as 


y” =(y',a)y —(y',y)a 
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Fig. 5.5 A unit speed elastic curve y can be described as the orbit of a charged particle moving 
in a linear magnetic field p + B(p) = a x p+ b. If the initial velocity T (a) is orthogonal to 
B(y(a)), the elastic curve lies in a plane 


and therefore the function 
g: [a,b] > R, g = (e3, y) 


satisfies the linear second order equation 


n 


g" = (y',a)g 


with the initial condition g(a) = g'(a) = 0. It follows that g vanishes identically 
and the image of y is contained in the plane R? cC R? given by E = {(x, y, z) € 
R? |z = 0}. Assuming that y has unit speed, T is the unit tangent of y and « its 
curvature, we can rewrite the equation in (vi) further as 


kJT =(T,a)y —(T,y)a 


_ (y, Ja) 
~ (a,a) 


(y, Ja) JT. 


({(T, a) Ja — (T, Jaja) 


The second of the above equalities can be verified by expanding y as 


y=ly sha ty Jalsa 
’ Jal / lal * fal/ > fal 
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while the third equality follows by a simillar expansion of T. This means that «x (x) 
is proportional to the distance of y(x) to the line in RŽ through the origin with 
direction fa: 


k = (y, Ja). 


Amazingly, this is exactly the description of planar elastic curves that had been 
given by Jakob Bernoulli in 1691 (see [25] for a historical survey, or [14]). 


5.3 Vortex Filament Flow 


Vortex filaments are curves of singularly concentrated vorticity in a moving fluid. 
Familiar examples are tornados and smoke rings. The mathematical theory of vortex 
filament motion started with Lord Kelvin, who in 1880 investigated the evolution of 
small pertubations of a straight vortex filament (cf. [39]). Later, these perturbations 
were called Kelvin waves. The full evolution equation for thin vortex filaments was 
found in 1906 by Tullio Levi-Civita and his student Luigi Sante da Rios. For a 
detailed history see [34]. Mathematically, the motion of a vortex filament can be 
described by a one-parameter family y; of curves and the da Rios equation says that 
this one-parameter family satisfies (Fig. 5.6) 


a ih EN 
= ds? ` 


Another breakthrough occurred in 1972 (cf. [15]) when Hidenori Hasimoto showed 
that the da Rios equation is equivalent to the non-linear Schrödinger equation: 


Fig. 5.6 A curve evolving according to the da Rios equation adapted from [20] with permission 
from William Irvine 
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Theorem 5.8 

Lett +> y, witht € [t), t2] be a smooth one-parameter family of unit speed 
curves yr: [a,b] > R? which solves the da Rios equation. Let T; be the unit 
tangent field of y+. Then there is a smooth family t +> W; of unit vectors W, € 
T;(a)+ such that the corresponding family t œ> «; of curvature functions for 
the curves y; satisfy 


é |r|? 
JKt + Kf SF 


Proof. We choose an arbitrary unit vector We Ti, (a)+ and define for t € [ż, t2] 
a family of unit vectors W, € T,(a)+ as the solution of the linear initial value 
problems 
Wi = WwW 
š . 1 
W; = —(W:, T:(a))Ti(a) — 5 (1), t, (a)) T; (a) x Wi. 


Let Z, be the parallel normal field along y; with Z; (a) = W;. By Theorem 5.4 


: 1 j ` 
(Ca T; x Zi) + zT 1) = (T, Ti x TH) + (T, T) 


= (T,, (T, x i) sy Ty + (T), eg 
= 0, 


where we used that the assumption that y; solves the da Rios equation for all t 
implies 


t=T xT". 


By construction we have W, = Z; (a), hence 


‘ 1 
(cz, T; x Zi) + rie m) (a)=0 


and therefore 


: 1 
(Zi, Ti xX Z) = a T} 
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Using the formulas for T” and T” from Sect. 4.3 and Theorem 5.5, it then follows 


that 
: 2 í 


(Ni) 


=-T'x TT! = i x pit 


= k| T! x THT x NCP — Kl) 
=T x (Nk) 
= N(Jx;') modT, 


where N is again the matrix of parallel normal fields which is used for the definition 
of k. This implies that t +> «; satisfies the nonlinear Schrödinger equation. Oo 


The nonlinear Schrödinger equation was known to be a so-called Soliton 
equation, and as a consequence also the da Rios equation admits infinitely many 
constants of the motion. Finally, in 1983 Marsden and Weinstein established (cf. 
[28]) vortex filament motion as a Hamiltonian mechanical system in its own right 
(see [10] for a survey article). 

The closed curve in Fig.3.4 is a critical point of bending energy under the 
constraint of fixed length and fixed enclosed area (cf. [2, Figure 8]) . It can be shown 
that this curve is the initial curve yo of a solution t +> yr of the da Rios equation 
that is defined for all times ¢ € R. In fact, it can also be proven that this solution 
is periodic in t. This solution (shown in Fig. 5.7) matches quite well the qualitative 
behavior of the vortex filament shown in Fig. 5.6. 


Definition 5.9 


A vector field X: R? — R? is called an infinitesimal rigid motion if there are 
vectors a, b € R3 such that for all pe R? we have 


X(p)=axp+b. 


The following is a reformulation of part (v) of Theorem 5.7: 
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Fig. 5.7 An initial curve which evolves according to the da Rios equation matches the qualitative 
behavior of a vortex filament (cf. Fig. 5.6). This is why the resulting flow is also referred to as 
“vortex filament flow”, or “smoke ring flow” 


Fig. 5.8 By Theorem 5.7, an elastic curve that evolves according to the da Rios equation will just 
undergo rigid motions. These rigid motions can be pure translations (/eft), pure rotations (right) or 
screw motions (middle) 


Theorem 5.10 

A curve y: [a,b] > R? is elastic if and only if there is an infinitesimal rigid 
motion X such that for every point of the curve the velocity y prescribed by 
the da Rios equation is given by evaluating X at that point: 


y = COW, 


Figure 5.8 shows closed elastic curves and their evolution under the da Rios 
equation. 
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5.4 Total Squared Torsion 


For pioneers of Differential Geometry like Jakob Bernoulli and Leonard Euler 
(cf. [25] for a historical survey), the motivation for studying elastic curves was to 
determine the shape y of a perfectly elastic thin wire (originally shaped as a straight 
line segment of fixed length when it came out of the factory). In a stable equilibrium 
position of such a wire the elastic energy stored in its deformation is minimized. 
Since Bernoulli (1691) and Euler (1744) focused on plane curves, bending energy 
as introduced in Sect. 1.3 was sufficient for modeling elastic energy. 

Later, Lagrange (1788, cf. [22]) and Binet (1844, cf. [5]) realized that for wires in 
space bending energy alone does not account for all relevant contributions to elastic 
energy. See also [40] and [4]. 

In R? one has to take account of the internal twisting of the wire: imagine a 
parallel normal field marked as a colored line on the surface on the wire in its 
original straight shape. If we bend the wire into space and want to know the elastic 
energy stored in the deformation, it is not enough to know the resulting shape y of 
the wire. We also have to know where the colored line goes on the deformed curve. 
The twisting of the wire made visible by the colored line contributes to the elastic 
energy, even if the curve is not changed at all (cf. Fig. 5.9). This means that elastic 
wires in IR? are more adequately modelled as a framed curve: 


Definition 5.11 


A framed curve in R? is a curve y: [a,b] > R? together with a unit normal 
field N along y. 


Fig. 5.9 Top: An elastic wire with no twist, as indicated by the red line. Bottom: The same wire 
in the same shape, only twisted (Reproduced from [14] with permission from Geoff Goss) 
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Instead of drawing many arrows, we will usually indicate the unit normal field N 
of a framed curve by marking a colored line on a slightly thickened version of the 
curve. 


Definition 5.12 


For a unit normal field N along a curve y: [a,b] —> IR, the function 
T: [a, b] > R given by 


is called the torsion of N. 
The torsion t measures the deviation of N from being a parallel normal field. 


After choosing unit vectors W, € T(a)* and Wp € T(b)+, we can assign a total 
torsion angle also to a framed curve: 


Twly, N) := p-a 
where the angles w, 8 € R/27z are defined by 


N (a) = cosa Wa + sing T (a) x Wa 
N (b) = cos 6 Wp + sin 8 T (b) x Wp. 


Tw (y, N) is related to the total torsion Tw (y) of the curve y itself as follows: 
Theorem 5.13 
Let N be a unit normal field with torsion t along a curve y: [a,b] > R? 


with unit tangent field T. Then, for any choice of unit vectors Wa € T(a)+ 
and W, € T (b)+ we have 


b 
Two N) =Twr)+ f E 


Proof. Let Z be the parallel normal field along y with Z (a) = Wa. Then there is a 
unique function 7: [a, b] > R with n(a) = œ such that 


N=cosn Z+sinn T x Z. 
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dN TSN 

T=(—,T x 
ds 
dn . dn 

= (—— sinn Z + —cosn T x Z,—sinn Z+cosyn T x Z 
ds ds 

a 

~ ds’ 

Furthermore, 


Z(b) = cos Tw(y) Wp + sinTw(y) T(b) x Wp 
and therefore 


cos Ê W,+sinB T(b) x Wp 
=N (b) 
=cos n(b) Z(b) + sinn(b) T (b) x Z(b) 
=cos (Tw (y) + n(b)) Wp + sin (Tw (y) + n(b)) T(b) x Wp. 


This means that 


Tw(y, N)+a =£ 
= Tw (y) + n(b) 


bq 
= T(y) + nla) + | = 
q as 


b 
=T) +a | T. 


Theorem 5.13 also explains why we use the terminology “total torsion”. 

Figure 5.9 (taken from [14]) shows on the bottom a configuration where the end 
points of the wire are still the same as in the relaxed configuration and therefore, 
in view of the fixed length, the curve y is still a straight line segment. However, 
additional energy has been stored in the twisting of the frame. In Fig. 5.10, moving 
the end points closer together has made it possible to for the wire to move away 
from the shape that would minimize bending energy in order to reduce its internal 
twisting. 
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Fig. 5.10 The same wire as in Fig. 5.9. The wire is still twisted by the same amount, but the 
endpoints are moved closer together (Reproduced from [14] with permission from Geoff Goss) 


One can show that in the limit of thin wires (where the thickness tends to zero) 
this additional energy is of the form 


cS(y, N), 
where c is a positive constant and S (y, N) is defined as follows: 


Definition 5.14 


Let N be a unit normal field with torsion t along a curve y: [a, b] > R?. Then 
the total squared torsion of the framed curve (y, N) is defined as 


1 b 
S(y,N) = J t? ds. 
, 2 


The constant c depends on material properties and on the thickness r of the wire. 
Following [37], we call c the twisting modulus. We work in units where the bending 
energy is given as in Definition 1.19. Starting from the formulas for the restoring 
torque and the bending stiffness, one finds that 


G 1 


“FE 2(1+v) 


where G is the shear modulus of the wire material, E is the Young modulus and 
v is the Poisson ratio. According to a table (cf. [9]) of Poisson ratios for common 
materials, the dimensionless constant c lies between 4 and 7: For example, copper 
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Fig. 5.11 Due to the different twisting moduli, different amounts of torsion are needed to form 
the same curve shape from a copper wire (top, c = 2), or a DNA-strand (bottom, c = 2) 


wires have c = Š. Also the twisting and bending of DNA strands (where there 
is no real “material”) can be modeled in the same way, see equation 4.1 of [37]. 
Depending on the ambient conditions, we have 4 < c < 2 (see Fig. 5.11). 
Fortunately, as we will see in Sect.5.5, the specific value of c is irrelevant for 
the possible shapes of elastic curves, i.e. of those curves y of a given length that 
are critical points of the total elastic energy B + S. The value of c only effects the 
normal field N that goes together with such a curve y, not the shape of y itself. 


5.5 Elastic Framed Curves 


In Sect. 5.4 we looked at the elastic energy (including the part that is due to internal 
twisting) stored in a perfectly elastic wire (modeled as a framed curve (y, N)) that 
came out of the factory as a straight line segment. Here we will show that the for 
an energetic equilibrium configuration of such a wire the curve y is an elastic curve 
(Definition 5.6) and the torsion t of the unit normal field N is constant. 


Definition 5.15 


Let y : [a, b] > R? be a curve and N a unit normal field along y. Then a smooth 
one-parameter family t > (yr, N;) of framed curves is called a variation with 
support in the interior of [a, b] if y(x) = y(x) and N; (x) = N(x) for all x 
near the end points of the interval [a, b]. 
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Definition 5.16 


A framed curve (y, N) is called an elastic framed curve with twisting modulus 
c > Oif 


dt 


(B) +eS(%, Ni)) =0 
=0 


t 


for all variations t > (yr, N;) of (y, N) with support in the interior of [a, b] 
which fix the length, i.e. for which 


dt 


L(y) = 0. 
=0 


t= 


Which curves y in R? can be supplemented by a unit normal field N in such a 
way that (y, N) is an elastic framed curve? It turns out that those curves are precisely 
the elastic curves: 


Theorem 5.17 
A framed curve (y, N) in R? is elastic with twisting modulus c if and only if 
its torsion t is constant and y is a critical point of 


BretT 


under the constraint of fixed length £. 


Proof. Let (y, N) be an elastic framed curve elastic with twisting modulus c. Let us 
first consider special variations t œ> (yr, N+) of (y, N) with support in the interior 
of [a, b] for which the curve itself does not move at all, i.e. for all t we have y; = y, 
so that for those variations we have 


d 
ae B = 0. 
dt = (vı) 


The normals that we consider are of the form 
N; = cos(ta)N + sin(ta)T x N. 


where a: [a,b] > R is a function with support in the interior of [a, b]. Then 
N=aT x N,T =Q, ds = 0 and 


A am T nv) 
t=(—,Tx 
ds 
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da dN 

= —(T x N,T x N)+a(—,T x (T x N) 
ds ds 

_ da 

~ ds” 


Therefore, for all such functions œ we have 


So 
I 
| 


(B) +¢S(%, Nr) 
t=0 
b _2 

Tis 

2 


d 
ee 
dt |o Ja 


b 
=c f tt ds 
a 


and therefore we must have t’ = 0. Let now t }> y; be an arbitrary variation of y 
with support in the interior of [a, b] for which 


ai L(y) = 0. 


t=0 


Then, for small t, we can define unit normal fields N, along y; (equal to N near the 
end points of the interval [a, b]) by projecting N(x) to T; (x)} where T; is the unit 
tangent field of y;: 


N —(N,T,)T, 


N; => ———— . 
IN — (N, T;)Tr| 
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Then, by Theorem 5.13 and with 4 


dt Pars) = 0, we have 


d d b 
2 Bran Bort f ids 
=0 a 


dt |, dt | 9 
d b T? 
=F _ Bante | ttds tez £ 


d 
F (Biv) + ¢S(%, N)) 
t t=0 


This proves the “only if” direction of our claim. We leave the “if” direction to the 
reader. H 


5.6 Frenet Normals 


Definition 5.18 


A unit normal field N : [a,b] —> R” along a curve y: [a,b] —> R” with unit 
tangent T is called a Frenet normal field if there is a function x: [a,b] > R 
such that 


If we ignore the effects of gravity, the unit vector pointing upward in the reference 
frame of an airplane like the one in Fig. 5.12 (which is lacking a rudder) will be a 
Frenet normal for its flight path (see Fig. 5.13). A curve y : [a, b] > IR? has exactly 
two unit normal fields, and both of them are Frenet. The one with N = —JT has 
kf = k Where « is the curvature function of y. 


Fig. 5.12 The vertical vector in the reference frame of an airplane with no rudder is a Frenet 
normal along its flight path 
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\ 


Fig. 5.13 A Frenet normal along a space curve 


Not every curve in IR? has a Frenet normal field. For example, any Frenet normal 
field N for the curve (cf. [36, Chapter 1]) 


(x, e7, 0), x <0 
y:[-1,1]> R, y@œ) = { (0,0,0), x=0 


(x, 0, ey, x>0 
would have to satisfy (see Fig. 5.14) 


E 0 
Nove e3 X < 
me? x> 0, 


which is impossible for a smooth map. Figure 5.14 shows a curve where four planar 
curves are stitched together in a smooth fashion, together with an attempt to define 
a Frenet normal field for this curve. 

Even if a Frenet normal field exists on an open dense set of [a, b] (which in 
general is not guaranteed), it can exhibit singularities that can be worse than the 
jump discontinuities from the previous example. For example, any Frenet normal 
field for the curve 


1 1 pe | 
Ex TI eS y IIs 0 
FLIS R3, yole (x, ex cos(.),e* sin()), x < 
(t, 0, 0), x>0. 
will have unbounded rotation speed, as is visible in Fig. 5.15. After Frenet normals 
were introduced in the middle of the nineteenth century, they quickly became a 
popular tool for studying curves. The second half of the nineteenth century saw 
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Fig. 5.14 Even when a Frenet normal field exists on an open dense set, on the whole curve there 
might be no such field 


Fig. 5.15 Away from a single point, this curve has a Frenet normal. However, its rotation speed t 
is unbounded 


the powerful appearance of Complex Analysis and Algebraic Geometry in the 
landscape of Mathematics, while Topology (and certainly Differential Topology) 
were still in their infancy. In those days it seemed natural to assume that the curves 
y under consideration were real analytic (locally representable as a power series). 
And every real analytic curve does indeed have a Frenet normal field: 


Theorem 5.19 
Ify: [a,b] > R” is real analytic, then y has a Frenet frame. 
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Proof. Without loss of generality we may assume that y has unit speed. If y 
parametrizes a piece of a straight line, then every unit normal field along y is Frenet 
and we are done. Otherwise, because of the real analyticity of y, there are only 
finitely many parameter values x1,...%m € [a,b] where y” vanishes. On each 
subinterval of [a, b] bounded by two of the points a, x1, ..., Xm, b there is a Frenet 
normal field, unique up to sign, which is obtained by setting N = fan . It is therefore 
sufficient to show that also in the neigborhood of each x; there is a Frenet normal 
field, unique up to sign. In the end, the signs can then easily be adjusted to yield 
a Frenet normal field on the whole interval [a, b]. By real analyticity, there is a 
neighborhood of x; where y can be expressed as 


CO 
v(x) = Do aa — xj)" 
k=0 
with ag € R”. Then 


(x) =J kk- Dax — xj)? 
k=2 


and there is an index £ € N such that ag = 0 for k = 2,...,€— 1 but ag Æ 0. Then 


y= (x — xj)? E+ DEK — Daa — xj =: & xja) 
k=0 


with n(x) 4 0 for all x in some neighborhood of x ;. In this neighborhood 


NO, 
|n| 


is the desired Frenet normal field. oO 


Nowadays, the standard assumption for curves is that they are smooth, i.e. 
infinitely often differentiable. Because for n > 3 not every C® curve in R” 
has a Frenet normal field, for n > 3 these fields cannot be used for studying 
global questions about smooth curves in R”. Moreover, when used in the context 
of numerical algorithms that operate on space curves, Frenet normals can cause 
unexpected behavior near curves that do not have a Frenet normal. 
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Surfaces and Riemannian Geometry 


The most simple quantity of a one-dimensional curve y : [a, b] > R” is its speed 
ly’|: [a, b] > R. The goal of this chapter is to arrive at the analogous statement for 
a two-dimensional surface in R”. Our first task will be to replace the interval [a, b] 
by a suitable domain of definition M C IR. For a surface f: M — R” the analog 
for the speed of a curve will be a Riemannian metric induced on M by f. 


6.1 Surfaces in R” 


Our investigations of curves in R” will be the guideline when we now start to 
study surfaces. For the most part we will focus on surfaces in R3. We will study 
the curvature of surfaces and the analog of the length of a curve (obviously the 
area of a surface) as well as the analog of the total squared curvature (called the 
Willmore functional). We will study the critical points of the area under variations 
with support in the interior (these surfaces are called minimal surfaces) and of the 
Willmore functional. We will prove a famous result that concerns the surface analog 
of f PK ds, the so-called Gauss-Bonnet theorem. We will investigate the analog 
(called the Euler characteristic) for the tangent winding number of a curve in R?. 

In our discussion of (non-closed) curves y in R”, y was always defined on a 
closed interval [a, b]. It would have made little difference if y would have been 
defined on the finite union of pairwise disjoint intervals. In the case of surfaces, it 
will be useful to allow for such disconnected domains. 
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Definition 6.1 


A subset M cC R? is called a connected compact domain with smooth 
boundary if 


M = Mo \ {M U... U My} 

where for each j € {0,..., k} 

Mj = 9;(D) 
is the image of the unit disk 

D:= {pe R* | |pl< 

under a diffeomorphism 

gj: D> R? 
and the M; are pairwise disjoint and contained in the interior of Mọ. A finite 
disjoint union of connected compact domains with smooth boundary is called a 


compact domain with smooth boundary (see Fig. 6.1). 


Notation Throughout the rest of the book, M will denote a compact domain with 
smooth boundary in R?. 


In order to avoid having to mention regularity constantly, we include regularity 
in the definition of a surface (see Fig. 6.2): 


M 


Fig. 6.1 A compact domain with smooth boundary 
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Go- 
SDN 


Fig. 6.2 A surface f: M > R? 


Definition 6.2 


A surface in R” is a smooth map f: M — R” whose derivative f’(p) is an 
(n x 2)-matrix of rank 2 for all p € M. 


We will denote the coordinates in IR? by u and v. Partial derivatives with respect 
to u or v will be denoted by subscripts, so for a surface f in IR” the matrix-valued 
function f’: M — R”*? is of the form 


| | 
f= tu fo 
| | 


with f,(p), fu(p) € R” linearly independent for all p € M. The following two 
definitions are special cases of the ones in Appendix A.1. 


Definition 6.3 


A map f: M — R” is called smooth if it is of the form f = flu for some 
smooth map f: U — R” where U C R? is an open set that contains M. 


It is easy to check that even at boundary points p € M the Jacobian matrix f (p) 
only depends on f, not on the specific way in which f extends f. We therefore can 
safely define f’(p) := f'(p). 


Definition 6.4 


If M,M c R? are two compact domains with smooth boundary, a bijective 
map yg: M —> M is called a diffeomorphism if both g and g~! are smooth. 
A diffeomorphism ¢ is called orientation-preserving if det y’(p) > 0 for all 
pemM. 
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cS C — 


Fig. 6.3 A reparametrization f of a surface f 


Definition 6.5 


If f: M —> R” and f: M —> R” are two surfaces, then f is called an 
(orientation-preserving) reparametrization of f if there is an (orientation- 
preserving) diffeomorphism g: M —> M such that 


f=fop, 
(see Fig. 6.3). 


As in the case of curves in the plane, it is not difficult to check that reparametriza- 
tion (as well as orientation-preserving reparametrization) defines an equivalence 
relation on the set of surfaces in R”. Although we will not formalize this, we are only 
interested in properties of surfaces that are invariant under orientation-preserving 
reparametrization, so the real objects of our study are the equivalence classes of 
surfaces under reparametrization. 


6.2 Tangent Spaces and Derivatives 
Let M c R? be a compact domain with smooth boundary and f: M > R* a 
smooth map. Then the directional derivative of f at a point p € M in the direction 


of a vector X € R? is given by 


df (p, 8) := f'(p)X. 
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This means that all these directional derivatives are encoded ina mapdf: M x 
R? > RÉ: 


Definition 6.6 


For a point p € M, a tangent vector to M at p is a pair X = (p, X) where 
X e€ R’. The set 


TM = {p} x R? 
of all these tangent vectors is called the tangent space to M at p. We make each 
T,M into a two-dimensional real vector space by defining for X = (p, X), Y = 
(p,Ê)andà eR 
X+Y=(p,X+Y) 
AX = (p, AX). 


The union TM = M x R? of all these tangent spaces is called the tangent 
bundle of M. The map 


z:TM—>M, (p, È) > p 
is called the projection map of the tangent bundle. 


One immediate benefit of this definition is a more concise notation for deriva- 
tives: 


Definition 6.7 


For a smooth map f: M — R” we define the derivative df: TM — R” of f 
by setting for X ¢ TM, X = (p, X) 


df (X) = f'(p)X, 
(see Fig. 6.4). 


The restriction of df to each tangent space T,M is a linear map from T,M to R”. 


Definition 6.8 


A smooth map X: M — TM is called a vector field if 7 o X = idm, which 
means that X (p) € TpM for all p € M. 
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Fig. 6.4 Two tangent vectors X, Y € T,M and their image under df 


If X: M > R? is a smooth map, then the assignment 
X: M >TM, X(p) = (p, X(p)) 


is a smooth vector field on M and all smooth vector fields are obtained in this way. 
Here is some convenient notation: 


Definition 6.9 


(i) The vector space of all smooth functions f: M — R is denoted by C% (M). 
(ii) The vector space of all smooth functions f: M — R” is denoted by 
C®(M, R”). 
Gii) The vector space of all smooth vector fields on M is denoted by r (T M). 


As is known from calculus class, for a smooth map f: M — R” the vector 
T (p)X €e R” can also be interpreted as the directional derivative of f at p in 
the direction of the vector X € R?. With this in mind we define the directional 
derivative of f € C (M, R”) in the direction of a vector field X € Ir (T M) by 


(dx f)(p) = dx(p) f = df (X (p)). 
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Definition 6.10 


The coordinate vector fields U, V € (TM) are defined as 


so) oe) 


The directional derivatives in the direction of U or V are just partial derivatives: 


du f = fu, dy f = fv. 


Definition 6.11 


Let M, M C R? be two compact domains with smooth boundary and ¢: M —> 
M a diffeomorphism. Then we define 


do: TM >TM, do(X) = (o), '(p)X) for X= (p. x). 


> Remark 6.12 Note that for X € T, M the vector dø (X) is an element of Typ) M, 
while for a surface f: M — R” the vector df (X) is just an element of R”, not 
an element of something like Tf(p) R”. We are relying here on the fact that in our 
situation we can naturally identify all such tangent spaces T4 R” with R” itself. This 
mild context-dependency of notation should lead to no confusion. It is very useful 
and common in Differential Geometry. 


With this notation in place, the chain rule now emerges in its most elegant form: 


Theorem 6.13 


(i) Suppose f: M —> R” is areparametrization of the surface f: M > R”, 
ie. f = f og for some diffeomorphism ọ: M —> M. Then 


df =df o dọ. 


(ii) If M, M, M c R? are compact domains with smooth boundary and 
p: M > M and ĝ: M > M are diffeomorphisms, then 


d(ọ o ©) = dọ o dẸ. 


Proof. The proof just involves spelling out our definitions and applying the ordinary 
chain rule. o 
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6.3 Riemannian Domains 


When it comes to investigating the geometry of a surface f: M — R”, the 
geometry of M as it sits in R? is completely irrelevant. Things like the length of 
a vector or the angle between vectors should be computed in the target space R” of 
f, not in R*. Accordingly, we endow each tangent space T,M with its own private 
Euclidean scalar product by defining 


()¢: U pM x TyM) > R, (X,Y) ¢ = (df (X), af (Y)). 
peM 


It is easy to check that for each p € M the restriction of (,) ¢ to T,M x TM is 
indeed a positive definite scalar product on Tp M. With respect to this scalar product, 
X € T,M is a unit vector if and only if df (X) € R” is a unit vector. 


Definition 6.14 
(,) ¢ as defined above is called the metric on M induced by f. 


> Remark 6.15 In older texts the induced metric is often called the first fundamen- 
tal form. We will not use this terminology. 


In general, objects like (, ) ¢ are interesting even when they are not induced by 
a map f: M — R”. That is, (,)¢ has the properties of a scalar product between 
two vectors, which allows us to measure lengths an angles. However, one can freely 
choose other ways to define such a metric, without explicit reference to a surface f, 
as long as the scalar product properties are satisfied. 


Definition 6.16 


Let M be a compact domain with smooth boundary in R?. Then: 
(i) A map 


(): LJ GM x TM) >R 
pEM 


is called a Riemannian metric on M if for each p € M the restriction of (, ) 
to TM x T,M is a positive definite scalar product and for any two smooth 
vector fields X, Y € T (T M) the function 

(X,Y): M>R 


is smooth. 
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(ii) M together with a Riemannian metric (,) on M is called a Riemannian 
domain. 


> Remark6.17 For brevity of the notation we will usually omit the index (-) f even 
for Riemannian metrics which are induced by some f: M —> RR. The inserted 
vectors should provide enough context to avoid confusion. 


A Riemannian metric (, ) gives rise to a function 
|-|: TM > R, |X| = /(X, X). 
The restriction of | - | to each tangent space is indeed a norm on T, M. One should 
note that the coordinate vector fields U and V are not necessarily orthonormal with 
respect to this induced metric. In fact, this is only true for special surfaces. We will 


elaborate more on this in Sect. 6.5. 


Example 6.18 
The norm corresponding to the metric (, ), induced by the inclusion map 


u M >RP, (u, v) => (‘) 
v 


satisfies 
|. |? = du? + dv’. 


The above equation should be read as a literal equality of two functions on T M. 


6.4 Linear Algebra on Riemannian Domains 


Even in the absence of a Riemannian metric, each single tangent space T,,M is a 
playing field for Linear Algebra. 


Definition 6.19 


A smooth map 
A:TM >TM 


is called an endomorphism field if its restriction to each tangent space T, M is 
a linear map 


Ap: TpM > TyM. 


Remember that TM = M x RÈ, so it is clear what we mean by a smooth 
map from TM to TM. In particular, it is clear what we mean by a smooth 
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endomorphism field. For every smooth endomorphism field A there are smooth 
functions a, b, c,d: M — R such that 


AU =a-U+c-V 
AV =b-U+d.-V. 


This means that a smooth endomorphism field basically is the same thing as a 
smooth map 


We always have the identity map as a canonical endomorphism field: 
Il: TM—>TM, 1X =X forallX «ETM. 
If A is an arbitrary smooth endomorphism field on M, for each p € M we can 


take the determinant or trace of the restriction of A to T,M and obtain smooth 
functions 


detA,trA: M— R. 


In the presence of a Riemannian metric we can define the adjoint of an endomor- 
phism field: 


Theorem 6.20 

Let (,) be a Riemannian metric on M and A a smooth endomorphism field 
on M. Then there is a unique smooth endomorphism field A* on M such that 
for all vector fields X, Y € T (T M) we have 


(AX, Y) = (X, A*Y). 


Proof. By definition of a Riemannian metric, the map 


G: M > R??, G= U) (U, a 
(V, U) (V, V) 
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is smooth and the matrix G (p) is invertible for all p € M. Now one can check that, 
given an endormorphism field A such that 


AU =aU+cV 
AV =bU +dV 


the endomorphism field A* defined by 


with 


TIE T 
ac —ı fab 
~~|=G G 
(; ‘) 6 n) 
is smooth and satisfies the desired identity. The uniqueness part of the claim is 
straightforward. o 


The endomorphism field 7 defined above is self-adjoint, which means 7* = I. The 
only structure on M we want to inherit from RR? is the notion of orientation: 


Definition 6.21 


Two vectors 
X = (p, Å) , Y = (p, Î) € TpM 


are said to form a positively oriented basis of T, M if X,Y c R?area positively 
oriented basis of R?, i.e. detg2 (X, Y) > 0. 


Each tangent space T, M of a Riemannian domain comes with its own determi- 
nant form: 
Theorem 6.22 


Let (,) be a Riemannian metric on M. Then there is a unique map 


det: (_) (T,M x TM) > R 
peM 


(continued) 
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Theorem 6.22 (continued) 
such that for every p € M the restriction 


det |7,MxT,M 
is a skew-symmetric bilinear form such that 
det(X, Y) = 1 


for every positively oriented orthonormal basis of T)M. The map det is called 
the area form of the Riemannian domain (M, (, )). 


Proof. The vector fields 


U 
(U, U) 
mS (U, U)V — (V, U)U 
© sU, UVU, UNV, V) — (U, V}? 


are orthonormal at each point p € M. Therefore, the function det we are looking 
for has to satisfy 


det(U, V) 
JU, U, XV, V) = (U,V)? 


1 = det(X, Y) = 


The skew-symmetric bilinear forms on T, M form a 1-dimensional vector space, so 
there is a unique such form det for which 


det(U, V) = JU, UV, V) — (U, VY, 
where the equality has to be read point-wise. This form already satisfies det(X, Y) = 
1. On the other hand, every positively oriented orthonormal basis of T, M is of the 


form 


X= cosa X(p) — sina Y(p) 


Y= sina X(p) + cosa Y(p) 


for some a € R. Therefore, we also have det(X oy j= 1; o 
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Theorem 6.23 


Let (,) be a Riemannian metric on M and det the area form defined in 
Theorem 6.22. Then there is a unique endomorphism field J on M such that 
for all p € M and all X, Y € T,M we have 


(JX, Y) = det(X, Y). 
In terms of the coordinate vector fields, J is given by 


_ WU, Z)V -—(V, Z)U 
JU, UV, V) — (U, Vy 


Hence J is smooth. For every positively oriented orthonormal basis of TpM we 
have 


JX=Y 
JY = —X. 


So, J operates in each tangent space Tp M as the 90°-rotation in the positive sense. 


The proof is straightforward and left to the reader. The theorem will be useful in 
several instances: 


Theorem 6.24 
For vectors X, Y, Z € TM the following identity holds: 


(X, ZYY — (Y, Z)X = det(X, Y)JZ. 


Proof. Both sides of the claimed identity are linear in X and Y and by Theorem 6.23 
the formula that we want to prove is true whenever X, Y € {U(p), V(p)}. Since 


U(p) and V (p) form a basis of T, M, the claimed identity then is true for all X, Y € 
T M. o 
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6.5 Isometric surfaces 


Definition 6.25 


Two surfaces f, f: M —> R” are called isometric if they induce the same 
Riemannian metric {, ) on M. 


Note that f and f are isometric if and only if 


(Jus fu) = Gus Fu 
(Fus fo) = (has fo) 
(Fo, fo) = Ris fo). 


The physical intuition concerning isometries is as follows: The deformation 
of the surface f to the surface f involves only bending, without any intrinsic 
deformation such as stretching within the surface. In the nineteenth century, 
geometers liked to demonstrate this using leather patches. By methods known for 
example to shoemakers, a leather patch can be brought into any initial shape. After 
the initial preparation, the leather can still be bent easily, but it will not allow 
stretching. In Fig. 6.5, the initial shape f of the patch is a planar ring. This ring 
can easily be fitted to a cone, assuming a shape f. 

It is clear that one can slide the ring freely around on the cone, in all directions. 
Few surfaces have the property that one can take a piece of the surface and slide it 
without distortion or stretching around on the surface. For example, the leather patch 
on the surface in Fig. 6.6 is clearly stuck in place. One famous surface on which 
such a patch can freely slide, already known in the nineteenth century and a popular 
tool for the leather demonstration, is the pseudosphere, that can be parametrized as 
follows: 


Fig. 6.5 The maps f and f are isometric, as the flat leather patch f can be placed onto a cone 
without tearing or stretching. It can even slide freely on the cone 


6.5 Isometric surfaces 101 


Fig. 6.6 The leather patch on this dodecahedron is stuck in place 


Fig. 6.7 A leather patch fitted to the pseudosphere (at first in one place) is able to slide around 
freely and even rotate freely 


Suppose that M C {(u, v) € R? | v > 1} and define 


cos(u) 


U 
sin(u) 


v 
log (Vo? =T + v) ie eal 


v 


f: M >R, fu, v) = 


Choose à > 1 and u € R and define f: M > R? by flu, v) = f (àu + u, àv). 
We leave it to the reader to verify 


(Jus tv) = Fes fo) =0 
ae ~ z 1 
(Jus fu) = (fos fo) = (fue fu) = (fo, fo) = z 


and that therefore f and f are isometric (see Fig. 6.7). 
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Here is another example, which will also be of interest later when we study 
minimal surfaces: Suppose that M C {(u,v) € R? | u > 0}. Letk,£ € Z be 
two integers with k + £ Æ —1. Then define the Enneper surface f: M — R? by 


u2k+! COCE Dv) _ utt! ea a 


fu, v) = u2k+! snk + y2etl sin GH) 
24 AT ep costet Do) 
KFI 


Again, we leave it to the reader to verify that 


(fas fo) =0 
| fuu, v)| = u% 


|fo(u, v)| = yokti 4. yer! 


u” 


and that for arbitrary à € R the analogous formulas also hold for 
f:M >R 
fu, v) = fu, v +A). 


This means that also here f and f are isometric, so a leather patch has at least 
one degree of freedom to slide on the surface without stretching (see Fig. 6.8). 


Fig. 6.8 A leather patch fitted to the Enneper surface is able to slide around, but only in one 
direction. It cannot rotate 
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Integration and Stokes’ Theorem 


For a curve y: [a, b] > R”, global quantities like the length or the bending energy 
were defined as integrals over arclength of certain functions on [a, b]. Before we 
can define similar quantities for a surface f: M — R”, for example, the area of 
f., we have to find a way to integrate functions g: M — R in a geometrically 
meaningful fashion. We will do this in terms of the area 2-form det of the metric 
induced by f. This means that first we have to develop the theory of differential 
forms on two-dimensional domains M, including the theorem of Stokes. 


7.1 Integration on Surfaces 


Let M, M C R? be two compact domains with smooth boundary, g: M —> M an 
orientation-preserving diffeomorphism and g: M — R. a smooth function. By the 
transformation formula for integrals, we have 


f sovag = | g. 
mM M 


Therefore, if one were to just use f m & as the definition for an integral of a function 
g over the surface f, this integral would not be invariant under reparametrization 
of f. On the other hand, we are now going to convince ourselves that it is perfectly 
possible to define the integral of an object like the area form det, as it was introduced 
in Theorem 6.22: 


Definition 7.1 
Let M be a compact domain with smooth boundary in R?. Then a map 


o: |_](TpM x TM) >R 
peM 
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is called a 2-form on M if for each p € M the restriction of o to TpM x TpM is 
a skew-symmetric bilinear form and the function 


o(U,V): M >R 
is smooth. 


We denote the set of all 2-forms on M by Q?(M). As a linear subspace of the 
vector space of all real-valued functions on some set, also Q?(M) is a real vector 
space. It is also clear how to define the product of a 2-form o with a function g € 
C™(M). One can say that 2-forms are similar to Riemannian metrics, only skew 
symmetric instead of symmetric and without any non-degeneracy assumptions. 


> Remark 7.2 A “model example” of a 2-form which we have already encountered 
is the well-known det € Q?7(M). 


2-forms are transported under diffeomorphisms by demanding that the trans- 
ported form applied to the transported tangent vectors yields the same value as 


before: 


Definition 7.3 


Let M, M C R? be two compact domains with smooth boundary, g : M—>Ma 
smooth map and o a 2-form on M. Then we define the pull-back of o under o 
as the 2-form y*o on M that for p € M and X, Y € T,M is given by 
(y*0)(X, Y) := o (dg(X), dp(Y)). 
Theorem 7.4 
In the situation of Definition 7.3, the map 
Q? (M) > 27(M), o B gto 


is linear and for g € C (M) satisfies 


g“ (go) = (go y)(y"*0). 
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The integral over M of a 2-formo € Q?(M) is defined as follows: 


Definition 7.5 


The integral of a 2-form o on a compact domain M with smooth boundary in 


R? is defined as 
1 o = ao (U,V) 
M M 


where U and V are the two vector fields on M introduced in Definition 6.10. 


The above definition is useful because f m © iS invariant under pull-back of o by 


an orientation-preserving diffeomorphism g: M —> M. 


Theorem 7.6 
Let M, M C R? be two compact domains with smooth boundary, o: M —> M 
an orientation-preserving diffeomorphism and o € Q7(M) a 2-form. Then 


[eo=] o 
M M 


Proof. Let us write 


which means 
dy(U) =aUog+cVog 
dy(V)=bUcg+dVog. 


Therefore, by the skew symmetry of o and the transformation formula, 


i gro= fe (y*o)(U, V) 
M M 


= j) OE, 


M 


= f ad- voy ww, V)og 
M 
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=f o (U, V) 
M 


afo 


In the context of surfaces f: M —> R”, we will never integrate functions g € 
C™(M) directly, but instead we will first make g into a 2-form by multiplying it 
with the area form det of the induced metric. Then we can be sure that 


/ g det 
M 


is a quantity that will stay the same if we reparametrize f as f = f o ọ (and, of 
course, simultaneously change g to g o g). Theorem 7.6 above makes it possible to 
define the area of a Riemannian domain in such a way that it does not change under 
isometries: 


oO 


Definition 7.7 


The area of a Riemannian domain (M, (, )) is defined as 


/ det 
M 


where det is the area form of (, ). 


7.2 Integration Over Curves 


In order to adequately deal with surfaces, we have found it necessary to add tangent 
spaces, Riemannian metrics and 2-forms to our toolbox. Let us investigate whether 
some of these notions might be useful already in the context of curves. For a curve 
y: [a, b] > R, the analog of the domain M of a surface f: M — R” is the interval 
[a, b]. The tangent bundle of [a, b] is 


T[a, b] = [a,b] x R 


and the tangent space at p € [a, b] is {p} x R. The analog of the vector fields U, V 
on M is the single vector field X € T ([a, b]) defined as 


X(p)= (p, 1). 
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Fig. 7.1 A 1-form œ € Q! (M) can be thought of as a smoothly varying ruler which “measures” 
a vector field X € I'(M). The spacing of the ruler-lines indicates the “strength” of w—the closer 
the spacing, the stronger is w 


The objects that can naturally be integrated over curves are the so-called 1-forms. 
We will need 1-forms also on planar domains (Fig. 7.1), so we take the opportunity 
to define also those. 


Definition 7.8 


Let [a,b] be a closed interval and M Cc R? a planar domain with smooth 
boundary. Smoothness of maps defined on T[a, b] or TM is to be understood 
in the sense of Definition A.1. Then 


(i) A smooth map æ: T[a, b] > R is called a 1-form if its restriction to each 
tangent space Tp[a, b] is linear. The space of all 1-forms on [a, b] is denoted 
by 2! (fa, bI). 

Gi) A smooth map w: TM — R is called a 1-form if its restriction to each 
tangent space 7, M is linear. The space of all 1-forms on M is denoted by 
Q! (M). 


A general theory of m-forms on domains in Ré is beyond the scope of this book, 
SO we just collect some special cases that we need: 


Definition 7.9 
Let [a, b] C R be a closed interval and M C R? a planar domain with smooth 


boundary. 


(i) If w is a 1-form on [a, b] and g: [a, b] > [a, b] is a smooth map, then we 
define the pull-back of w under ¢ as the 1-form g*w € Q! ({a, b]) which is 
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for all Y € r(T[ã, b]) given by 
(y*@)(Y) = o(dg(Y)). 
Gi) If œw is a 1-form on M and ¢: M — M is a smooth map, then we define 
the pull-back of w under ¢ as the 1-form g*w € Q'(M) which is for all 
Y e T(T M) given by 
(y*@)(Y) = w(dg(Y)). 
(iii) If æ is a 1-form on M and y: [a,b] —> M is a smooth map, then we define 
the pull-back of w under y as the 1-form y*w € Q! ({a, b]) which is for all 
Y eT (7[a, b]) given by 


(y*o)(Y) = o(dy(¥)). 


Definition 7.10 


For œ € &!([a, b]) we define the integral of a 1-form w over [a, b] as 


b 
f o= f w(X). 
[a,b] a 


Theorem 7.11 
Let g: [a,b] — [a,b] be an orientation-preserving diffeomorphism, i.e. a 
bijective smooth map with gy! > 0. Then 


/ g“w = / w 
[a,b] [a,b] 
Proof. By the substitution rule and dg(X) = ọ' - X oy we have 


b 
Í , go = f (y*@)(X) 


b ~~ 
7 I oid 


i 
=f wl -X og) 
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b 
al g'-w(X)og 


b 


Theorem 7.12 

Let M C R? be a compact domain with smooth boundary and w € Q'(M) 
a l-form. Let y: [a, bD] > M bea reparametrization of a smooth map 
y: [a,b] > M, soy = y o ọ for an orientation-preserving diffeomorphism 


gy: [ã, b] > [a, b]. Then 
/ y*o =|! y*o. 
[ã,b] [a,b] 


Proof. By the chain rule, we have dy = dy o dọ and therefore 
“o = ¢*(y*o). 


Therefore, Theorem 7.11 gives us 


f g w= : ve = | yo. 
[a,b] [a,b] [a,b] 


oO 


As a consequence, we can define the integral of a 1-form w on M over a curve in M 
in a way that is invariant under reparametrization: 
Definition 7.13 


Let M C R? be a compact domain with smooth boundary and w € Q!(M). Let 
y: [a,b] — M be a curve. Then we define 


fox] y*o. 
Y [a,b] 
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In the context of a regular curve y: [a, b] > R”, what is the analog of the area 
form det of a surface f: M > R”? 
Definition 7.14 


The arclength 1-form ds € Q!({a, b]) of a curve y : [a, b]  R” is defined as 
ds(X) = |dy(X)| = |y". 


If we define the arclength function s as in Definition 1.13, the arclength 1-form ds 
is indeed the derivative of s, which explains the notation. Moreover, if we interpret 
the left-hand side according to Definition 7.13 and the right-hand side according to 
Definition 1.14, for a function g: [a, b] > R we have 


b 
fl eds= | gds. 
[a,b] a 


7.3 Stokes’ Theorem 


When dealing with curves, we frequently used the fundamental theorem of calculus, 
for example in the form of integration by parts. Also in surface theory we would 
no get very far without the surface analog of this theorem, which is the so-called 
Stokes’ theorem. 

Let M C R? be a compact connected domain with smooth boundary. The 
boundary 0M of M can be parametrized by a finite collection of n closed curves 


Yj: laj, bj] => R? 


where j € {1,...n}. We assume that each y; is oriented in such a way that for any 
vector Y € R? which at yj; (x) points out of M, we have 


det(Y, yj(x)) > 0. 


Given a 1-form œw € Q!(M), we make use of Definition 7.13 in order to define the 
integral of w over 0M: 
fue 
aM 


n 


Djo 


j=1 
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Theorem 7.15 (Stokes Theorem) 
Let M C R? be a compact domain with smooth boundary and w € Q'(M). 
Then there is a unique 2-form dw € Q7(M) such that for all subdomains 


M C M we have 
f w=) w. 
M aM 


In fact, dw is the unique 2-form on M that satisfies 


dæ(U, V) = @(V)u —@(U)y. 


Proof. M could be an arbitrarily small disk around an arbitrary point p in the 
interior of M, so there can be at most one 2-form o with the property that for all 


subdomains M 
Í. pa l p 
M aM 


This proves the uniqueness part of the claim. If we write 


we have 
dyj(X) =ajUoyj + Bj Voy; 
and therefore 
w(dyj(X)) =aj aU) oy; + Bj a(V) 0 y;. 
Let us define o € 2?(M) as the unique 2-form for which 
o(U,V) =a(V)y — OU )y. 


Now we apply Green’s theorem from vector calculus to the map 


Y: M > RÊ, Y = e 
w(V) 
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and obtain 


pal o(U, V) 
M M 


=f (@(V)u — @(U)v) 
M 


n bj 
=o] @joU)oyj+BjoV)oyj) 


fare 


n bj 
=> f vay) 
jai" 4i 


We can apply this argument also to any subdomain M c M, which proves the 
existence part of the claim. o 


Theorem 7.16 
If 9: M —> M is a smooth map and w € Q! (M), then 


g“ (dw) = d(g* a). 


Proof. The proof is easy if ọ is an orientation-preserving diffeomorphism: If Mc 
M is any subdomain, then by Theorems 7.6, 7.11 and 7.15 we have 


[ede = f : ao = f : w= | oe 
M (M) ayg(M) aM 


By the uniqueness part of Theorem 7.15 we then must have y* (dw) = d(g*a). 
Unfortunately, here we only assume that g is a smooth map, so we have to rely 
on the coordinate formula provided in Theorem 7.15. We use the notation from the 
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proof of Theorem 7.6 together with the equalities a; = bg and cz = dj that follow 
from the commutativity of partial derivatives of the component functions of ø to 
compute: 
d(y*w)(U, V) 
= w(dg(V))z — o(dgO))i 
=(b-wU)og+d-w(V)og)i—(a-oU)ogt+c-#(V) og)s 
= b(law(U) og + cw), 09) +d(a-aV),o~tc-aV) og 
—a(b-oU),og+d- aU), 0g)—c(b:o(V) op +d: oV) 09) 
= (ad — bc) (@(V)u — œ (U )v) 0g 
= (ad — bc)dæ(U, V) o 
=do(a-Uog+c:-Vog,b-Uogt+d-Vog) 
= dw(dg(U), dp(V)) 
= (y*dw)(U, V), 


which proves the claim. o 
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Curvature 


From this chapter on we will focus attention on surfaces f: M —> R3. The most 
fundamental tool for analysing such a surface is its unit normal field N: M —> S? 
which is a map to the unit sphere S? C RÌ. The derivative of N reveals information 
about the curvature of f. In particular, the area covered by N on S$? provides us with 
a geometric interpretation of the so-called Gaussian curvature of f. 


8.1 Unit Normal of a Surface in R3 
Most of the material in the Chaps. 6 and 7 was concerned with the intrinsic geometry 


of Riemannian domains or with surfaces in R”. From now on we will focus on 
surfaces f: M > R?. 


Definition 8.1 


Let M C R? be a domain with smooth boundary and f: M —> R? a surface. 
Then there is a unique smooth map N: M —> R? with (N, N) = 1 such that 


(i) For all p € M and all X € T,M we have 
(N(p), df (X)) = 0. 
(ii) For all p € M and every positively oriented basis X, Y of T,M we have 


det(N (p), df (X), df (Y)) > 0. 


N is called the unit normal of f (see Fig. 8.1). 
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Fig. 8.1 The normal vector N (p) of a surface f at a point p 


In terms of the coordinate vector fields U and V we can express N as 


_ fu X fo 


N= ; 
|fu x fol 


Theorem 8.2 
Forall p € M andall X,Y € T,M we have 


df (JX) = N(p) x df (X) 
dets (X, Y) = det(N(p), df (X), df (Y)). 


For the area of f we get 


apel dy = f dets (U, v=f det(N, ja [fu X fol. 
M M M M 


Similar as for a surface f: M > R3, we can consider the derivative d N of the unit 
normal N: M -> R?. In the case of plane curves the derivative of the normal N 
gave us the curvature « via the equation 


/ 


N'=ky’. 
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In order to find the analogous equation for surfaces, let us consider a vector field 
X e T(TM) and take the derivative in the direction of X of the equation | = 
(N, N): 


0 = dx(N, N) = 2(dxN, N). 


This means that for all X € TpM the vector dN(X) lies in the image of the 
restriction of df to Tp M. Therefore, there is a vector Y € T ,M such that d N (X) = 
df (Y). Obviously, the dependence of Y on X is linear, so there is a linear map 
Ap: TpM — T,M such that for all X € Tp M we have 


dN(X) = df (AX). 
We leave it to the reader to check that A is a smooth endomorphism field on M. 


Definition 8.3 


The smooth endomorphism field A is called the shape operator of f. 


Theorem 8.4 
The shape operator A is a self-adjoint endomorphism field with respect to the 
induced metric, i.e. for all X, Y € T (T M) we have 


(AX, Y) = (X, AY). 


Proof. Since at each point p € M the two vectors U (p), V(p) form a basis of Tp M, 
it is sufficient to prove the theorem in the special case X = U, Y = V. Using the 
fact that 


(N,df(U)) = (N,df(V)) =0 
we obtain 


(AU, V) = (df (AU), df(V)) 
= (dN), df(V)) 
= dy(N, df(V)) — (N, dudf (V)) 
= -(N, fou) 
= -(N, fuv) 
= dy(N,df(U)) — (N, dydf(U)) 
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= (dN(V), df (U)) 
= (df (AV), df (U)) 
= (AV, U), 


where we used that the partial derivatives commute, i.e. fuv = fvu- oO 


8.2 Curvature of a Surface 
The shape operator A of a surface f: M — R? captures all the information about 


how the surface is curved. In fact it measures deviation from being planar: 


Theorem 8.5 

Let M C R? be a connected compact domain with smooth boundary and 
f: M > R? a surface with shape operator A. Then A vanishes identically if 
and only if there is a plane E C R? with f (M) C E. 


Proof. If f(M) C E with 
E = {pe R? | (Ñ, p) = c} 
for some unit vector Ñ € R? and c € R, then 


(Ñ, df (X)) = dx (Ñ, f) =0 


for all X € TM, so the unit normal of f satisfies N (p) = +N for all peM.īin 
particular, dN = 0 and therefore A = 0. 

Conversely, by the connectedness of M, A = 0 implies that N is constant, i.e. 
N(p) = Ñ for some Ñ € R? and all p € M. Then d(N, f) = O and (by the 
connectedness of M) there is c € R such that (Ñ, f(p)) =c forall p e M. o 


At a given point, a surface can be curved by a different amount in different 
directions. We call a vector X € TM a direction if (X, X) = 1. 


Definition 8.6 


For a direction X € TM we define the directional curvature «x (X) of f in the 
direction of X as 


K(X) := (AX, X). 
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Fig. 8.2 For every unit vector X € T,M a surface f: M > R? has a different directional 
curvature 


If X1, X2 is an orthonormal basis of T, M then we can parametrize all unit vectors 
in T,,M as 


X(0) = cos0 Xı + sin@ X2. 


Figure 8.2 contains a plot of the function 0 > «(X (0)). 
By Theorem 8.4, for all p € M the linear map 


Ap := Alr,M: T,M > TM 


is self-adjoint, so there is an orthonormal basis X1, X2 in TpM such that X; and X2 
are eigenvectors of Ap: 

AX, =k (p)X1 

AX? = k2(p)X2. 
If we assume x; > «2 the eigenvalue functions «1, «2: M — R are well-defined 
and continuous. They arise from solving the characteristic equation of A p, in which 
a square root is involved. This means that in general (if there are points where «1 (p) 
and «2(p) coincide) they are not smooth functions. 

Definition 8.7 


For p € M the numbers «ı(p) and «2(p) are called the principal curvatures 
of f at p. A vector X € TpM with (X, X) = 1 is called a principal direction 
corresponding to the principal curvature « j if 


AX =«;(p)X. 
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If we parametrize directions X (0) at p based on principal directions X1, X2 as 
above we obtain 
K(0) = (A(cos 6 X; + sin@ X2), cos@ X; + sin@ X2) 
= k\(p) cos? 0 + K2 (p) sin? 0 


= k1(p) + k2(p) i k1(p) — k2(p) 


5 5 cos(20). 


Definition 8.8 


The mean value 


1 20 
H(p) := =f k(0)d0 


is called the mean curvature of f at the point p. 


We have 


1 
o= Kı (p) Ha = st(Ap), 


so the function H: M — R is smooth. 


Definition 8.9 


The smooth function 
K: M > R, K(p) = det Ap = kı (p)k2(p) 
is called the Gaussian curvature of f. 


If K(p) > O then the directional curvatures at p are either all positive or all 
negative. In the first case, the surface looks convex when viewed from “outside” 
(when we think of N as pointing “outward”). Otherwise it looks concave. Figure 8.3 
shows surfaces whose Gaussian curvature is positive everywhere on M. 

If K(p) < 0 Then the surface bends towards N (p) is some directions and away 
from N in other directions. Figure 8.4 shows surfaces whose Gaussian curvature is 
negative everywhere on M. 
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A 


Fig. 8.3 Three surfaces with positive Gaussian curvature 


Fig. 8.4 Three surfaces with negative Gaussian curvature 


Points where the principal curvatures coincide (and therefore all directions are 
principal directions) are special and we give them a name: 

Definition 8.10 

A point p € M is called an umbilic point of the surface f if at p the surface has 

the same curvature in all directions, i.e. for all directions X € T,M we have 


k(X) = H(p). 


The most interesting theorems in Differential Geometry lead from local assump- 
tions (curvature properties at each given point) to conclusions about global shape. 
Here is our first theorem of this kind in the context of surfaces: 
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Definition 8.11 
A subset S$ C R? of the form 


S={pe€ R? | (p—m,p—m) =7’} 


with m € R? andr > 0 is called a round sphere. 


Theorem 8.12 (Umbilic Point Theorem) 
Let M C R? be a connected compact domain with smooth boundary and 
f: M —> R? a surface. Then the following are equivalent: 


(i) All points p € M are umbilic points. 


(ii) Either f(M) C E for some plane E C R? or f (M) C S for some round 
sphere 


s={peR | (P-m,p-m) =°}. 


with center m and radius r > 0. 


Proof. If f (M) is contained in a plane, we already know that A = 0 and therefore 
all points are umbilic points. If f (M) is contained in a round sphere, then there is a 
point m € R? and a radius r > 0 such that 


(f-m, f —m) =r’. 


Clearly then, f — m Æ 0 for all p e M. Differentiating the above equation reveals 
that for all p € M and all X € T,M we have 


(df (X), f — m) = 0. 
Therefore, at each p € M the unit normal of f must be given by 


1 
N(p) = +-(f (p) — m). 


r 


By the connectedness of M this implies 


1 
N = +-(f —m) 


a 
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and therefore all points are umbilic points: 
dN= +i f. 
Conversely, assume that all points p € M are umbilic points of f. Then 
Ay fu + H fuv = Nuv = Now = Hu fo + H fou 
and therefore Hy fa — Hu fy = 0. By the connectedness of M, this means that H 


is constant. In the case H = 0 we have A = 0 and by Theorem 8.5 we know that 
Ff (M) is contained in a plane. Otherwise, there is a constant r > 0 such that 


1 
H=+-. 
r 
The function 
m: M > RŽ, m(p) = f(p) £rN(p) 


then satifies dm = 0 and, by the connectedness of M, must be constant. This means 
that f (M) lies on a sphere around m with radius r. o 


8.3 Area of Maps Into the Plane or the Sphere 


Recall the second formula from Theorem 8.2: The area form det of a surface 
f: M — R? with unit normal N is given on X, Y € T M by 


det ¢(X, Y) = det(N (p), df (X), df (Y)). 


There are situations where we know what N should be, even if f is not a surface 
but just a smooth map whose derivative dp f : T, M > R? might fail to have a two- 
dimensional image for some p € M: Define the Euclidean plane £? as the subset 
of IR? where the third component is zero. Then at any point p € E? we consider the 
third basis vector e3 as the unit normal vector of E? at p. Define the unit two-sphere 
S? as the set of all p € R? with |p| = 1. Then at any point p € S* we consider p 
itself as the unit normal vector of S? at p. 


Definition 8.13 


Let M C RÊ be a compact domain with smooth boundary and g a smooth map 
defined on M with values in either E* or S?. Then we define the covered area 
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form o, € Q?(M) on X,Y € T,M as follows: 
(i) For a smooth map g: M —> E? we define 

Og(X, Y) = det (e3, dg(X), dg(Y)). 
(ii) For a smooth map g: M —> S? we define 


Og(X, Y) = det (g(p), dg(X), dg(Y)). 


If we identify Æ? with IR? in the obvious way and use the standard determinant 
det on RÊ, the first part of the above definition becomes: 


Definition 8.14 


Let M C RŽ be a compact domain with smooth boundary and g: M —> R? a 
smooth map. Then we define the covered area form og € Q? (M yon X,Y € 
T,M as 

P 


0g(X, Y) = det(dg (X), dg (Y )). 


Definition 8.15 


Let M cC R? be a compact domain with smooth boundary and g a smooth 
map defined on M with values in either E”, S* or R?. Then we define the area 


covered by a map g as 
| 
M 


Given a smooth map g: M — R? from the unit disk M into R?, we obtain a 
loop y : [0, 277] — IR? defined by 


ros) 


Conversely, every loop y : [0, 277] —> IR? arises in this way: 


Theorem 8.16 

Let y: R — R? be a loop and M C R°? the unit disk. Then there is a smooth 
map g: M — R? such that the Fig. 8.5 becomes a commutative diagram, i.e. 
Y=sgoS. 
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Fig. 8.5 By Theorem 8.17, 
the area covered by the map g 
equals the sector area of the 
boundary loop y, i.e. the blue 
region is counted positively, 
the orange region is counted 
negatively and the region with 
mixed color is not counted 

at all 


Proof. Let gy: [0, 1] ~ R be a smooth function such that 
gl) = 1 


g(x) =0 for x< 7 


Then we can define g: M — IR? as the unique map such that for all r € [0, 1] and 


all t € R we have 
t 
g (: ee )) = g(r)y(t). 
sin t 


Theorem 8.17 

Let y: [0,27] — R? be a loop, M C R? the unit disk and g: M —> RR? 
any smooth map such that Fig. 8.5 is a commutative diagram. Then the area 
covered by g equals the sector area of y. 


Proof. Define a 1-form w € Q!(M) by setting for X € TpM 


1 
o(X) = z det(ge (p), dg(X)). 
Then 


2dw(U, V) = dy det(g, dg(V)) — dy det(g, dg(U)) 
= det(dy g, dv g) + det(g, dudy 8) 
— det(dy g, dyg) — det(g, dydy 8) 
= 20, (U, V). 
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Our claim now follows from Stokes Theorem. With 


-n (cost 
PUNS ag 


we obtain 


1 2n 1 27 
> f det(y, y’) = J det(g o Ÿ, (g'o fy’) 
0 0 


oO 


By Definition 8.13, we obtain a similar interpretation for the area covered by a 
map g: M — S?. For us, the most important case is g = N where N is the unit 


normal of a surface f: M —> R? (see Fig. 8.6): 
Fig. 8.6 If the Gaussian curvature K is positive in some subregion M C M, the normal map M 


will be orientation-preserving in M (left). If K is negative, N will be orientation-reversing (middle). 


On the right we see a situation where K changes sign in M 
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Theorem 8.18 
Let f: M — R? be a surface with unit normal N and Gaussian curvature 
K. Then the covered area form of N is 


on = K det. 


Proof. For vector fields X, Y € r (T M) we have 


on (X, Y) = det(N, dN(X), dN(Y)) 
= det(N, df (AX), df (AY)) 
= det (AX, AY) 
= det A dets (X, Y) = K dets (X, Y). 
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Levi-Civita Connection 


The Levi-Civita connection of a surface f: M —> R? provides a geometrically 
meaningful way to take directional derivatives of a vector field Y on M. Based on 
the Levi-Civita connection we derive two important equations that are satisfied by 
the curvature of a surface in R°: the Gauss equation and the Codazzi equation. 


9.1 Derivatives of Vector Fields 
Let f: M — R? be a regular surface in R? with unit normal N: M — Rĉ. Let Y € 
T (T M) be a smooth vector field on M, p € M and X € Tp M. Then, differentiating 
0 = (N,df(Y)) in the direction of X, we obtain 

0 = dx(N,df(Y)) = (df (AX), df (Y )) + (N, dxd f (Y)). 
Fixing a point p € M, we can decompose every vector v € R” uniquely as 

v=AN +df(Z) 

for some à € R and some Z € TpM. The vector AN is called the normal part 
of v and df (Z) is called the tangential part of v. In our case, the normal part of 
dxdf (Y) equals —(AX, Y)N. The tangential part is of the form df(Z) for some 
vector Z € T,M that we denote by (VY)(X) or also by VxY. This gives us 


dxdf (Y) = —(AX, Y)N + df (VxY). 


We leave it to the reader to show that the map Y +> VY is linear. 
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Fig. 9.1 Two vector fields X, Y € r(T M) 


Definition 9.1 
The linear map V: F(TM) — T (End T M) that assigns to a vector field Y the 
endomorphism field VY is called the Levi-Civita connection of f. 


VxY can be interpreted as the directional derivative of the vector field Y in the 
direction of X (Fig. 9.1). 
Here are some useful properties of the Levi-Civita connection: 


Theorem 9.2 
Let X,Y, Z be vector fields on M and à: M — R. a smooth function. Then 


(i) Vx(AY) = (dxA)¥ +AVxY 
(ii) dx{Y, Z) = (VxY, Z) + (Y, Vx Z) 


Proof. Equation (i) is left as an exercise. For equation (ii) we have 


dx (Y, Z) = dx (df (Y), df (Z)) 
= (dxdf (Y), df(Z)) + (df (Y), dxdf(Z)) 
= (df (VxY), df(Z)) + (df (Y), df (Vx Z)) 
= (VxY, Z) + (Y, VxZ). 
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We can also use V to define directional derivatives of endomorphism fields (as 
defined in Definition 6.19): 


Theorem 9.3 

If B is a smooth endomorphism field on M and X € T (T M) is a vector field, 
then there is a unique smooth endomorphism field (Vx B) on M such that for 
all Y € I(T M) the following Leibniz rule holds: 


Vx (BY) = (VxB)Y + B(VxY). 


Proof. Define C as the unique endomorphism field on M for which 


CU = Vx(BU) — BVxU 
CV = Vx(BV) — BVyV. 
Clearly, if the endomorphism field Vy B exists, it has to be equal to C. This proves 
the uniqueness part of the theorem. To prove existence, we show that C has the 
property we claim for Vx B. Let us write Y as a linear combination of U and V: 
Y=aU + DV. 
Then, by Theorem 9.2, 


Vx (BY) = Vx(aBU + bBV) 
= (dya) BU + aVx(BU) + (dxb)BV + bVx(BV) 
= (dya)BU +a(CU + BVxU) + (dxb)BV + b(CV + BVxV) 
= CY + BVxY. 


A smooth endomorphism field B for which VB = 0 is said to be parallel. 
Theorem 9.4 
Let X,Y, Z be vector fields on M and à: M — R a smooth function. Then 


(i) Vx(JY) = JVxY 
(ii) dy det(Y, Z) = det(VxY, Z) + det(Y, Vx Z). 
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Proof. For Equation (i) we choose a positively oriented orthonormal vector fields 
X,Y c T(TM). Then JX = Y and JY = —X and any vector field W € (TM) 
can be written as W = BX + êY. Then for Z €e F(T M) 


Vz(JW) = Vz(BY — 6X) = (dzB)Y + BVzY — (dz8)X — 6VzX 
where we use equation (ii) of Theorem 9.2 for the last equality. Moreover, 


JIVzW = J ((dzB)X + BVzX + (dzô)Y + 6VzY) 
= (dzB)Y + BJVzX — (dzô)X + 6JVzY. 
We define w(Z) := (VzX, Y) and conclude that 
Vz(JW) = (dzB)Y — Bw(Z)X — (dz5)X —da(Z)Y = JVzW. 
In physics, the 1-form w is known as the angular velocity of the orthonormal basis 
Xp, Yp € T M in the direction Zp € T,M. 


Now equation (ii) follows: 


dx det(Y, Z) = dx (JY, Z) 
= (Vx(JY), Z) + (JY, Vx Z) 
= (J(VxY), Z) + (JY, Vx Z) 
= det(VxY, Z) + det(Y, Vx Z). 


oO 


> Remark 9.5 Note that the equations from 9.2 and the second equation from 9.4 
have the flavor of a Leibniz rule. 


Theorem 9.6 
For the coordinate vector fields we have 


VuV = WU. 
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Proof. This can be seen by looking at the tangential component of 


— (AU, VIN +df(VuV) = fou = fuv = —(AV, U)N + df (Vy). 


9.2 Equations of Gauss and Codazzi 


The derivative VA (defined in Theorem 9.3) of the shape operator A of a surface 
f: M — R has an important symmetry property: 


Theorem 9.7 
For all vector fields X, Y € T (T M) the Codazzi equation holds: 


(VxA)Y = (Vy A)X. 


Proof. We can write X and Y as linear combinations (with functions as coefficients) 
of U and V. If we expand both sides of the equation in question accordingly, we see 
that it is sufficient to consider the special case X = U and Y = V. In this case, our 
claim follows from the fact that partial derivatives of N commute: Equality for the 
normal part of 
— (AU, AV)N + df (Vu(AV)) = dydf (AV) 

= dydy(N) 

= dydy (N) 

= dydf (AU) 

= —(AV, AU)N + df (Vy(AU)) 


is automatically satisfied, while the tangential part gives us what we want to prove. 
oO 


There is another important relation between the shape operator A and the Levi- 
Civita connection V, the so-called Gauss equation: 


Ifh: M — RÝ is a smooth function, then the partial derivatives of h commute, i.e. 


dydyh — dydyh = 0. 
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For vector fields this is not true, and the failure of “partial derivatives” of vector 
fields to commute is determined by the Gaussian curvature of f: 


Theorem 9.8 
For any vector field Z € T (T M) the Gauss equation holds: 
Vu Vy Z — Vy Vy Z = —K det(U, V)JZ 


where K = det A is the Gaussian curvature of f. 


Proof. Collecting only the terms that are orthogonal to N in 


dy(—(AV, Z)N + df (Vv Z)) = dydydf (Z) 
= dydydf (Z) 
= dy (— (AU, Z)N + df (Vu Z)) 


we obtain 
— (AV, Z)AU + VyVyZ = —(AU, Z)AV + VyVuZ. 
Substituting in Theorem 6.24 AU for X, AV for Y and using 
det(AU, AV) = det A det(U, V). 


we arrive at the equality this we wanted to prove. o 


9.3  Theorema Egregium 


The following theorem is due to Gauss. He called it the “Theorema Egregium”, 
which means “most excellent theorem”. 


Theorem 9.9 (Theorema Egregium) 
Suppose that the surfaces f, f: M > R? induce the same Riemannian metric 
on M. Then f and f have the same Gaussian curvature K: M > R. 
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Proof. By the Gauss equation (Theorem 9.8), it is sufficient to prove that if f and 
f induce the same Riemannian metric on M, they also induce the same Levi-Civita 
connection. This in turn follows from Theorem 9.10. oO 


Theorem 9.10 
Suppose that the surfaces f, f: M —> R? induce the same Riemannian metric 
on M. Then the Levi-Civita connections of f and f are identical. 


Proof. We show that the Levi-Civita connection V induced on M by f is already 
completely determined by the induced metric (, )¢. By Theorem 9.6 and the second 
equation of Theorem 9.2 


1 


(VuU, V) = du (U, V) — (U, Vu V) 
= dy (U, V} — (U, Vy U) 


1 


(Vu V, U) = (VyU, U) 


1 
-dy (U, U 
5 v( ) 


1 
Hence VyU and VyV = VyU are completely determined, as well as (by a 
similar calculation) Vy V. Therefore, VU and VV are completely determined by 
the knowledge of (, )¢ alone. By the first equation of Theorem 9.2, then also VY is 
determined for an arbitrary vector field Y = bU + boV. oO 


Figure 9.2 reveals the reason why the leather patch on the smoothed dodecahedron 
in Fig. 6.6 was stuck. 
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Fig. 9.2 By the Theorema 
Egregium, any isometric 
motion of the patch has to 
preserve Gaussian curvature 
(indicated by color), so the 
patch cannot slide 
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Total Gaussian Curvature 1 0 


If we know a plane curve y: [a,b] — IR? near its end points, we know its total 
curvature f H k ds up to an integer multiple of 27x. This follows from the results in 
Chap. 3. Here we prove a similar result for surfaces f: M — R? in three-space: If 
we know f near the boundary of M, we know its total Gaussian curvature f m K det 
up to an integer multiple of 27. However, unlike the situation for plane curves, the 
integer in question is already completely determined by the topology of M. 


10.1 Curves on Surfaces 


Definition 10.1 


Let f: M — R? be a surface with unit normal field N and y: [a,b] > M a 
curve in M. Then the pair (y, f) is called a curve on the surface f. The space 
curve 


y=foy 


is called the trace of (y, f). The velocity of (y, f) is defined as |y’| and 
accordingly the derivative with respect to arclength of a function g: [a, b] > R* 
is to be interpreted as 


dg g 


ds |p 


The unit tangent T of y is called the unit tangent of (y, f) and the unit normal 
field 


N=Noy 
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along y is called the surface normal of (y, f). The unit normal field 
B=TxN 
along y is called the binormal of (y, f). 


If (y, f) is a curve on the surface f, then (y, N ) defined as above will be a 
framed curve according to Definition 5.11. 


Definition 10.2 


If (y, f) is a curve on the surface f and T, Ñ, B are defined as in Definition 10.1, 
then 


(i) The normal curvature of (y, f) is defined as 
kn SN T). 
Gi) The geodesic curvature of (y, f) is defined as 
kg = (B',T). 
(Gii) The geodesic torsion of (y, f) is defined as 
t = (N’, B). 


Traditionally, curves on a surface f for which one of these quantities vanishes 
are designated by special names (see Fig. 10.1): 


Definition 10.3 
Let (y, f) be a curve on the surface f: M —> IR?. Then 


(i) (y, f) is called an asymptotic line if its normal curvature x, vanishes. 
(ii) (y, f) is called a geodesic if its geodesic curvature «, vanishes. 
Gii) (y, f) is called a curvature line if its geodesic torsion t vanishes. 


> Remark 10.4 The geodesic in Fig. 10.1 illustrates nicely that geodesics are 
locally length minimizing, but globally they are not necessarily the shortest path 
between two points. 
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Fig. 10.1 An asymptotic line (top), a geodesic (middle) and a curvature line (bottom) on a torus 


10.2 Theorem of Gauss and Bonnet 


Let M C R? be a compact domain with smooth boundary. By Definition 6.1 and the 
arguments surrounding Fig. 8.5, each of the n components of the boundary 0M can 
be parametrized by a closed curve y;: [0, 27] —> M. Given a surface f: M —> R, 
we define the total geodesic curvature of the boundary 0M by summing up the 
integrals of the geodesic curvature Kg over the corresponding curves (yj, f) on the 
surface f (Definitions 10.1 and 10.2): 


n 


Kg i= Ko ds. 
[ead] « 


j=l 
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Definition 10.5 


Let M cC R? be a domain with smooth boundary having k components and n 
boundary curves. Then 


xX(M) := 2k —n 


is called the Euler characteristic of M. 


Theorem 10.6 (Gauss-Bonnet Theorem) 

Let M C R? be a compact domain with smooth boundary having k connected 
components. Assume that the boundary 0M has n components. Let f: M > 
R? be a surface, K: M — R its Gaussian curvature and det its area form. 
Then 


f Kas f Kg ds = 2x (M). 
M aM 


Before we proof the theorem we note that we may always choose a vector field 
Z €e IT (TM) with (Z, Z} = 1, for example one could take Z = mU . Moreover, 
we will make use of the following helpful observations. 


Lemma 10.7 

Let M C R? be a compact domain with smooth boundary, y : [—x, n] > R? 
be a parametrization of a boundary curve of M and Z € T (T M) be a vector 
field with (Z, Z) = 1. Then, there is a smooth function a: R —> R. with 
a(x + 27) = a(x) + 2x £ for some £ € Z such that 


T =cos(a)Z oy +sin(a)JZ o y, 
where T = sea is the unit tangent vector of gamma. 
v¥ 
Proof. For each x € [—z, 1] we can choose an a(x) € R such that 


T(x) = cos(a(x))Z o y (x) + sin(a(x))JZ o y (x). 
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Due to the smoothness of the boundary curve, locally this choice of œ can be made 

in a smooth fashion. Therefore, w := a’ is well defined, so that we can safely define 
x 


a(x) = a0) + f w. 
0 


In particular, 
T 
f w =a(mT)— a(r) = 27k 
=T 
for some £ € Z, which is exactly the tangent winding number of the curve 


Wn Zoy) \ op 2 
Co a i el 


which we will denote by €(Z, (-, -)). oO 


Lemma 10.8 
L(Z, (-, -)) is independent of the chosen metric. 


Proof. Let (-,-)~ be any other metric. Then for any £ € [0,1] also (-,-); := 
(1 —t)(-,-) +t, -)~ is again a Riemannian metric. As £ is an integer and depends 
continuously on £ we conclude that it is constant and therefore @(Z, (-,-)) = 
L(Z, (aey J: oO 


Theorem 10.9 
Let M C R? be a compact domain with smooth boundary having k 


connected components and n boundary curves. Let y,..., Yn be closed 
parametrizations of the n boundary curves of M and £, . . . , £n be the tangent 
winding numbers of V1, . . . , Yn with respect to a unit vector field Z € T (T M). 
Then 


XZ) = 27 x (M). 


= 


Proof. Without loss of generality let (-,-) = (-, +)}g2 . For j = 1,..., denote the 
tangent winding number of y; by €;(Z). Then either ¢;(Z) = 1 Gf yj parametrizes 
the outer boundary of one of the components of M) or ¢;(Z) = —1. Since there are 
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k components and n — k interior components, we have 
n 
Pes =k—(n—k)=2k—n=x(M). 
j=l 
| 


Proof of Theorem 10.6 Choose a vector field Z € F(T M) with (Z, Z) = 1 and 
define a 1-form ņ € Q'(M) by 


n(X) = (VxZ, JZ). 
Think of ņn(X) as the rotation speed of Z in the direction of X. Because of 
(VxZ, Z) = 0 (which follows from differentiating (Z, Z} = 1) and and since 
Z, JZ is a positively oriented basis of T, M we must have 


VxZ = n(X)JZ. 


Using this, (ii) of Theorem 9.2, the Gauss equation (Theorem 9.8) and (i) of 
Theorem 9.4 we find 


dn(U, V) = dun(V) — dyn(U) 
= (Vu VvZ, JZ) + (Vv Z, JVuZ) 
— (Vy VuZ, JZ) — (VuZ, JVyZ) 
= (Vy VyZ, JZ) — (VvyVuZ, JZ) 
= —K det(U, V) 


and therefore 
dn = —K det. 


In particular, this means that dn does not depend on our choice of Z. Therefore 


Stokes’ theorem implies 
— f K det = f n. 
M aM 


For each of the boundary components, the geodesic curvature of y; := f o yj can 
be expressed as Kj ‘= (Ti, Nx T, Then, with 


(df(Zoy))' =—(Ay', Zoy)N+df (VyZoy) 
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we have 

kj = (Ti, N xT) 
(df (T;))', N x df (T)) 
=a'+n(y’). 


II 


Putting everything together we obtain 


= aD HZ | n 


j=l am 
= 27 x (M) -f K det. 
M 
0 
It is quite striking that the total amount of Gaussian curvature (in the sense of 


f m K det) is completely determined by the geometry of f near the boundary of 
M (see Fig. 10.2). 


AA 


Fig. 10.2 Even if we do not know the shape of a rounded cone near its tip (only revealed under a 


microscope), the integral of the Gaussian curvature can be deduced from the opening angle of the 
cone 
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Example 10.10 

Suppose M is a disk which is mapped to the top half of a round sphere with radius r > 0 by 
f. Then the boundary curve lies on the equator which is known to be a geodesic, i.e. kg = 0. 
Therefore, the Gauss-Bonnet theorem yields 


/ K det = 5, A(f(M)) = 27. 
M 


Example 10.11 
If M is an annulus and f maps it onto a cylinder, then x (M) = 0 and K = 0, so the Gauß-Bonnet 


formula yields 
f Kg ds = 0. 
aM 


10.3 Parallel Transport on Surfaces 


In Sect.5.1 we studied the normal transport P : T(a)+ —> T(b)+ of a curve 
y: [a,b] > R? with unit tangent T. A closer look reveals that in order to define 
P only the smooth map T: [a,b] —> S? is needed. Therefore, given a surface 
f: M — R? and a smooth map y: [a,b] — M, we can we can use the same 
strategy in order to transport tangent vectors W € Ty, q)M to tangent vectors 
PW) € TywM: 

Definition 10.12 

Let f: M > R? be a surface with unit normal field N and y: [a,b] > Ma 

smooth map. Define Ñ : [a, b] > S? by 

N:=No y 


and for W € T,q)M define the parallel transport map P, (W) € Tyœ)M in 
such a way that 


df (P(W)) := Z(b) 
where Z: [a, b] + R? solves the initial value problem 


Z(a) =df(W) 
Z! =—(Z,N’)N. 
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N plays exactly the same role here as T did in Sect.5.1. Hence, for the same 
reasons as in Sect. 5.1, we have (Z, N} = 0 and indeed for all x € [a, b] the vector 
Z(x) is an element of df (Ty (x) M). Furthermore, 


Py: Tya) M -=> Tya; M 


is an orientation-preserving orthogonal map with respect to the metrics induced by 
f on Tya M and Tyo) M. 

The derivative Z'(x) is a multiple of N(y(x)), so it has no component in 
df (T,(x)M). In the spirit of the Sect. 9.1 (Levi-Civita connection), where a deriva- 
tive VxY of a vector field Y was defined in terms of the tangential component of 
dx (df (Y)), this means that P can be viewed as parallel transport along y. 

Imagine a pendulum swinging at a point of a surface f: M — R3 subject to 
gravity pointing away form the unit normal of the surface. Suppose we transport 
the swinging pendulum along a path f o y where y: [a, b] —> M is a smooth map 
and that the plane in which the pendulum swings initially is given as df (W) where 
W € T, a)M is a unit vector with respect to the induced metric. Then Physics tells 
us that the plane in which the pendulum swings once it arrives at f(y(b)) will be 
given by the unit vector df (P(W)). 

In the special case where f parametrizes the surface of the earth and the move- 
ment y corresponds to the rotation of the earth, this effect can be experimentally 
verified and is known under the name of Foucault’s pendulum (see Fig. 10.3). 

As in Sect.5.1, if we choose unit vectors (with respect to the induced metric) 
Wa € Tja M and Wa € Tœ M, we can measure the parallel transport along y 


Fig. 10.3 In 1851, Léon Foucault build a Foucault pendulum to demonstrate the rotation of the 
earth (left). The phenomenon can be understood with the concept of parallel transport, where a 
tangent vector is transported along a circle of latitude (right) 
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Fig. 10.4 Parallel transport of a tangent vector along a closed curve on a surface with positive 
Gaussian curvature (/eft) and along the boundary of a surface with negative Gaussian curvature 


(right) 


by an angle Pw € R/27z. For closed curves y this angle does not depend on the 
choice of Wa and Wp as long as we make sure that Wa = Wp. In the special case 
where y parametrizes the boundary 0M of M, this angle can be expressed in terms 
of the total Gaussian curvature of f (see Fig. 10.4). 


Theorem 10.13 

Suppose that f : M — R? is a surface and that M has only a single boundary 
component parametrized by a curve y: [a,b] — M. Then the monodromy 
angle of y satisfies 


M(y) =f K det mod 2r Z 
M 


where K is the Gaussian curvature of f. 


Proof. Let us assume that y has unit speed with respect to the induced metric and 
therefore T := 7’ is the unit tangent field of 7 := f o y. Define Ñ := N o y where 
N is the unit normal field of f. Let W and Z be defined as in Definition 10.12. Then 
there is a smooth function a: [a, b] —> R such that 


Z = cosa T + sing Ñ x T. 
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We denote by kg = (T', Ñ x T) the binormal curvature of the framed curve Q, Ñ). 
Because Z’ is normal, we have 


0 = (Z', N x Z) =a' + kg. 


Finally, by the Gauss-Bonnet Theorem 10.6 we have 


b b 
MP) = a(b)- aa) = f a=- f x= f K det mod 2x2. 
a a M 
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Closed Surfaces 1 1 


We define a closed surface as a surface f: M — R? whose boundary components 
have been matched in pairs in such a way that f as well as its unit normal N are 
continuous across the boundary. This allows us to prove an analog of the fact that the 
tangent winding number of a closed plane curve is an integer: The total Gaussian 
curvature f m K det of a closed surface f: M > IR? is equal to 27 x (M) where 
x (M) is the Euler characteristic. 


11.1 History of Closed Surfaces 


Our goal here is to define “closed surfaces” in such a way that we are able to prove 
an analog of Theorem 3.8, which says that the turning number of a plane curve is an 
integer. Furthermore, in Sect. 13.1 we want to discuss for closed surfaces the analog 
of the total squared curvature of a curve. 

Our approach will be based on the very idea that was already at the heart of the 
1845 paper by Mobius where closed surfaces were studied for the first time: By 
cutting them into horizontal slices, Möbius decomposed closed surfaces into pieces 
each of which can be parametrized by a compact domain with smooth boundary in 
R?. Figure 11.1 is adapted from the paper by Möbius. This very idea was already 
the motivation for us to allow for disconnected domains in the case of surfaces and 
will be formalized in Sect. 11.2. 

More details on the early history of surface theory can be found in an article by 
Peter Dombrowski [11]. 

A more advanced way to define closed surfaces in IR” (that would not need to 
cut the surface into pieces that can be parametrized by planar domains) would be 
to define them in terms of smooth maps f: M — R” defined on 2-dimensional 
compact manifolds M. Such manifolds were first defined in 1910 by Hermann Wey] 
in a famous book with the title “Die Idee der Riemannschen Fläche” [44]. 
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Fig. 11.1 Möbius decomposed closed surfaces into pieces that can be parametrized by compact 
domains in R? with smooth boundary (modeled after Mobius’ original sketch in [29]) 


On the other hand, the fully developed version of the Gauss-Bonnet theorem 
(which we will prove in the next chapter) is already contained in the 1903 thesis of 
Werner Boy [8], that he did under the supervision of David Hilbert. 

Modern treatments of Differential Topology (like the books by Andrew Wallace 
[42] and Morris Hirsch [16]) often discuss surface topology in their last chapters. 
The main work there goes into proving (with the help of Morse theory) that indeed 
every compact 2-dimensional manifold can be decomposed into pieces each of 
which can be parametrized by a compact domain with smooth boundary in R?. 
Therefore, the work that will be done in the next two chapters would not become 
obsolete even if we had manifolds at our disposal. 


11.2 Defining Closed Surfaces 


Suppose that for a surface f: M —> R? the boundary components of M match up 
in pairs in such a way that, given suitable parametrizations of the boundary curves, 
corresponding points of 3M are mapped to the same points in IR. If in addition also 
the unit normals of f fit together up to sign on 0M, we consider f (together with a 
specification of the boundary matching) as a closed surface (Fig. 11.2): 
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Fig. 11.2 A closed surface f 


Definition 11.1 


Let M C R? be a compact domain with smooth boundary and f: M > R? a 
surface with unit normal N. We parametrize the boundary curves of M by closed 
curves 
VYis-++9 Yn: [-x, 7] => R? 
and define curves yj,..., Yn: [-7, m] > R3 by 
y j= f ° Y, j: 


As in Definition 10.1, we equip the closed space curves y; with unit normal fields 


Nj := N o yj. Let 
p: {1,... n} > {1,... 7} 
a bijective map such that 
(pop) =j 


for all j. Then the pair (f, p) is called a closed surface if there are signs 
€1,---,€, E {—1, 1} such that for all j € {1,..., n} we have: 


G) If p(j) # j then 


Yog) Œ) = Fj (ejx) 


Noo Œ) = —e; Ñ; (ejx). 
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Fig. 11.3 The surface in Fig. 11.2 made into a non-closed surface by applying a small translation 
to each piece 


Gi) If o(j) = j then €j = 1 and 


š Yx +r) forx e[-x,0) 
ya) = į 


yj(x — x) forx e [0,7] 


ss —Ni(x+z) forx € [—z, 0) 
Nœ S= 
—N;(x— x) forx € [0,7]. 
It is easy to see that such €1, . . . , €n are uniquely determined by f and p. We say 
that a closed surface is oriented if €e; = —1 for all j € {1,..., n}. 


Figure 11.3 shows the shape of the individual pieces that are being glued in 
Fig. 11.2. It has k = 6 components and n = 18 boundary curves. 

Here is another example: M now consists of a disk with boundary yı and an 
annulus with boundary curves yz and y3. First, we tentatively define f on the disk 
bounded by yı and obtain the cap on the upper right of Fig. 11.4. Postponing for the 
moment the task (indicated by the double-arrow on the right) of gluing yı to y2, we 
first glue y3 to itself and obtain a Möbius band (on the bottom of the lower right of 
Fig. 11.4): 

By growing the Möbius band (see Fig. 11.5) we finally obtain the closed surface 
we wanted to construct: 

This surface (fully closed in Fig. 11.6) was found by Werner Boy in 1903 and is 
called the Boy surface. 

Figure 11.7 shows two surfaces which are obtained by gluing the boundary curve 
of an annulus to itself appropriately. Even though both compact domains have k = 1 
components and n = 2 boundary loops, the distinct maps f, f lead to distinct closed 
surfaces. In particular, although the map p is the same, they have opposite sign €. 
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k=2 
n=3 
p(3) =3 
ep=-l 
€&3 = 1 


Fig. 11.4 The annulus part of the domain on the left has its red boundary component glued to 
itself. After growing the resulting Möbius strip, the other boundary component can be glued to the 
image of the disk part of the domain. The result is the so-called Boy surface 


11.3 Boy’s Theorem 


Definition 11.2 


We say that a surface f: M — R? closes up if there is p such that (f, p) is a 
closed surface in the sense of Definition 11.1. 


Recall that for every closed plane curve y: [a,b] — IR? there was an integer 


n € Z such that 
b 
f k ds = 27n. 
a 


Surprisingly, the analog of this fact in the context of surfaces (cf. Theorem 11.3) 
does not involve any information about the specific way in which f closes up, but 
only depends on properties of the domain M. The theorem is a variant of the Gauss- 
Bonnet Theorem 10.6. Usually, it would be called by the same name. However, 
historically this is not quite correct. This theorem was in fact the main result of the 
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Fig. 11.5 A growing Möbius strip can be capped off to form a Boy surface 


thesis of Werner Boy [8], written in 1903 under the supervision of David Hilbert. 
For this reason, we name it after Boy: 


Theorem 11.3 (Boy's Theorem) 
Let f: M — R? be a surface that closes up. Then the Gaussian curvature K 


of f satisfies 


/ K det = 2x x(M). 
M 
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Fig. 11.6 The Boy surface is a closed, non-oriented surface 


Fig. 11.7 After pushing the two boundary curves together, we obtain a closed surface which is 
oriented—a torus (left), or a closed surface that is not oriented—a so-called Klein bottle (right) 


Before we give the proof, we introduce the notion of an orientation cover of a 
closed surface. Given a closed surface (f, 0) with f: M > R3, we can define an 
oriented closed surface (f, 6) in the following way: 

Let us use M_, as another name for M and, using an orientation-reversing 
isometry g: R? — RŽ, we place a second copy Mı = g(M) into R? in such a 
way that M_; and M; are disjoint. Then we define 


M = Mı UM, 
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and 
ee A fp) if peM_, 
f: M > R?, f(p)= T . 
(fog )(p) if peMi. 
We can label the boundary curves of M by the elements of {—1, 1} x {1,..., n} 


and parametrize them by maps 


5 yj if i=-—l1 
vap: R> OM, Yip = 47 oe 
Xr goyj(—-x) if i=l. 


Finally, we define 
p: {-l, 1} x {1,..., 2} = {-1, 1} x {1, iath pi, j) = (-€; i, p(j)). 


We now leave it to the reader to check that ( Ô, P) is an oriented closed surface, i.e. 
we obtain a closed surface by setting čq, j = —1 forall (i, j) € {-1, 1}x{],..., n}. 


Definition 11.4 


The closed surface (f, 6) constructed above is called an orientation cover of f. 


Proof of Theorem 11.3—Boy’s Theorem If p has no fixed points (no boundary 
component is glued to itself), one just has to note that the existence of p (making 
(f, e) into a closed surface) implies that in Theorem 10.6 the total geodesic 
curvatures of the individual boundary curves cancel in pairs. If po has fixed points, 
we note that the ¢ of the orientation cover has no fixed points and therefore our 
theorem holds for f . Dividing both sides of the resulting equation by two, we see 
that our theorem also holds for f. o 


11.4 The Genus of a Closed Surface 


The Euler characteristic of a closed surface was solely a property of its domain M, 
the specific way the various boundary curves are glued is irrelevant for the Euler 
characteristic. There is another number associated with a closed surface (f, o), the 
so-called genus, that depends on the gluing correspondence p: 

Suppose M cC R? is a domain with k components and n boundary curves. 
Consider the map that assigns to each j € {1, ..., n} the index c(j) € {1,..., k} of 
the component of M to which the jth boundary component belongs. Let us consider 
the graph G whose vertex set is {1, ..., k} and in which two vertices £, Ê with £ # i 
are connected by an edge if and only if there is an index j € {1,...,n} for which 
c(j) = £ and c(p(J)) = £, which means that the components of M with indices j 
and j are glued via one (or more) of their respective boundary curves. We say that 
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two vertices £ and £ of G are connectable in G if it is possible to travel from £ to é 
by following edges. Connectability is an equivalence relation and the corresponding 
equivalence classes are called the connected components of G. 

Definition 11.5 

If {€1,..., zh is a component of the graph G, then 

f= F Me, U...UMe, 
closes up with boundary gluing p read off from (f, o). We call the resulting 
closed surface (f, 0) a component of (f, o). We call (f, p) connected if it has 


only one component. 


So the components of a closed surface are in one-to-one correspondence with the 
components of its associated graph G. 


Definition 11.6 


Let M be a compact domain with k components and n boundary curves. Let 
(f, e) be a closed surface with f: M — R?. If (f, e) has m connected 
components, we define the genus of ( f, p) as 


n y 5 
=o m. 
a 
In terms of the genus, the Gauss-Bonnet formula takes the form 
f K det = 47r (m — g8). 
M 
The first surface featured in Sect. 11.2 has genus g = 4, the Klein bottle has 


genus g = 1 und the Boy surface has genus g = 5. The two surfaces in Fig. 11.8 
have genus g = 3 and genus g = 2 respectively. 
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Fig. 11.8 Non-oriented surfaces of genus g = 3 (left) and g = 2 (right). They are obtained by 
smoothly gluing handles onto a Boy surface or respectively a Klein bottle 
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Variations of Surfaces 1 2 


We derive the basics of Vector Calculus on surfaces and explore variations of 
surfaces. In particular, we compute the variational derivative of the area form det 
and of the shape operator A. We show that the critical points of the area functional 
are the surfaces with mean curvature H = 0. If we constrain the enclosed volume, 
the critical points of area are the surface with constant mean curvature. These results 
mirror the situation for plane curves, where the analogous variational problems lead 
to straight lines (x = 0) or circles («x = const). 


12.1 Vector Calculus on Surfaces 


Throughout this section, M C IR? is a Riemannian domain, f: M —> R? a surface 
and (,) its induced metric. We will only use the area form det, the 90°-rotation 
J and the Levi-Civita-connection V, which by Theorems 6.22, 6.23 and 9.10 are 
already determined by the induced metric. This means that this section is dealing 
only with intrinsic geometry. 

If g € C™(M) is a smooth function, then for each p € M the restriction 
(d Bir; m is a linear map on TpM and the restriction (, Iz, MxTpM is a Euclidean 
scalar product. Therefore, there is a unique vector Y (p) € T,M such that dg (X) = 
(Y (p), X) for all X € T M. The smoothness of the vector field Y defined in this 
way follows in the usual way, see for example the proof of Theorem 6.20. This leads 
us to the following: 


Definition 12.1 


For g € C (M) there is a unique vector field 


gradg € T (TM) 
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grad(zo f) 


Fig. 12.1 The gradient vector field of the function z o f (z being the third coordinate function on 
R3) for a surface f: M — R?. On the left, the value of z o f is indicated by color-coding 


characterized by the fact that for all vector fields X € r (TM) we have 
dg(X) = (grad g, X). 
The vector field grad g € y (T M) is called the gradient of g. 
So a function g € C™(M) gives us a vector field gradg € T'(TM) (see 


Fig. 12.1). On the other hand, by taking the trace of the endomorphism field VY, 
a vector field Y € '(7M) gives us a function div Y € C% (M). 


Definition 12.2 
For a vector field Y € Ir (T M) the function 


div Y: M—>R, divY = tr(VY) 
is called the divergence of Y. 
The following theorem from Linear Algebra is useful for calculating the trace of 


an endomorphism field. 


Theorem 12.3 
Let W be a 2-dimensional vector space with a determinant function det and 
A: W — W a linear map. Then for any two vectors X, Y € W we have 


det(AX, Y) — det(AY, X) = tr A det(X, Y). 
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Proof. If X and Y are linearly dependent, both sides of the equation vanish. 
Otherwise, X and Y form a basis of W and we can write 


AX =aX +cY 
AY = bX + dY. 
Our claim now follows from 
trA=a-+d. 


oO 
For the divergence of the product of a function and a vector field we have a Leibniz 
formula: 
Theorem 12.4 


For g € C®(M) and Z € T (TM) we have 


div(gZ) = (grad g, Z) + gdiv(Z). 


Proof. With the notation G := grad g and with the help of Theorems 6.24 and 12.3, 
for X, Y € T (T M) we have 
div(gZ) det(X, Y) = det(Vx(gZ), Y) — det(Vy (8 Z), X) 
= det((X, G) Z + Vx Z, Y) — det((Y, G)Z + gVyZ, X) 
= — det((Y, G) X — (X, G)Y, Z) + g div(Z) det(X, Y) 
= (det(—J G, Z) + g div(Z)) det(X, Y) 
= ((G, Z) + g div(Z)) det(X, Y). 


Definition 12.5 


The divergence of the gradient of a function g € C% (M) 
Ag := div grad g 


is called the Laplacian of g. 
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The divergence of a 90°-rotated gradient vanishes: 
Theorem 12.6 
For every g € C% (M) we have 


div(J grad g) = 0. 


Proof. Using again the notation G := grad g, by Theorems 9.6 and 12.3 we obtain 


div(J grad g) det(U, V) = det(Vy (JG), V) — det(Vy (JG), U) 
= (VyG, U) — (VuG, V} 
= dy (G, U) — dy (G, V) 
= dydyg — dydy g 
=0. 


oO 


The theorem below is a reformulation of Stokes Theorem in terms of vector fields 
instead of 1-forms. The integral 
f gds 
aM 


of a function g: 0M — R is defined in the same way as for total geodesic 
curvature—as the sum of integrals over the boundary loops. 


Theorem 12.7 (Divergence Theorem) 
Let Y € T (TM) be a vector field and B the outward-pointing unit normal 
field on the boundary 0M. Then 


f axy act= [ (Y, B) ds. 
M aM 


Proof. Define a 1-form w € Q!(M) by setting for X € TpM 


w(X) = (JY (p), X). 
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Then, by Theorems 9.2, 9.4, 9.6 and Lemma 12.3, 


dw(U, V) = dyw(V) — dya(U) 
= (JVuY, V) + (JY, Vu V) — (JVvY, U) — (JY, Vy U) 
= det(Vy Y, V) — det(Vy Y, U) 
= tr(VY) det(U, V). 


Therefore dw = div Y det. Using again the notation of the proof of Theorem 10.6 
and applying Stokes Theorem 7.15 we obtain 


f awy det = [ao 
M 

I 7 
ba 


oe 


(JY, T)d 


12.2 One-Parameter Families of Surfaces 


Throughout this chapter M C RR? will be a compact domain with smooth boundary 
and [fo, t1] C R a closed interval. 


Definition 12.8 


Let gr: M —> R” a smooth map, defined for each t € [to, t1]. Then the one- 
parameter family of maps [tọ, t1] Ə f + g; is called smooth if the map 


M x [to, ti] > R”, (p, 1) > g:(p) 
is smooth (as always, in the sense of Remark 1.2). 
> Remark 12.9 The variable ¢ is also referred to as the time. 
Given a smooth one-parameter family 


the (g: M —> R”), te [t,t] 
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of maps and a vector field X € Ir (T M), also 
t œ> dyg 


is a smooth one-parameter family of maps dxg;: M — R”. The same holds for 
t œ> & where g;: M — R” is defined as 


: 4 
& (p) := TE &:(p). 


The following fact will be used many times in upcoming chapters: 


Theorem 12.10 

For a smooth one-parameter family of maps t +> g; from M to R”, the 
directional derivative in the direction of a vector field X € T (T M) commutes 
with the time derivative: 


(dy gry = dx &. 


Proof. In the special case where X is one of the coordinate vector fields U and 
V, this is just the fact that partial derivatives of the smooth map (p,t) => g;(p) 
commute. In the general case, we can write 


X=aU+bv 
where a, b € C®(M) are independent of t. Then 


(dx gr)’ = (a dugi + b dyg) 
=ady& +bdy& 
= dx &. 


Definition 12.11 


A smooth one-parameter family ¢ > g, of maps from M to R” is called a 
variation of a smooth map g: M — R” if 


tp <O0<t 
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and 


In this context, we will also use the notation 
& = ĝo- 
One should compare the arguments below with our reasoning in Sect. 2.4. 


Definition 12.12 


A variation of a surface f: M — R” is a smooth one-parameter family of 
surfaces 


fi: M > R”, tel[-e6e] 


such that 


The map f: M —> R” defined as 
Ê = fo 
is called the variational vector field of the variation t > fi. 


Definition 12.13 
Let M c RÊ? be a compact domain with smooth boundary. Suppose we have a 
way to assign to each surface f: M — R” areal number €(f). Then € is called 
a smooth functional if for every smooth one-parameter family 

tre fe t€ [t,t] 
of surfaces f: M — R” the function 


[to, tı] = R, the E( ft) 


is smooth. 
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In many circumstances, we want to consider only variations of f: M — R” that 
keep the surface fixed near the boundary 0M: 


Definition 12.14 


Let M C R? be a compact domain with smooth boundary and f: M > R” a 
surface. Then a variation 


tre fi, te[-e,€] 


of f is said to have support in the interior of M if there is a compact set Mo C 
M such that for all p € M, p ¢ Mo we have 


fi(p) = f(p) forall t¢[-e,€]. 


Definition 12.15 


Let M C RŽ be a compact domain with smooth boundary and € be a smooth 
functional defined on the space of surfaces f: M — R”. Then a surface 
f: M — R” is called a critical point of € if for all variations t œ> f, of f 
with support in the interior of M we have 


d 
oF Lee =0. 


Definition 12.15 spells out the notion of an equilibrium of a variational energy €, 
to which we will refer to in later sections. Moreover, one should note that, as already 
explained in the beginning of Sect. 2.4, we will work with a definition of a critical 
point under constraints that is slightly stronger than the standard one. 


Definition 12.16 


Let M C R? be a compact domain with smooth boundary, f : M —> R” a surface 
and €,€ two smooth functionals on the space of all surfaces f: M > R”. 
Then f is called a critical point of € under the constraint of fixed € if for all 
variations t +> f; of f with support in the interior of M 


2) 326 

dtl o 
implies 

d 

=> E=0. 

dt | 9 
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Using the Linear Algebra Therorem 2.21 in the same way as we used it in 
Sect. 2.4, we obtain 


Theorem 12.17 

Let M C R? be a compact domain with smooth boundary and €E, È two 
smooth functionals on the space of all surfaces f: M — R”. Suppose we 
have a way to associate to each surface f: M — R” smooth maps 


Gr, Gr: M => R” 


such that for all variations t œ> f, of f with support in the interior of M we 
have 


d : 

Be E 
dt 1=0 Pa f 

d 5 EA 

= é= f Gee. 
dt m Me f 


Then f is a critical point of E under the constraint of fixed E if and only if 
there is a constant à € R such that 


Gf = AG p. 


For reasons already explained in Sect. 2.4, we call à a Lagrange multiplier for 
the constraint of fixed E. 


12.3 Variation of Curvature 


Given a smooth variation t +» f; of a surface f, we are mainly interested in the time 
derivative at time zero of quantities like the area form det, or the shape operator A; 
associated with the surfaces f;. In situations where it clear with which variation 
t œ> fi we are dealing, we will usually drop the index zero when we mean the time 
derivative at time zero. So, for example, we will write 


A= Ap. 
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Theorem 12.18 

Let f: M — R? be a surface with unit normal N, shape operator A and 
Levi-Civita connection V. Let t +> f; be a variation of f whose variational 
vector field 


f=$N +adf(Z) 
is described in terms of a function o € C°(M) and a vector field Z € 
I'(TM). Denote by N, and A, the unit normals and the shape operators of 
the surfaces fi. Define vector fields G, W € T (TM) as 


G := grad ġ 
W := AZ — G. 


Then 
df (X) = —(W, X)N + df (¢AX + VxZ) 
N = df(W) 
det = (26H + div Z) det 
Å = VzA — VG — $A? + A(VZ)—(VZ)A 


‘ 1 
Å = dzH — z4vG = OH? R 


Proof. The proof of the first equation is straightforward: 


df (X) = do(X)N + ¢ df (AX) — (AX, Z)N + df (VxZ) 
= (G — AZ, X)N + df (AX + VxZ) 


Differentiating (N, df(X)) = 0 with respect to time we obtain 
(N, df (X)) = —(N, d f (X)) = (W, X) = (df (W), df (X). 


This holds for all X € TM and this implies the second equation. For X, Y € T,M 
we know that Ñ (which is orthogonal to N), df(X) and df(Y) are linearly 
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dependent. Using this and Theorem 12.3 we obtain 


det(X, Y) = det(N, df (X), df(Y))* 
= det(N, df (X), df (YY) + det(N, df (PAX + Vx Z), df(Y)) 
+ det(N, df (X), df (QAY + VyZ)) 
= det(@AX + Vx Z, Y) + det(X, PAY + VyZ) 
= tr(pA + VZ) det(X, Y) 
= (26H + div Z) det(X, Y). 


This proves the third equation. For the fourth equation, consider the directional 
derivative of the second equation in the direction of X and make use of the first: 


— (AX, W)N +.df(VxW) = dN(X) 
= (dN(X)) 
= (df (AX)Y 
= df (AX) +df (AX) 
= —(W, AX)N +df(@A?X + VayZ) + df (AX) 


The normal part of this equation is satisfied automatically. The tangential part, 
together with the Codazzi equation (Theorem 9.7) gives us 


AX =VxW — ¢A°X — VaxZ 
= Vx(AZ) — VxG — p4? X — VayZ 


= (VyA)Z + AVxZ — VxG — $A°X — VaxZ 


= (VzA)X + A(VZ)(X) — (VG)(X) — A2X — (VZ)(AX). 


This proves the fourth equation. For the fifth we take the trace of the fourth and 
multiply by 5. The last two terms in the fourth equation do not contribute because 
we see here the commutator of two endomorphisms A and V Z, which always has 
zero trace. Regarding the first term, one can verify (for example by taking the 
directional derivative of the equation in Theorem 12.3 in the direction of Z) that 
indeed for any endomorphism field A 


tr(VzA) = dz(tr A). 
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Finally, by diagonalizing A one can easily check the equality 


1 
JTA? =2H*— K. 


12.4 Variation of Area 


Variations of surfaces (as defined in Definition 12.12) are needed in order to define 
and determine those surfaces that represent equilibria of geometrically interesting 
variational functionals. 

Examples of smooth functionals of surfaces are the Willmore functional (to be 
introduced in Sect. 13.1) and the cone volume that will be defined in Sect. 12.5. In 
this chapter we will focus on the area functional 


A= | det ;. 
M 


Theorem 12.19 (First Variation Formula of Area) 
As in Theorem 12.18, suppose the variational vector field of a variation t > 
fi: of a surface f: M — R? is written as 


f =6N +df(Z) 
with p € C (M) and Z € F(T M). Then 


d 


dt 


Ah) =2f oH dert | (Z, B) ds 
(=i M ƏM 


where B is the outward pointing unit normal on ðM. 


Proof. By Theorem 12.18 and the Divergence Theorem 12.7, 


d 


dt 


A= f det 
M 


=2 f oH det+ f div Z det 
M M 


=2 f oH det+ f (Z, B) ds. 
M ƏM 


t=0 
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Fig. 12.2 The Schwarz-P 
minimal surface 


Definition 12.20 


A surface f: M — R? is called a minimal surface if it is a critical point of the 
area functional A. 


Figure 12.2 shows a minimal surface whose six boundary curves are all mapped onto 
prescribed circles. In fact, it is here a solution of the so-called Plateau problem, 
which means that it minimizes area among all surfaces whose boundary is mapped 
onto a prescribed set of curves. 


> Remark 12.21 The Plateau problem was first solved by Jesse Douglas [12] and 
Tibor Rado [32] independently. 


Theorem 12.22 
A surface f: M —> R? is a minimal surface if and only if is mean curvature 
H vanishes. 


Proof. If H = 0 and t > f is a variation of f with support in the interior of M, 
then Z vanishes near the boundary of M and by Theorem 12.19 the variation of area 
is zero. Conversely, suppose that f is a minimal surface but there is a point p € M 
for which H (p) # 0. Then there is such a p also in the interior of M, so we assume 
pE M. Let us treat the case H (p) > 0, the case H(p) < 0 being similar. Then 
we can construct a bump function g € C™(M) such that g vanishes outside of a 
compact set contained in the interior of M and 


g(p) = 1 
H(q) <0 = g(q)=0. 
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Fig. 12.3 An Enneper 
surface is a minimal surface 


Then, for small enough e€ > 0, 


tre fi, tel-e,€] 
fr=ft+t-g-N 


(N being the unit normal of f) will be a smooth variation of f with support in the 
interior of M and 


d 
— am = f gH > 0, 
dt 1-0 M 
which contradicts our assumption that f is minimal. o 


As the reader may verify, the Enneper surfaces defined in Sect. 6.5 have mean 
curvature H = 0, so by Theorem 12.22 they are minimal surfaces. Figure 12.3 
shows one of these Enneper surfaces. 


12.5 Variation of Volume 
Definition 12.23 


Let M C R? be a compact domain with smooth boundary and f: M > R? a 
surface. Then the cone volume of f is defined as 


1 
VA = 5 / SRC Ys fd. 
M 


12.5 Variation of Volume 175 


V(f) can be interpreted as the volume covered by the map 
F: [0,1] x M > R?, F(s, p)=s- f(p). 
Here the “volume covered” should not be understood as the volume of the image 
F([0, 1] x M), but rather in the spirit of Theorem 8.17. At first sight, the cone 
volume does not look like an honorable geometric functional. For example, the 
version f = f +a of f that has been translated by a vector a € IR? in general 


does not have the same cone volume as f. On the other hand, for closed surfaces 
the cone volume is invariant under translations: 


Theorem 12.24 
If (f, p) is an oriented closed surface (Definition 11.1) and a € R?, then 


VAC +a) = V(f). 


Proof. Define a 1-form w € Q! (M) by 


o(X) = É deta, f,df(X)). 


Then 
1 1 
dæ(U, V) = zta, Í, fo)u — det(a, Í, fu)v = 3 det(a, tus fo) 
and therefore, by Stokes Theorem 7.15, 
vsta -VA = | do= | eai 
M aM 


The last equality follows from the fact that (f, o) is oriented, and therefore the 
integrals of w over the various boundary curves of M cancel in pairs. o 
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Fig. 12.4 The cone volume 
of f (left) and of a variation 
fi of f (right) 


Moreover, by almost the same reasoning as in the above proof one can show: 


Theorem 12.25 
Let M C R? be a compact domain with smooth boundary, 


the fr, te [-e,€] 


a variation with support in the interior of M of a surface f: M — R? and 
a € R°. Then 


d d 
T ili V(fı +a) = T S VC fi). 


Theorem 12.25 implies that for the purposes of variational calculus the cone 
volume V behaves in the same way as a translationally invariant functional (see 
Fig. 12.4). 

We can view df as an Rĉ?-valued 1-form on M. Given smooth maps f, f : M > 
R? we then obtain a scalar valued 1-form 


= 5 def, $ df) € ËM). 


Theorem 12.26 (First Variation of Cone Volume) 
Let f: M — R? be a surface. Then for every variation t > f; of f we have 


d 


A 1 R 
aE l (FN) det-+ f ee 
dt |0 M am 3 
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Proof. We have 


dt 


1 : : . 
Vf) = 3 1 (det(f, fus fo) + det(f, fus fo) + det f, fur A)) 
M 


t=0 
1 : . ; 
= J (derc/, fas fo t det, f, fo)u — det(fu, f, fv) 

M 


— det(f, f: fou) + det( f, fu, Êh — det( fv, fu, f) 
—det(f, fav, f)) 


al, (f, det(U, von) + f dw 
M 


M 


= f| imat fo. 


It is easy to see that, on its own, the cone volume functional does not have any 
critical points. However, we can use it in the context of variational problems under 
a volume constraint. Here is our first application of Theorem 12.17: 


oO 


Theorem 12.27 

Let M C R? be a compact domain with smooth boundary. Then a surface 
f: M > R? is a critical point of the area A under the constraint of fixed 
cone volume V if and only if the mean curvature H of f is constant. 


Proof. By Theorems 12.19, 12.26, and 12.17, f is a critical point of area under 
fixed cone volume if and only if there is a constant A € R such that 


HN=MN. 
oO 
The surface in Fig. 12.5 minimizes area among all surfaces that are bounded by the 


same six circles as the first surface shown in Sect. 12.4 and have a certain prescribed 
volume: 
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Fig. 12.5 A surface with the same boundary as the surface in Fig. 12.2. It is a critical point of the 
area functional under the constraint of having a prescribed cone volume 


BaD 


Fig. 12.6 A Wente torus—a closed surface of genus g = 1 with constant mean curvature H 4 0 


> Remark 12.28 In 1984 Henry Wente found a counterexample to a conjecture 
by Heinz Hopf which stated that every closed surface in R? with constant mean 
curvature is round sphere [43]. In Fig. 12.6 it is shown how the Wente torus can be 
build from a fundamental piece. 

Nevertheless, the conjecture is true if one demands that the surface is embedded in 
R3, or has genus g = 0. These results are due to Alexandrov [1] and Hopf [17]. 
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Willmore Surfaces 1 3 


The analog for a surface f: M — R? of the bending energy f. ? k?ds is the 
Willmore functional W ( f) = f m” 2 det. There are several versions of the Willmore 
functional, all of which are equivalent for the purposes of Variational Calculus. One 
of these versions is unchanged if we transform the surface by inversion in a sphere. 
The analogs of elastic curves are called Willmore Surfaces. 


13.1 The Willmore Functional 


In the context of curves y: [a,b] —> R? we studied in detail the total squared 
curvature f «ds (notation from the end of Sect. 7.2). What is the analog of this 
energy in the context of surfaces? 

One might say that k = 0 characterizes straight lines, which minimize length 
among all curves with the same end points. So Si a.b] k? ds measures the deviation 
from being length-minimizing. The analog of length-minimizing curves are area- 
minimizing surfaces, i.e. minimal surfaces, surfaces with mean curvature H = 0. 
So a natural analog of Ji abi k? ds can be defined as follows: 


Definition 13.1 
If f: Ma R? is a surface, then 


Wf) := f H? det 
M 


is called the Willmore functional of f. 
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Surfaces f that are critical points of the Willmore functional are characterized 
by the property that they are “as minimal as possible”, given that they are held fixed 
near the boundary of M. 

Alternatively, one might say that k = 0 only happens for straight line segments, 
so for surfaces we want to measure the deviation of being planar. Parametrizations 
of pieces of the plane are characterized by the fact that both principal curvatures 
vanish, so we want to measure the deviation of both «xı and x2 (not just their average 
H) from being zero. This reasoning leads to a different analog for Ji PEN k? ds: 


Definition 13.2 
Iff: M> R? is a surface, then 


1 2 2 2_K 
E(f) CAG ~ 5) ae 


is called the bending energy of f. 


Surfaces f that are critical points of the bending energy are characterized by the 
property that they are “as planar as possible”, given that they are held fixed near the 
boundary. 

Finally, one might formulate a different wish and ask for surfaces that are “as 
round as possible” which means they are “as spherical as possible”. In view of the 
Umbillic Point Theorem 8.12 this motivates the following definition: 


Definition 13.3 
If f: M — R? is a surface, then 


WF) := if (a Wa)? det = f (H? — K) det 
4 Ju M 


is called the conformally invariant Willmore functional of f. 


Note that the integrands in all three of the above energies differ only by a 
term proportional to K det, so the Gauss-Bonnet Theorem 10.6 tells us that for 
the purposes of Variational Calculus, (cf. Definition A.4) all three energies are 
equivalent to a large extent: 


Theorem 13.4 N 
Let f, f: M —> R? be two surfaces such that f (p) = f (p) for p outside of 
some compact set contained in the interior of M. Then 


W(f) — Wf) = E(/) — Ef) = Wf) - Wp). 
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For surfaces that close up we see that the difference between the three functionals 
only depends on the genus: 


Theorem 13.5 
Let f: M — R? be a surface that closes up with genus g. Then 


E(f) = W(f) + 2x(g — 1) 
W(f) = Wf) +4m(8 - 1). 


Theorem 13.6 
The estimates below are sharp, i.e. in each case there is a surface that closes 
up with the prescribed genus and which realizes the lower bound: 


(i) If M is connected and a surface f: M — R? closes up with genus 0, 
then 


W(f) = 4r. 


(ii) If M is connected and a surface f: M —> R? closes up with genus + 
then 


W(f) > 12r. 


(iii) If M is connected and a surface f: M —> R? closes up with genus 1, 
then 


W(f) = 27°. 


We will not prove this theorem. Part (i) of Theorem 13.6 was proved by Tom 
Willmore in [46] in 1965. The minimum is attained for a round sphere. Part (ii) was 
proved by Rob Kusner in [21] where he also proved that the Boy surface shown 
on the right of Fig. 13.1 realizes the minimum 127x. The two surfaces on the right 
of Fig. 13.2 are Lawson surfaces which were found by Blaine Lawson [23] and are 
possible candidates for minimizing the Willmore functional among all surfaces with 
genus g = 2 and g = 3 respectively [18]. 

In the paper already mentioned above, Willmore also formulated (iii) as a 
conjecture and demonstrated that the value 27° is realized by the torus obtained 
by rotating a circle of radius one around an axis in such a way that its center has 
distance /2 from the axis (Fig. 13.2, left). This Willmore conjecture remained 
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Fig. 13.1 A Boy surface with the minimal possible Willmore functional 127r 


P - ™, x a S 


Fig. 13.2 The Torus on the left has Willmore functional 277, which is optimal among surfaces 
with genus g = 1. The two surfaces on the right are possible candidates for minimizing the 
Willmore functional among all surfaces with genus g = 2 and g = 3 respectively 


a famous open problem in Differential Geometry for a long time, until in 2012 
Fernando Marques and André Neves proved the conjecture [26]. 


> Remark 13.7 The question of critical points of the Willmore functional acquired 
greater importance starting from the 1960s, initiated by T. Willmore and his paper 
[46]. It was later found that parts of the theory were already known to Wilhelm 
Blaschke [6] and his student Gerhard Thomsen in the 1920s [38]. For an historic 
overview of contributions which were made to the problem see the last chapter of 
[47], or [27] for a more recent survey. 
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13.2 Variation of the Willmore Functional 


According to the discussion in Sect. 13.1, the Willmore functional W has alternative 
versions which measure how “non-flat” or how “not round” a surface is. It was also 
explained that for the purposes of Variational Calculus all these different versions 
of the Willmore functional are equivalent. Being a critical point of the Willmore 
functional (which version we take does not matter) means that the surface (at least 
locally) is “optimally round”. It also means that the total amount of curvature of 
the surface cannot be decreased by modifying f only in a small neighborhood of a 
given point, while leaving the rest of the surface unchanged. 


Definition 13.8 


A surface f: M — R? is called a Willmore surface if it is a critical point of the 
Willmore functional W. 


Let us first compute for W the rate of change under a general variation, not 
necessarily with support in the interior: 


Theorem 13.9 (First Variation Formula for the Willmore Functional) 

Let f: M — R? be a surface with unit normal N and with binormal field 
B along the boundary 0M. Lett œ> f, be a variation of f with variational 
vector field 


f =ON +df(Z) 
where ¢ € C (M) and Z € F(T M). Then 


d 


dt 


w=- f o (AH + 2H(H? — K)) det 
M 


t=0 


+f (B, H?Z — H grado + ¢ grad H) ds. 
aM 


Proof. Using Theorem 12.18 as well as the notation G := grad @ borrowed from 
there we obtain 


: 1 
(n? det) = (21 (azu - 5div G - 62H? — K) + H?(QHọ + div z)) det 
= (2H (grad H, Z) — HdivG — 26H (H? — K) + H?div Zz) det 


= (diviH?z) — H div G — 26H(H? — K)) det 
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= (diveH?z — HG) + (grad H, grad $) — 26H (H? — K)) det 

= (aiv(H?z — HG) + div(¢ grad H) — PAH — 26H (H? — K)) det 

= (aiv (az — HG + ġ grad H) ang (AH +2H(H? — K))) det. 
Together with the Divergence Theorem 12.7, this proves our claim. o 


As an immediate consequence, we obtain (cf. [38]) 


Theorem 13.10 
A surface f: M —> R? is a Willmore surface if and only if 


AH +2H(H? — K)=0. 


Round spheres are Willmore, because for them all points are umbilic points (so 
H? — K = 0) and H is constant (so AH = 0). Moreover, all surfaces with H = 0 
(minimal surfaces) are Willmore. Here is another example: 


Example 13.11 


Take a unit speed curve y : [0, L] > IR? c R3, where R? is realized as those points in IR? where 
the last coordinate is zero. Now for a compact domain with smooth boundary M C [0, L] x R 
define the cylinder f: M — R? over y by 


fu, v) = y(u) + ves. 


It is easy to check that the Levi-Civita connection of f is given by VU = VV = 0, the Gaussian 
curvature K of f vanishes and the mean curvature H of f satisfies 


H(u, v) = w 
(grad H) (u, v) = £ w U (u, v) 
(AH)(u, v) = aa 


This means that the cylinder f over y is Willmore if and only if y is freely elastic, i.e. 


The cylinder over a freely elastic curve is seen in Fig. 13.3. 
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f 


A 


Pá 


\ 
\ OA 
A” A 


> 
J J 


Fig. 13.3 The cylinder over a free elastic plane curve is a Willmore surface 


Fig. 13.4 Another Willmore surface 


There are many other ways to construct Willmore surfaces, most of which are 
beyond the scope of this book. The surface in Fig. 13.4 is from the 2019 paper [7]. 
13.3 Willmore Functional Under Inversions 


For a surface f: M > R?, the Willmore functional 
W(f) = i H? det 
M 


is clearly unchanged if we postcompose f by an isometry g: R? > R?. It is also 
invariant under scaling: For 4 Æ 0 the surface f = Af has the same Willmore 
functional. This is because, under such a scaling, det acquires a factor of i while 
H gets a factor of L, As its name indicates, if we consider the Möbius-invariant 
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Fig. 13.5 The images under an inversion of the surfaces in Figs. 13.3 and 6.8 respectively are also 
Willmore surfaces 


Willmore functional 
win = f (H? — K) det 
M 


a similar statement is true for a more general class of transformations, that can be 
written as compositions of isometries, scalings and inversions in spheres, the so- 
called Mobius transformations (Fig. 13.5). 

Let f: M — R? bea surface such that the origin o of R? is not in the image of 
jf. Then we can postcompose f with the so-called inversion in the unit sphere 


: R? R3, g) = >— 
8: R? \ {0} > 8 (p) atl 


and obtain a new surface 


: | a 
: R?, = : 
PMR T 


Computing the derivative of f is straightforward and yields 


z_ af (df, f) f 
df= =2 
d (F f) (F SY 
1 


=n T 
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where for each p € M the orthogonal (3 x 3)-matrix R(p) € O(3) acts onv € R? 
as 
(f (p), V) 
R(pv =v -2—1 f(p). 
(F(p), f(P)) 

For each p € M the matrix R(p) is a reflection and hence orientation-reversing. 
The sign of the unit normal depends on orientation, which is why the unit normal 
field of f is given by 


N = -RN 


Il 
as 
z 


Theorem 13.12 
In the situation above, the induced metric (,)~, the area form det and the 
shape operator A of f are given by 


na 1 
(a) Tp pe! 
aS = ae 
GEDA 


A = -(f, fyA+2(N, fYI. 


Proof. The first two formulas follow directly from our calculations above. The third 
follows from 


- (XAN, f) ain, dh) 2. fae 
( Ff) ar J) Gp 


=2(N, fidf —(f, fidf oA. 


Theorem 13.13 
If f arises from f by inversion in the unit sphere, then 


Wf) =W). 
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Proof. By Theorem 13.12, the principal curvatures K1, K2 of f satisfy 
Ka — Kı = —(f, f)(k2 — K1). 


AS a consequence, 


a poe le eee a _ py 
(H Kdet= zl K1) det = ee kı) det = (H K)det. 


Theorem 13.14 
If f is a Willmore surface, then so is its image f under inversion in the unit 
sphere. 


Proof. By Theorem 13.4, W has the same critical points as W and by Theo- 
rem 13.13 inversion in the unit sphere maps critical points of W to critical points of 
W. o 


The surface on the right of Fig. 13.6 shows the image under an inversion of a 
minimal surface already known to Euler (shown on the left), the so-called catenoid 
f: M > R? given by 


1+u?+v? 


u 
utu 


2 
fv) =] mény 


log(u* + v?) 


So, even more Willmore surfaces can be obtained by inverting surfaces which we 
have already encountered (see Fig. 13.5). 


ZO 


Fig. 13.6 A catenoid (left) and its image under a sphere inversion (right) 
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Some Technicalities 


A.1 Smooth Maps 


The standard definition of a differentiable map f: U — R” requires U to be an 
open subset of IR. Such an f is called smooth if all higher order partial derivatives 


a” fi 
Xj- OX jy 
of all its component functions exist. 
On the other hand, we want to define curves in R” as certain smooth maps 
y: [a,b] — R” defined on a closed interval [a,b] C R. Similarly, we want to 


use a certain kind of compact subsets M C R? as the domain of definition for 
surfaces f: M — R”. We therefore have to work with more general domains: 


Definition A.1 


Let U C RE be an open set and M := U its closure. Then a function f: M —> 
R” is called smooth if there is an open set U C R” with M C U and a smooth 
function f: U —> R” such that 


fim = f. 
For x € M we define 


fo PÀ 
DX jy... OX jm i DX jy OX jm 
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In order for this definition to make sense, we have to verify that the higher partial 
derivatives of f are well-defined: 


Theorem A.2 
The higher partial derivatives of f defined in Definition A.1 are independent 
of the choice U and the extension f of f to U. 


Proof. Because every point x € M is a limit point of points y € U and the partial 
derivatives of f are continuous, we have 


qm rA : qm 
ia a a ae) = a a — 0) 
Xj, +. OX jm ae Xj, OX jm 
g” 
= lim l w) 
y> OX j, OX jn 


oO 


When we discuss reparametrizations of curves or surfaces, we make use of the 
following notion: 
Definition A.3 


Let M, M c RÉ be two subsets which are closures of open subsets U, Ü c R” 
respectively. Then a smooth map f: M — M is called a diffeomorphism if it 
is bijective and its inverse f~!: M —> M is also smooth. 


A.2 Function Toolbox 
On several occasions in this book the need arises to construct a so-called bump 


function, i.e. a non-negative smooth function on R that vanishes outside of a small 
neighborhood of a given point, but not at this point. 


Definition A.4 


Let f: A — R be a function defined on some subset A C R* of R*. Then the 
support of f is defined as 


supp f := {x € R” | Ye > 0 Jy €e A with |y — x| < € and f(y) 4 0}. 


The following theorem is easy to prove. 
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Theorem A.5 
If f : U + Ris a smooth function on an open set U C RÝ and supp f C U, 
then f can be extended to a smooth function f : R* — R by setting 


f(x) forxeU 


ez fF forx € U. 


The basic ingredient for constructing functions with support in a given open set 
U C R* is the function f: R —> R given by 


0 for x <0 
= 1 


for x > 0. 


=] 


e 


Clearly, f is smooth at all points x Æ 0. It is not hard to check that it is smooth also 
at x = 0. Figure A.1 shows the graph of f. 

The second function in our toolbox is the so-called bump function g: R —> R 
given by 


2 0 for |x| > 1 
gœ) = fA- _ 1 


e i=? forjx|< 1. 


As a composition of smooth functions, g is also smooth (see Fig. A.2). 
Other versions of g like 


&(x) = g (ex — xo)) 


Fig. A.1 A smooth fuction f : R —> R with f(x) = 0 for x <0 
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Fig. A.2 The bump function g: R > R 


Fig. A.3 A bump function g on R? 


can be adapted to be non-zero only within an arbitrarily prescribed interval. Another 
tool in our toolbox is the function h: R — R given by 


ney = f g 
-1 


and variants of it that are shifted and scaled in a similar way as the function g above. 
Finally, we need bump functions of several variables (see Fig. A.3), like the 
function : R —> R given by 


a(x) = fA- Ixl®. 


Timeline 


Table B.1 is not meant as a comprehensive overview of the history of the whole 
field that deals with the Differential Geometry of curves and surfaces. Only those 
milestones are listed that are explicitly mentioned in the preceding chapters. 


Table B.1 The milestones of 
differential geometry of 
curves and surfaces which are 


Year Milestone Section 


1673 Newton defines the curvature of curves in R? 3.1 


covered in this book 1691 Jacob Bernoulli defines elastic curves in R? 2.4 
1744 Euler classifies elastic curves in R? 25 
1744 Euler shows that the catenoid minimizes area 133 
1760 Euler defines the principal curvatures of a 8.2 
surface 

1827 Gauss proves the Theorema Egregium 93 

1844 Binet derives the equation of elastic curvesin 5.4 
R3 

1845 Möbius investigates the topology of closed ikl 
surfaces 

1848 Bonnet proves the Gauss-Bonnet theorem 10.2 


1859 Kirchhoff proves that the tangent of an elastic 5.2 
curve follows the motion of the axis of a 
spinning top 


1903 Boy proves the Gauss-Bonnet theorem for 13 
closed surfaces 
1906 Da Rios defines the filament equation 33 


1923 Thomsen defines Willmore surfaces, then called 13.1 
Konformminimalflächen 


1931 Douglas and Rado independently prove the 12.4 
existence of a minimal surface with prescribed 
boundary curve 


(continued) 
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Table B.1 (continued) 


B Timeline 
Year Milestone Section 
1937 Whitney and Graustein proof their theorem 3.6 
1956 Hopf proves that round spheres are the only PS 
constant mean curvature surfaces in R? of 
genus zero 
1965 Willmore states his conjecture BA 
1970 Lawson finds closed minimal surfaces in S? tsil 
with any genus 
1972 Hasimoto shows that the filament equation isa 5.3 
soliton equation 
1984 Wente finds the first constant mean curvature 12.5 
torus 
2012 Marques and Neves prove the Willmore ISl 


conjecture 
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Symbols Connected compact domain with smooth 
90°-rotation, 29 boundary, 88 

1-form, 109 Covered area form, 125 

2-form, 106 Critical point, 20, 168 


constraint, 22 
under constraints, 168 


A Curvature, 30 
Absolute curvature, 50 function, 53 
Angular velocity, 134 line, 140 
Arclength, 8 total curvature, 38 
Arclength 1-form, 112 Curve, 3 
Arclength coordinate, 8 closed, 4 
Arclentgh parametrization, 9 framed, 73 
Area covered by a map, 126 on a surface, 139 
Area form, 98 oriented, 5 
Area of a Riemannian domain, 108 parallel curve, 49 
Asymptotic line, 140 Cylinder, 186 

Cylinder over a free elastic plane curve, 187 
B 
Bending energy, 10, 182 D 
Binormal, 140 Derivative, 91 
Boy’s Theorem, 156 Derivative with respect to arclength, 8 
Boy surface, 154, 155 Determinant, 96 
Bump function, 195 Diffeomorphism, 89, 194 

orientation-preserving, 89 

Direction, 120 
C Directional curvature, 120 
Cartesian leaf, 5 Directional derivative, 92 
Catenoid, 190, 191 Divergence, 162 
Cauchy-Schwarz inequality, 50 Divergence Theorem, 164 


Codazzi equation, 135 
Compact domain with smooth boundary, 88 


Component, 159 E 

Cone volume, 174 Elastic curve, 65 

Congruent, 7 free elastic curve, 21 

Connectable, 159 torsion-free elastic curve, 25 

Connected, 159 Elastic framed curve with twisting modulus, 77 
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Endomorphism field, 95 Lawson surfaces, 183 
Enneper surface, 102, 174 Length, 6 
Euclidean plane, 125 Levi-Civita connection, 132 


Euler characteristic, 142 


M 

F Möbius transformations, 188 
First Variation Formula of Area, 172 Mean curvature, 122 
First Variation Formula of the Willmore Metric, 94 

Functional, 185 Minimal, 173 
First Variation of Cone Volume, 176 Minimal surface, 173 
Foucault’s pendulum, 147 Monodromy angle, 148 
Framed curve, 73 
Frenet normal, 80 
Fundamental Lemma of the Calculus of N 

Variations, 20 Neil’s parabola, 5 
Fundamental Theorem of Plane Curves, 31 Norm, 95 


Normal curvature, 140 
Normal field, 47 


G Normal part, 131 
Gauss-Bonnet theorem, 142 Normal space, 47 
Gauss equation, 136 Normal transport, 59 


Gaussian curvature, 122 
Genus, 158, 159 


Geodesic, 140 O 

Geodesic curvature, 140 Optimization under constraints, 22 
Geodesic torsion, 140 Orientation cover, 157, 158 
Gradient, 23, 162 Orientation-preserving, 6 


Gravity vector, 26 


P 
H Parallel, 48 
Hasimoto curvature, 51 endomorphism field, 133 
Helix, 5 normal field, 48 
Homotopic, 39 transport, 49 
through closed curves, 42 transport map, 146 


Parametrized by arclength, 9 
Pendulum equation, 26 


I Plateau problem, 173 
Infinitesimal rigid motion, 71 Poisson ratio, 76 
Integral, 107 Positively oriented basis, 97 
of a l-form, 110 Principal curvatures, 121 
of a 2-form, 107 Projection map, 91 
over arclength, 8 Pseudosphere, 101 
Isometric, 100 Pull-back, 106, 109 
K R 
Klein bottle, 157 Regular homotopy, 39 


Regular homotopy through closed curves, 42 
Reparametrization, 6, 90 

L Riemannian domain, 95 

Lagrange multiplier, 23, 169 Riemannian metric, 94 

Laplacian, 163 
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Rigid motion, 7 
Round sphere, 124 


S 
Schwarz-P, 173 
Sector area, 33 
Shape operator, 119 
Shear modulus, 76 
Smoke ring flow, 72 
Smooth, 3, 193 
functional, 19, 167 
map, 3 
one-parameter family of maps, 13 
Speed, 6 
Spherical pendulum, 26 
Support, 194 
Support in the interior, 19 
Surface, 88 
closed surface, 153 
closes up, 155 
normal, 140 
oriented closed surface, 154 


T 

Tangent bundle, 91 
Tangential part, 131 

Tangent space, 91 

Tangent vector, 91 

Tangent winding number, 38 
Tension, 25 

Theorema Egregium, 136 
Time, 165 

Torsion, 74 

Total geodesic curvature, 141 
Total squared curvature, 36 
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Total squared torsion, 76 
Total torsion, 61 

Total torsion angle, 74 
Trace, 96, 139 

Twisting modulus, 76 


U 

Umbilic point, 123 
Unit normal, 117 
Unit tangent field, 9 


v 
Variation, 14 

of a smooth map, 14, 166 

of a surface, 167 

support, 77 

support in the interior, 168 
Variational vector field, 14, 167 
Vector field, 91 

coordinate vector fields, 93 
Vortex filament, 69 
Vortex filament flow, 72 
Vorticity, 69 


WwW 

Wente torus, 178 

Willmore conjecture, 183 

Willmore functional, 181 
conformally invariant, 182 

Willmore surface, 185 


Y 
Young modulus, 76 


